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The  most  productive  approaches  to  the  problem  of  the  gamma- ray  laser  have 
focused  upon  upconversion  techniques  in  which  metastable  nuclei  are  pumped 
with  long  wavelength  radiation.  At  the  nuclear  level  the  storage  of  energy 
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broad  interdisciplinary  nature  requiring  the  fusion  of  concepts  taken  from 
relatively  unrelated  fields  of  physics.  Our  research  group  has  described 
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20.  Abstract  (continued) 

several  means  through  which  this  energy  might  be  coupled  to  radiation  fields 
with  cross  sections  for  stimulated  emission  that  could  reach  10-17  cm2.  Such 
a  stimulated  release  could  lead  to  output  powers  as  great  as  3  x  1021 
Watts/liter.  Since  1978  we  have  pursued  an  approach  for  the  upconversion  of 
longer  wavelength  radiation  incident  upon  isomeric  nuclear  populations  that 
can  avoid  many  of  the  difficulties  encountered  with  traditional  concepts  of 
single-photon  pumping.  Experiments  have  confirmed  the  general  theory  and 
have  indicated  that  a  gamma-ray  laser  is  feasible  if  the  right  combination 
of  energy  levels  and  branching  ratios  exists  in  some  real  material.  Of  the 
1,886  distinguishable  nuclear  materials,  the  present  state-of-the-art  has 
been  adequate  to  identify  29  first-class  candidates,  but  further  evaluation 
cannot  proceed  without  remeasurements  of  nuclear  properties  with  higher 
precision.  A  laser-grade  database  of  nuclear  properties  does  not  yet  exist, 
but  the  techniques  for  constructing  one  have  been  developed  and  utilized 
under  this  contract.  Resolution  of  the  question  of  the  feasibility  of  a 
gamma-ray  laser  now  rests  upon  the  determination  of:  1)  the  identity  of  the 
best  candidate,  2)  the  threshold  level  of  laser  output,  and  3)  the  upconver¬ 
sion  driver  for  that  material. 

Previous  work  focused  upon  our  approach  that  is  the  nuclear  analog  to  the 
ruby  laser.  It  embodies  the  simplest  concepts  for  a  gamma-ray  laser  and  not 
surprisingly,  the  greatest  rate  of  achievement  in  the  quest  for  a  subAng- 
strom  laser  was  realized  in  that  direction.  For  ruby  the  identification  and 
exploitation  of  a  bandwidth  funnel  were  the  critical  keys  in  the  development 
of  the  first  laser.  There  was  a  broad  absorption  band  linked  through 
efficient  cascading  to  the  narrow  laser  level. 

In  1987  we  reported  a  major  milestone  which  showed  that  comparable  structure 
existed  at  the  nuclear  scale  in  the  first  of  the  29  candidate  isomers 
available  for  testing,  180Tani.  Populations  of  the  isomer  were  successfully 
pumped  down  with  flashes  of  x  rays  absorbed  through  an  astonishingly  large 
cross  section  of  40,000  on  the  usual  scale  (~  10~28  cm2  keV)  where  10 
describes  a  fully-allowed  process.  This  corresponded  to  a  partial  width  for 
useful  absorption  of  0.5  eV,  even  better  than  what  had  been  assumed  for 
idealized  nuclei.  In  the  course  of  successive  work  we  discovered  that  the 
giant  pumping  resonances  occurred  with  a  gratifying  frequency  throughout  the 
table  of  nuclides,  reaching  optimal  size  and  strength  in  the  mass  region 
where  the  better  candidates  lie. 

At  the  conclusion  of  that  work  we  focused  upon  the  extension  of  our  achieve¬ 
ments  in  pumping  nuclei  to  those  having  laser- like  transitions  with  micro¬ 
second  lifetimes.  For  the  test  cases  of  181Ta  and  176Hf  giant  resonances  were 
found  for  the  pumping  of  fluorescence  with  the  Texas -X  linear  accelerator 
installed  in  our  facility  for  dedicated  use  in  gamma-ray- laser  research. 
The  importance  of  those  first  results  is  that  the  particular  nuclei  studied 
are  very  close  neighbors  in  the  sense  of  nuclear  structure  to  the  still 
unavailable  first-class  candidate  nuclei  for  a  gamma  ray  laser.  Such 
favorable  results  with  those  close  simulations  argued  strongly  for  analogous 
successes  with  the  actual  candidate  nuclei. 

During  the  most  recent  reporting  period  efforts  focused  upon  identifi¬ 
cation  and  study  of  the  best  of  the  29  candidate  isomers  for  a  gamma-ray 
laser,  the  31  year  isomer  of  hafnium,  178Hfm2.  It  stores  the  greatest  energy 
density  available  for  release  in  a  gamma- ray  laser,  1.3  kilojoules  per 
microgram.  Because  of  the  expected  difficulty  in  obtaining  samples  of  this 
exceedingly  rare  material,  efforts  this  past  year  were  concentrated  upon 
systematic  studies  and  new  technologies  that  will  reduce  the  amount  of 
isomeric  material  needed  for  the  critical  experiments  planned. 
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Status  and  issuas  in  tha  development 
of  a  Y"r*y  luu; 

ZZI .  Preferred  nudidas 

By  C.  B.  Collins  and  J.  J.  Carroll 

University  of  Texas  at  Dallas, 

Center  for  Quantum  Electronics 
P.0.  Box  830688,  Richardson,  TX  75083-0688 

Abstract 

A  gamma-ray  laser  would  stimulate  the  emission  of  radiation  at 
wavelengths  below  1  A  from  excited  states  of  nuclei.  However,  the 
anticipation  of  a  need  for  very  high  pump  powers  tended  to 
discourage  early  research  and  the  difficulties  in  demonstrating  a 
device  were  first  assumed  to  be  insurmountable.  Over  the  past 
decade  advances  in  pulsed  power  technology  have  changed  these 
perceptions  and  studies  have  built  a  strong  momentum.  The  concept 
for  a  nuclear  analog  of  the  ruby  laser  now  appears  in  the  textbooks 
and  many  of  the  component  steps  for  pumping  the  nuclei  have  been 
demonstrated  experimentally.  A  quantitative  model  based  upon  the 
new  data  and  concepts  has  shown  the  gamma-ray  laser  to  be  feasible 
if  some  real  isotope  has  its  properties  sufficiently  close  to  the 
ideals  modeled.  The  greatest  positive  impact  upon  feasibility  has 
come  from  the  discovery  of  giant  resonances  for  pumping  nuclei  that 
greatly  reduce  the  levels  of  pump  power  needed.  Most  recently 
attention  has  been  focused  upon  efforts  to  demonstrate  prelasing 
levels  of  fluorescence  from  nuclides  which  are  actual  candidates 
for  use  in  a  gamma-ray  laser.  Problems  being  addressed  are  the 
acquisition  of  macroscopic  samples  for  testing  and  the 
demonstration  of  appropriate  instrumentation. 


1 


1 .  Introduction 

A  decade  ago  the  difficulties  in  realizing  a  gamma-ray  laser 
were  considered  almost  insurmountable.  However,  in  1982  we  began 
to  emphasize  an  interdisciplinary  concept  of  upconversion  which 
ultimately  launched  a  renaissance  in  the  field  (Collins,  et  al. 
1982).  The  essential  concept  was  the  "optical"  pumping  of  nuclei. 
In  this  case  optical  meant  x  rays,  but  the  fundamentals  were  the 
same.  Useful,  resonant  absorption  of  pump  power  would  occur  over 
short  distances  to  produce  high  concentrations  of  excited  nuclei 
while  wasted  wavelengths  would  be  degraded  to  heat  in  much  large 
volumes.  Of  all  the  cases  considered,  the  nuclear  analog  of  the 
ruby  laser  embodied  the  simplest  concepts  for  a  gamma-ray  laser. 
Not  surprisingly,  the  greatest  rate  of  achievement  has  been 
realized  in  that  direction. 

For  ruby,  the  identification  and  exploitation  of  a  bandwidth 
funnel  were  the  critical  keys  in  the  development  of  the  first 
laser.  There  was  a  broad  absorption  band  linked  through  efficient 
cascading  to  the  narrow  laser  level.  Our  concept  called  for  a 
nuclear  analogy  of  this  structure  which  was  unknown  in  1986  when 
the  first  phase  of  experiments  was  started.  Moreover,  we  had 
envisioned  pumping  not  only  the  ground  states  of  nuclei,  but  also 
isomeric  levels.  Such  metastable  nuclei  store  keV  and  even  MeV  of 
energy  for  years  so  most  of  the  pump  power  could  be  input  long 
before  the  time  of  use  and  triggering  requirements  would  be 
minimized.  The  model  has  been  confirmed  in  principle  and  now 
appears  in  the  textbooks  (Collins  1991) . 


The  use  of  intense  flashes  of  x  rays  to  dump  the  stored  energy 
of  an  isomer  was  first  shown  (Collins,  et  al.  1988}  for  the  rare 
nuclide  l*°Ta"  which  exists  naturally  as  100%  inverted.  Then 
experiments  with  the  unique  superconducting  linac  at  Darmstadt 
showed  that  populations  of  the  isomer  1,#Ta“  are  pumped  down  through 
two  isolated  resonances  of  extraordinary  strength,  about  10*  times 
larger  than  expected  for  the  excitation  of  nuclei  by  x  rays 
(Collins,  et  al.  1990).  A  quantitative  model  based  upon  the  new 
data  showed  the  gamma-ray  laser  to  be  feasible  if  one  of  the  real 
nuclear  isomers  has  its  properties  sufficiently  close  to  the  ideals 
modeled  (Collins,  et  al.  1992;  Collins,  et  al.  1993). 

The  greatest  impediment  to  research  on  the  feasibility  of  a 
gamma-ray  laser  has  come  from  the  lack  of  suitable  materials. 
There  are  1886  distinguishable  nuclear  isotopes,  but  most  do  not 
occur  naturally  and  there  have  not  been  any  reasons  to  produce  them 
artificially.  In  1986  a  search  of  all  nuclear  properties 
identified  29  isotopes  as  first-class  candidates  for  a  gamma-ray 
laser  and  three  are  shown  in  Table  I  as  examples.  Only  the  two 
poorest,  180Tam  and  wsTe  could  be  acquired  for  testing.  Those  are  the 
ones  upon  which  the  breakthrough  discovery  of  the  giant  pumping 
resonance  was  made. 

The  better  candidates  for  a  gamma-ray  laser  have  never  been 
available  in  sufficient  quantities  for  evaluation  and  research  had 
to  focus  instead  upon  19  simulation  materials  that  served  to  teach 
the  systematics  for  using  the  giant  pumping  resonances  to  excite 
nuclei  (Carroll,  et  al.  1991).  The  results  on  those  simulated 
candidates  were  extremely  encouraging  as  the  giant  resonance  for 


pumping  with  bremsstrahlung  were  found  in  many  of  them.  Moreover, 
the  results  have  considerable  significance  as  they  served  to 
identify  the  best  candidate  material  for  pumping  a  gamma-ray  laser, 
the  31-year  isomer  of  l7,Hf. 

2.  Technical  Background 

2 . 1  Pump  Model 

The  sequence  for  triggering  the  release  of  the  energy  stored 
in  an  isomeric  state  (Collins,  et.  al.  1982;  Collins  1991)  is 
described  by  an  integrated  cross  section  for  the  pumping  of  that 
level  with  a  continuum  of  x  rays.  The  population  is  transferred 
through  a  gateway  state  where  the  selection  rules  that  would 
otherwise  limit  the  process  do  not  apply.  The  normal  decay  from 
the  gateway  is  accompanied  by  the  emission  of  immediate 
fluorescence  that  leads  to  the  principal  laser  level  from  which  the 
sustained  output  of  power  will  be  emitted.  Such  processes  fall 
into  the  general  category  of  (Y/Y')  reactions  shown  schematically 
in  Fig.  1.  Population  can  be  accumulated  in  the  output  level  by 
continuing  to  run  the  pump  cycle  for  a  time  comparable  to  the 
lifetime  of  the  output  state.  In  the  case  of  the  17,Hf  that  lifetime 
is  an  almost  ideal  value  of  about  70  sec. 

Whether  or  not  the  initial  state  being  pumped  in  a  iy,y’) 
reaction  is  the  ground  state  or  a  long-lived  isomeric  level,  the 
principal  figure  of  merit  is  the  integrated  cross  section  for 
pumping  an  output  transition  through  such  a  broad  band,  Jib.boroo/2. 
Constituent  parameters  are  identified  in  Fig.  1  where  it  can  be 
seen  that  the  branching  ratios  b«  and  b„  specify  the  probabilities 
that  a  population  pumped  by  absorption  into  the  j-th  broad  level 


will  decay  back  into  the  initial  state  or  by  cascade  to  the  output 
level,  respectively. 

For  a  sample  which  is  optically  thin  at  the  pump  wavelength, 
a  computation  of  the  number  of  nuclei  pumped  into  a  fluorescence 
level  in  the  scheme  of  Fig.  1  is  straightforward  (Collins,  e£  al. 
1982;  Collins  1991) .  Most  intense  x-ray  sources  emit  continua, 
either  because  bremsstrahlung  is  initially  produced  or  because 
spectral  lines  are  degraded  by  Compton  scattering  in  the  immediate 
environment.  The  time-integrated  yield  of  final-state  nuclei,  N, 
obtained  by  irradiating  initial  targets  with  an  x-ray  flux,  in 
photons  cm*2  delivered  in  a  continuum  is, 

Nf  =  N^0  f  0  (E)  F(E,  E0)  dE  ,  (1) 

0 

where  F(E,E0)  is  the  distribution  of  intensities  up  to  an  endpoint 
energy  E0  within  the  input  spectrum  normalized  so  that 

So 

f  F(E,E0)dE  =  1  ,  (2) 

o 

and  a(E)  is  the  effective  cross  section  for  the  excitation  of  the 
final  state  from  the  initial. 

All  (Y/Y')  reactions  occurring  at  energies  below  the  threshold 
for  particle  evaporation  excite  discrete  pump  bands,  or  gateways, 
as  shown  in  Fig.  1.  Although  only  one  gateway  appears  in  the 
figure,  there  could  be  more.  Each  would  be  excited  at  a  different 
pump  energy  but  all  would  branch  to  some  extent  into  the  same 


fluorescence  level,  f.  The  j-th  gateway  is  shown  in  Fig.  1  as 
typical. 

Although  the  width  of  level  j  is  broad  on  a  nuclear  scale,  it 
is  narrow  in  comparison  to  the  scale  of  energies,  E  over  which 
F(E,E0)  varies.  Then,  the  final-state  yield,  expressed  as  the 
normalized  activation  per  unit  photon  flux,  A^Eo)  can  be  written 
from  Eq.  (1)  as, 


l£t0  =  £  <°r)«  *  (3) 

In  this  expression  (or)fj  is  the  integrated  cross  section  for  the 
production  of  final-state  N,  as  a  result  of  the  excitation  of  the 
gateway  at  Ej  ,  so  that 

(aT)fj  =  f  a  (E)  dE  ,  (4) 

JTj-A 

where  A  is  an  energy  small  compared  to  the  spacing  between  gateways 
and  large  in  comparison  to  their  widths.  Levels  of  this  type  are 
sometimes  called  intermediate  states  or  doorways. 

It  is  straightforward  to  show  that, 

(or)fJ  =  (nbabJ'aQ/2)  fj  ,  (5) 


where  o0/2  is  the  peak  of  the  Breit-Wigner  cross  section  for  the 
absorption  step,  and 


°0 


XI  2I,+1  i 
2n  2Ig+l  ap*l 


(6) 


where  X  is  the  wavelength  of  the  gamma  ray  at  the  resonant  energy, 
E.,;  I.  and  Ig  are  the  nuclear  $pips  of  .the  excited  and  ground  states, 
respectively;  and  ap  is  the  total  internal  conversion  coefficient 
for  the  system  shown  in  Fig.  1. 

2.2  Experimental  Confirmation 

In  experimental  work  of  the  last  five  years  the  bremsstrahlung 
from  five  accelerators  in  different  experimental  environments  was 
used  to  verify  the  model  of  Eqs.  (1)  -  (4)  for  pumping  nuclear 
fluorescence  and  to  cross-check  the  accelerator  intensities 
(Carroll,  et  al.  1991) .  Spectral  intensities  were  calculated  with 
the  EGS4  coupled  electron/photon  transport  code  (Nelson,  et  al. 
1985)  adapted  for  each  individual  configuration  from  closely 
monitored  values  of  accelerator  currents.  In  this  way  both  F(E,E0) 
and  $0  were  obtained.  In  some  cases  ♦,  was  verified  by  in-line 
dosimetry. 

During  these  experiments  samples  with  typical  masses  of  grams 
were  exposed  to  the  bremsstrahlung  from  the  five  accelerators  for 
times  ranging  from  seconds  to  hours  for  the  continuously  operating 
machines  and  to  single  flashes  from  the  pulsed  devices.  Results 
were  in  close  agreement  (Anderson  &  Collins  1987;  Anderson,  et  al. 
1988a;  Anderson  &  Collins,  1988b)  with  the  predictions  of  Eq.  (3) 
used  with  literature  values  of  (or)tj.  The  most  detailed 
confirmation  of  theory  was  obtained  with  the  reaction  *7Sr  ( y ,  Y  *  ),7Sr" 


as  shown  in  Fig.  2.  It  served  to  completely  confirm  the  model  and 


to  validate  the  us*  of  the  EGS4  code  for  calculating  bremsstrahlung 
intensities  from  measured  accelerator  parameters. 

Many  of  the  29  candidates  belong  to  the  class  of  nuclei 
deformed  from  the  normally  spherical  shape  (Collins  1991) .  For 
those  systems  there  is  a  quantum  number  of  dominant  importance,  K 
which  is  the  projection  of  individual  nucleonic  angular  momenta 
upon  the  axis  of  elongation.  To  this  is  added  the  collective 
rotation  of  the  nucleus  to  obtain  the  total  angular  momentum  J. 
The  resulting  system  of  energy  levels  resembles  that  of  a  diatomic 
molecule. 

In  most  cases  an  isomeric  state  has  a  large  lifetime  because 
its  value  of  K  differs  considerably  from  those  of  lower  levels  to 
which  it  would  otherwise  be  radiatively  connected.  As  a 
consequence,  bandwidth  funneling  processes  such  as  shown  in  Fig.  1 
that  start  from  isomeric  levels  must  span  substantial  changes  in  K 
and  component  transitions  have  been  expected  to  have  large,  and 
hence  unlikely,  multipolarities. 

From  this  perspective  the  third  of  the  isomers  listed  in  Table 
I,  1,0Ta  was  initially  unattractive  as  it  had  one  of  the  largest 
changes  of  angular  momentum  between  isomer  and  ground  state,  8K. 
However,  because  a  macroscopic  sample  was  readily  available,  1,0Ta 
became  the  first  isomeric  material  to  be  optically  pumped  to  a 
fluorescent  level. 

This  particular  isomer,  “"Ta*  carries  a  dual  distinction.  It 
is  the  rarest  stable  isotope  occurring  in  nature  and  it  is  the  only 
naturally  occurring  isomer.  The  actual  ground  state  of  1,0Ta  is  1* 
with  a  halflife  of  8.1  hours  while  the  tantalum  nucleus  of  mass  180 


occurring  with  0.012%  natural  abundance  is  the  9_  isomer,  “°Ta".  It 
has  an  adopted  excitation  energy  of  75.3  keV  and  a  half  life  in 
excess  of  1.2  x  10“  years. 

In  an  experiment  conducted  in  1987,  1.2  mg  of  “°Ta"  was  exposed 
to  the  bremsstrahlung  from  a  6  MeV  linac  and  a  large  fluorescence 
yield  was  obtained  (Collins,  et  al.  1988)  .  This  was  the  first  time 
a  (Y/Y')  reaction  had  been  excited  from  an  isomeric  target  as 
needed  for  a  gamma-ray  laser  and  was  the  first  evidence  of  the 
existence  of  giant  pumping  resonances.  Simply  the  observation  of 
fluorescence  from  a  milligram  sized  target  proved  that  an 
unexpected  reaction  channel  had  opened.  Usually  grams  of  material 
have  been  required  in  this  type  of  experiment.  Analyses  (Collins, 
et  al.  1988;  Collins,  et  al.  1990)  of  the  data  indicated  that  the 
partial  width,  b.bj’  for  the  dumping  of  “°Ta“  was  around  0.5  eV. 

To  determine  t^he  transition  energy,  Ej  from  the  “°TaB  isomer  to 

i 

the  gateway  level,  a  series  of  irradiations  (Collins,  et  al.  1990) 
was  made  at  the  S-DALINAC  facility  using  fourteen  different 
endpoints  in  the  range  from  2.0  to  6.0  MeV.  The  existence  of  an 
activation  edge  was  clearly  seen  in  the  data  shown  in  Fig.  3.  The 
fitting  of  such  data  to  the  expression  of  Eq.  (3)  by  adjusting 
trial  values  of  (or)fj  provided  the  integrated  cross  sections  for 
the  dumping  of  “°TaB  isomeric  populations  into  freely-radiating 
states.  Reported  values  (Collins,  et  al.  1990)  were  12,000  and 
35,000  in  the  usual  units  of  10'29  cm2  keV  for  gateways  at  2.8  and 
3.6  MeV,  respectively.  These  are  enormous  values  exceeding 
anything  previously  reported  for  transfer  through  a  gateway  by  two 


orders  of  magnitude.  In  fact  they  are  10,000  times  larger  than  the 
values  usually  measured  for  nuclei. 

A  survey  of  19  isotopes  (Carroll,  et  al.  1991) conducted  with 
the  four  U.S.  accelerators  over  a  fairly  coarse  mesh  of 
bremsstrahlung  endpoints  confirmed  the  existence  of  giant 
resonances  for  transferring  K  in  the  region  of  masses  near  180  as 
shown  in  Fig.  4.  Activation  edges  continued  to  support  the 
identifications  of  integrated  cross  sections  for  pumping  and 
dumping  of  isomers  that  were  of  the  order  of  10,000  times  greater 
than  usual  values.  The  details  of  these  studies  of  systematics 
have  been  reported  in  recent  reviews  (Collins,  et  al.  1992; 
Collins,  et  al.  1993)  together  with  the  analyses  of  the  nuclear 
structures  to  which  these  giant  pumping  resonances  may  be 
attributed.  The  conclusion  of  these  modeling  studies  is  that  the 
best  of  the  29  candidate  materials  for  a  gamma-ray  laser  is  the 
first  listed  in  Table  I.  Because  it  should  be  the  easiest  to 
trigger,  the  31-year  isomer  17,Hf**  is  the  nuclide  of  choice  for  the 
pumping  of  a  gamma-ray  laser. 

3.  Candidate  Isomer  for  a  Gamma-Bay  Laser 

3.1  Candidates  and  Simulations 

While  an  exact  ranking  of  the  29  candidate  isomers  for  a 
gamma-ray  laser  depends  upon  a  complex  weighting  of  combinations  of 
nuclear  parameters,  the  potential  importance  of  several  nuclides  is 
magnified  by  some  very  pragmatic  considerations.  These  issues  are 
typified  by  comparisons  of  the  three  examples  presented  in  Table  I. 

From  the  perspective  of  shelf  life  and  availability,  the  third 
entry  in  Table  I,  1,0Ta"  is  far  superior.  It  is  a  naturally 


occurring  material  composing  0.01%  of  all  Ta  and  can  be  prepared 
simply  by  separating  natural  Ta  by  atomic  masses.  Samples  of 
milligram  weight  exist  and  one  such  specimen  was  used  in  the 
breakthrough  experiment  proving  the  use  of  x  rays  to  dump  the 
energy  stored  in  isomeric  populations.  In  contrast  the  entire 
world  inventory  of  the  first  entry,  "•Hf1*2  appears  to  be  about  1015 
nuclei  (Oganessian,  et  al.  1992). 

However,  from  the  perspective  of  triggering,  1,0TaB  appears  to 
be  much  less  practical.  Although  the  energy  storage  is  still 
impressively  high,  the  isomeric  level  lies  quite  low  in  energy  when 
compared  to  the  value  for  the  K-breaking  gateway  at  2.8  MeV.  For 
178Hf' «2  the  energy  of  the  isomer  at  2.45  MeV  is  the  highest  known  that 
still  lies  below  the  likely  gateway  energy  between  2.5  -  3.0  MeV. 
The  isomer  17,Hf  is  even  higher  at  3.312  MeV,  but  the  availability 
(Walker,  et  al.  1990)  of  a  spontaneous  transition  down  to  the  K- 
mixing  level  at  2.685  MeV  reduces  the  shelf  half life  to  only  3.7 
ysec. 

The  ideal  energy  for  an  isomer  to  store  for  ease  in  triggering 
would  be  as  much  as  possible  without  exceeding  the  energy  of  the  K- 
breaking  level  for  that  nuclide.  As  long  as  the  transition  energy 
needed  for  triggering  is  positive,  the  isomer  cannot  dump 
spontaneously  as  happens  with  17,Hfm,  so  the  lifetime  of  the  initial 
population  will  be  as  long  as  possible  and  problems  of  storage  will 
be  minimized.  Guided  to  17,Hf,“2  by  energetics,  a  reasonable  next 
concern  from  the  pragmatic  viewpoint  is  for  the  duration  of  the 
trigger  pulse. 


Even  in  the  best  cases  the  requirement  for  the  energy  in  a 
trigger  pulse  is  large.  The  problem  is  compounded  if  the 
fluorescence  lifetime  of  the  level  into  which  the  isomeric 
population  is  to  be  dumped  by  the  trigger  is  too  short.  Large 
scale  pulsed-power  devices  typically  deliver  their  outputs  in 
pulses  with  duration  of  nanoseconds  to  a  few  microseconds,  so  it 
would  be  fortunate  to  find  a  laser  candidate  with  strong 
fluorescence  lifetimes  of  comparable  magnitudes.  In  this  case  as 
well,  the  ""Hf-2  is  favorable  having  transitions  with  lifetimes  of 
both  tens  of  nsec  and  tens  of  psec,  as  will  be  seen  in  the 
following  section. 

Since  the  experiments  needed  to  confirm  our  pragmatic  scoring 
of  the  better  of  the  29  candidates  for  a  gamma-ray  laser  have  been 
postponed  by  the  difficulties  in  obtaining  samples,  it  has  seemed 
useful  to  identify  simulations  which  can  serve  as  vehicles  with 
which  to  develop  systems  and  instrumentations  for  triggering.  To 
be  useful  a  simulation  nuclide  needs  to  have  a  strong  output 
fluorescence  transition  that  can  be  reached  by  cascade  from  a  K- 
breaking  level  pumped  from  the  ground  state  by  intense  x-ray 
pulses.  Moreover  the  levels  radiating  such  transitions  should  have 
lifetimes  of  tens  of  nsec  to  tens  of  psec.  Four  convenient 
simulation  nuclides  are  shown  in  Table  II. 

Recently  we  have  studied  the  pumping  of  the  4  simulation 
nuclides  shown  in  Table  II  (Hong  1993) .  Excitation  was  provided  by 
our  4  MeV  linac  which  produces  150  mA  pulses  of  electrons  of  3  ps 
duration  at  a  repetition  rate  of  360  Hz.  These  are  routinely 
converted  to  bremsstrahlung  with  a  cooled  tungsten  converter.  The 


Compton  scattering  of  this  bremsstrahlung  from  a  milligram-sized 
target  support  usually  produces  the  main  background  noise  tending 
to  obscure  gamma  fluorescence  from  the  nuclei  being  irradiated. 

Our  present  system  has  used  a  gated  detector  to  reduce  the 
sensitivity  during  the  pump  period.  A  combination  of  issues  limits 
the  recovery  of  the  system  to  times  after  a  certain  minimum  product 
of  photon  energy  and  delay  time  for  the  fluorescence  signal.  As 
currently  implemented  this  limit  is  150  keV  x  2  ps.  Shown  in  Fig. 
5  is  data  from  the  delayed  fluorescence  of  mTa  detected  at  482  keV 
with  a  recovery  time  from  the  end  of  the  x-ray  pulse  of  7.22  ps. 
The  literature  value  of  the  lifetime  of  the  fluorescent  state  being 
pimped  is  18  ps  which  agrees  well  with  the  decay  constant  seen. 
The  energetics  of  the  process  are  summarized  in  the  diagram  in  Fig. 
6.  Similar  results  are  shown  in  Fig.  7  for  the  nuclide  l7,Hf  for 
which  the  relevant  energy  levels  are  summarized  in  Fig.  8. 
Analyses  of  the  fluorescence  data  using  Eqs.  (1)  -  (4)  yielded  new 
values  for  the  integrated  cross  sections  for  pumping  short-lived 
fluorescence.  Figure  9  shows  these  results  included  in  the  summary 
plot  of  Fig.  4.  It  is  encouraging  to  observe  that  the  systematics 
for  the  pumping  of  these  simulated  laser  levels  closely  follows  the 
trends  already  established  for  the  excitation  of  longer  lived 
levels  that  were  easier  to  study.  As  yet  the  complementary  data  on 
the  excitation  energies  needed  to  reach  the  pumping  gateways  has 
not  been  measured  in  these  simulations,  but  they  might  reasonably 
be  expected  to  also  continue  the  established  systematics. 


3.2  First-Ranked  Candidate 

The  first  priority  candidate  for  a  gamma  ray  laser  is  the  31- 
year  isomer  of  17,Hf  shown  in  the  first  line  of  Table  I  and  it  is 
superior  to  the  next  possibilities  by  orders-of-magnitude.  The 
difficulty  has  been  that  macroscopic  samples  have  •  *  been 
available. 

The  energy  level  diagram  for  l,,Hf  is  shown  in  Fig.  10  with 
prominent  fluorescent  transitions  indicated.  The  fundamental 
question  is  whether  the  K-breaking  band  will  be  found  at  the  i  vel 
shown  in  the  figure  that  is  predicted  by  systematics. 

In  case  of  a  weak  success  in  which  K-breaking  is  tending  to 
fail  we  would  expect  the  preferential  population  of  the  band  to  the 
right  of  the  figure  leading  to  fluorescence  at  437  keV  with  35% 
efficiency  from  the  bandhead  having  68  ps  lifetime.  This  is  well 
within  current  capabilities  demonstrated  by  the  simulations. 

In  case  of  a  strong  success  in  which  K-breaking  is  complete  we 
would  expect  to  preferentially  populate  the  left-side  band  which 
leads  to  fluorescence  at  922  keV  with  65%  efficiency  and  1247  keV 
with  30%  efficiency  from  the  bandhead  having  78  nsec  lifetime.  In 
this  case  demonstrations  of  dumping  the  populations  would  need  to 
use  an  accelerator  with  shorter  pulses  such  as  APEX-I  in  our 
facility  which  produces  30  kA  of  electrons  in  pulses  of  30  nsec 
having  endpoint  energies  variable  from  0.7  to  1.2  MeV.  Although  a 
single  shot  device,  it  served  well  in  our  first  investigations  of 
gamma  fluorescence  (Anderson  &  Collins  1988b) .  There  is  a  great 
advantage  in  the  use  of  this  device  since  it  can  be  configured  to 
pump  the  left-side  band  with  bremsstrahlung  having  0.7  to  0.8  MeV 


endpoint  that  will  give  signal  fluorescence  at  1.25  MeV.  It  should 
be  straightforward  to  set  discrimination  levels  in  the  electronics 
to  reject  the  scattered  pump  radiation  since  it  consists  of  such 
lower  energy  photons. 

Using  the  model  of  the  preceding  section  we  can  calculate  the 
expected  number  of  fluorescent  photons  with  two  assumptions  about 
the  overall  integrated  cross  section  for  pumping  to  fluorescence, 
(oDsy,  =  5  X  10**s  cm2  keV  ,  and  (7a) 

( or =  2.5  X  10«  cm*  keV,  (7b) 

where  the  first  is  a  number  consistent  with  the  measurements  of 
Fig.  9  and  the  second  is  scaled  for  the  dependence  on  pump  energy 
given  by  Eq.  (6)  .  The  range  of  excitation  energies  spanned  by  the 
measurements  of  Fig.  9  is  insufficient  to  determine  any  correlation 
between  (or)  and  the  energy  needed.  As  can  be  seen  from  the 
results  summarized  in  Table  III,  there  is  an  ample  safety  margin 
for  the  arrangement  of  realistic  geometries  for  the  collection  and 
detection  of  the  fluorescent  photons  if  samples  of  the  sizes  shown 
can  be  acquired. 

4 .  Conclusions 

In  the  most  recent  years  a  breakthrough  in  research  on  the 
feasibility  of  a  gamma-ray  laser  has  proven  the  pervasive 
occurrence  of  K-breaking  levels  which  can  be  used  to  dump  the 
energies  stored  in  isomeric  nuclei  by  pumping  them  with  x  rays. 
Studies  of  systematics  have  shown  these  levels  to  lie  between  2.5 
-  3.0  MeV  for  nuclei  lying  in  the  range  of  masses  spanned  by  *‘7Er 
and  1,sPt.  An  extension  of  those  studies  being  reported  here  has  now 
shown  that  the  favorable  occurrence  of  such  giant  resonances  for 


pumping  nuclear  fluorescence  extends  to  the  excitation  of  laserlike 
transitions  having  lifetimes  between  ns  and  ps.  These  latest 
results  have  provided  further  guidance  in  the  rankings  of  the  29 
candidate  isomers  for  a  gamma-ray  laser  that  have  served  to 
identify  the  one  most  likely  to  be  successfully  pumped,  the  31-year 
isomer,  "'Hf"2. 

While  the  paucity  of  actual  samples  has  inhibited  the  final 
demonstration  of  feasibility,  estimates  of  likely  yields  from  such 
experiments  have  shown  that  strong  levels  of  fluorescence  from 
laserlike  levels  of  the  31-year  isomer  l7,Hf“2  should  be  observable 
when  samples  become  available.  These  conclusions  would  seem  to 
provide  a  strong  motivation  for  the  production  of  this  first  ranked 
candidate  for  use  in  a  gamma-ray  laser. 


References 


Anderson,  J.  A.  &  Collins,  C.  B.  1987  Rev.  Sci.  Instrum.  58, 
2157. 

Anderson,  J.  A.  et  al.  1988a  Phys.  Rev.  C  38,  2833. 

Anderson,  J.  A.  &  Collins,  C.  B.  1988b  Rev.  Sci.  Instr  59,  414. 

Carroll,  J.  J.  et  al.  1991  Phys.  Rev.  C  43 

Collins,  C.  B.  et  al.  1982  J.  Appl.  Phys.  53,  4645. 

Collins,  C.  B.  et  al.  1988  Phys.  Rev.  C  37,  2267. 

Collins,  C.  B.  et  al.  1990  Phys.  Rev.  C  42,  R1813. 

Collins,  C.  B.  1991  Handbook  of  Laser  Science  and  Technology  (CRC 

Press,  Boca  Raton,  FL) ,  pp.  561-567. 

Collins,  C.  B.  et  al.  1992  Laser  Interaction  and  Related 

Plasma  Phenomena,  G.  H.  Miley  and  H.  Hora,  eds.  (Plenum,  New 
York)  Vol .  10,  pp.  151-166. 

Collins,  C.  B.  et  al.  1993  Laser  and  Particle  Beams  11,  43. 

Hong,  C.  Ph.D.  Dissertation  1993  University  of  Texas  at  Dallas. 

Nelson,  W.  R.  et  al.  1985  The  EGS4  Code  System,  Stanford  Linear 
Accelerator  Center  Report  No.  SLAC  265  (Stanford  Linear 
Accelerator  Center,  Stanford,  CA) . 

Oganessian,  Yu.  Ts.  et  al.  1992  J.  Phys.  G  18,  393. 

Walker,  P.  M.  et  al.  1990  Phys  Rev.  Lett.  65,  416. 


Table  I:  Examples  of  nuclear  properties  for  three  of  the  29 
candidate  isomers  for  a  gamma-ray  laser  that  are  of  pragmatic 
importance  in  estimating  the  probabilities  for  the  successful  use 
of  that  material. 


ISOMER 

ENERGY 
DENSITY 
( Joules/pg) 

SHELF 

HALFLIFE 

TRIGGER 

PHOTON 

(MeV) 

l7*Hf>2 

1300 

31  years 

<0.5 

179HfM2 

600 

25  days 

-1.5 

1,0Ta“ 

40 

101S  years 

-2.5 

Table  II:  Summary  of  properties  of  nuclides  useful  in  simulating 
the  pumping  of  the  better  of  the  29  candidates  for  a  gamma-ray 
laser. 


Nuclide 

Abundance 

(%) 

Output 

transitions 

(keV) 

Fluorescence 

lifetime 

(ps) 

•lBr 

49.31 

260 

37.6 

276 

17<Hf 

5.206 

736 

9.9,  9.6 

1043 

mTa 

99.988 

346 

18 

482 

201Hg 

13.1 

521 

92 

Table  III:  Summary  of  expected  numbers  of  excitations  from 
microweight  samples  of  the  31-year  isomer  of  l7,Hf  pumped  with  the 
UTD  linac  in  a  100  %  coupling  geometry. 


Isomeric  (17,Hf) 
Nuclei 

(or) 

Linac  (3.3  x  10“ 
photons/cm2  per 
hour) 

3  x  10“ 

Systematic 

2.3  x  104 

Best 

1.1  x  10‘ 

10“ 

Systematic 

7.6  x  10‘ 

Best 

3.8  x  10‘ 

3  x  10“ 

Systematic 

2.3  x  10* 

Best 

1.1  x  107 

10“ 

Systematic 

7.6  x  10s 

Best 

3.8  x  107 

Captions 

Figure  1 :  Schematic  representation  of  the  decay  modes  of  a  gateway 
state  of  width  T  sufficiently  large  to  promote  bandwidth  funneling. 
The  initial  state  from  which  population  is  excited  with  an 
absorption  cross  section  o0  can  be  either  ground  or  isomeric. 

Figure  2:  Activation  efficiencies.  A,  for  the  reaction 
,7Sr  (y, y' )  ,1Sr-  are  shown  as  functions  of  the  endpoint,  E0  of  the 
brerasstrahlung  used  for  excitation.  The  solid  curve  plots  values 
computed  from  Eq.  (3)  using  gateway  parameters  found  in  the 
literature  together  with  calculated  photon  spectra.  The  figure 
compares  these  Af  values  with  measurements  obtained  from  four 
different  accelerators. 

Figure  3:  Activation  efficiencies.  A*  for  the  reaction 
»°Ta“(Y,  Y'  )180Ta  are  shown  as  functions  of  the  endpoint,  E0  of  the 
bremsstrahlung  used  for  excitation.  Gateway  parameters  were  found 
from  calculated  photon  spectra  using  Eq.  (3) . 

Figure  4:  Integrated  cross  sections  for  pumping  isomeric  nuclei 
obtained  in  the  survey  of  Carroll  et  al.  1991.  The  groupings  of 
pumping  strengths  seen  in  the  figure  correspond  to  mass  islands 
between  magic  numbers  for  neutrons  and  protons.  The  best 
candidates  for  a  gamma-ray  laser  lie  within  the  island  containing 
the  largest  values  of  integrated  cross  sections  corresponding  to 
giant  pumping  resonances.  Within  each  island  the  integrated  cross 


sections,  and  the  gateway  excitation  energies  plotted  with  the 
right-hand  axis,  vary  only  slowly  with  changing  mass  number,  A. 

Figure  5:  Fluorescence  detected  from  decay  of  the  615  keV  isomer 
of  mTa.  The  data  were  obtained  using  a  gated  detector  system  and 
correspond  to  nearly  10*  acquisition  cycles.  Each  cycle  began  when 
the  detector  was  gated  on  7.22  ps  (t  *  0  on  the  horizontal  axis) 
after  the  end  of  a  3  ps  pulse  of  x  rays  from  the  linac,  as 
illustrated  in  schematic  form.  The  solid  line  shows  a  fit  to  the 
data  which  agrees  well  with  the  literature  value  of  18  ps. 

Figure  6:  Schematic  diagram  showing  the  energy  levels  important 
for  the  excitation  of  the  615  keV,  18  ps  laserlike  level  in  mTa. 
The  large  upward  arrow  indicates  the  absorption  transition  to  a 
gateway  located  near  2.8  MeV  which  populates  the  laserlike  level  by 
a  branch  of  its  decay  cascade  (shown  by  the  downward  dashed  arrow) . 
Transitions  providing  photons  used  to  detect  the  decay  of  the  18  ps 
state  are  indicated  with  their  energies. 

Figure  7:  Fluorescence  detected  from  decay  of  the  1333  keV  isomer 
of  l7*Hf .  The  data  were  obtained  using  a  gated  detector  system  and 
correspond  to  nearly  10s  acquisition  cycles.  Each  cycle  began  when 
the  detector  was  gated  on  7.22  ps  (t  =  0  on  the  horizontal  axis) 
after  the  end  of  a  3  ps  pulse  of  x  rays  from  the  linac.  The  solid 
line  shows  a  fit  to  the  data  which  agrees  well  with  the  literature 
value  of  9.6  ps. 


Figure  8:  Schematic  diagram  showing  the  energy  levels  important 
for  the  excitation  of  the  1333  keV,  9.6  ps  laserlike  level  in  n‘Hf. 
The  large  upward  arrow  indicates  the  absorption  transition  to  a 
gateway  located  near  2.8  MeV  which  populates  the  laserlike  levels 
by  branches  of  its  decay  cascade  (shown  by  the  downward  dashed 
arrows) .  Transitions  providing  photons  used  to  detect  the  decay  of 
the  9.6  ps  state  are  indicated  with  their  energies. 

Figure  9:  Integrated  cross  sections  for  pumping  laserlike  levels 
along  with  those  of  Fig.  4  for  the  pumping  of  long-lived  isomers. 
Again  the  groupings  correspond  to  mass  islands  between  magic 
numbers  and  the  best  candidates  are  identified  as  neighbors  of  the 
nuclide  180Ta. 

Figure  10:  Energy  level  diagram  of  l1*Hf.  The  spontaneous  decay  of 
the  31-year  isomeric  state  does  not  excite  transitions  in  the  left¬ 
side  band  nor  is  it  able  to  excite  the  K  =  14  bandhead  which 
cascades  through  the  right-side  band  during  its  68  ps  half life. 
The  K-mixing  level  inferred  from  systematics  is  shown  by  the  wide 
shaded  band  near  2.8  MeV  and  the  large  upward  arrow  symbolizes  the 
pumping  step  from  the  31-year  isomer.  Transitions  providing 
photons  useful  for  detection  of  the  dumping  reaction  are  identified 
along  with  their  energies. 
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Continuing  the  systematic  investigation  into  the  photoexcitation  of  isomers  over  wide  mass  and  energy 
ranges,  the  production  of  77Se'",  ’’Br'",  and  1  J7Ba*"  was  studied  with  the  bremssirahlung  facility  at  the  su¬ 
perconducting  Darmstadt  linear  accelerator.  These  isomers  have  half-lives  on  the  order  of  seconds  Ex¬ 
citation  functions  were  measured  for  the  (y.y'l  reactions  populating  the  metastable  states  for  energies  of 
2-7  MeV  and  the  important  intermediate  states  were  identified.  Nuclear  structure  calculations  with  the 
quasiparticle-phonon  model  for  7,Br  and  the  particle-  (hole-)  core  coupling  approach  for  "  Ba  gave  satis¬ 
factory  descriptions  for  the  strength  and  position  of  the  dominant  mediating  levels  Admixtures  of  frag¬ 
mented  outershell  single-particle  strength  shifted  to  low  energies  were  identified  as  essential  features  of 
the  wave  functions  of  those  states.  Intermediate  states  in  77Se  displayed  very  large  strengths  compared  to 
other  isomers  in  the  same  mass  region,  providing  further  support  for  the  correlation  between  integrated 
cross  sections  and  ground  state  deformations  recently  discovered  in  the  A  =  160-200  mass  region.  Such 
an  enhancement  would  considerably  improve  the  feasibility  of  a  gamma-ray  laser  based  on  the  sudden 
deexcitation  of  isomeric  populations  in  deformed  nuclei. 

PACS  number(s):  25.20.Lj,  27.50.  +  e.  27.60. +j,  21.90.  +  f 


I.  INTRODUCTION 

It  has  been  known  for  more  than  fifty  years  that  nu¬ 
clear  isomers  could  be  populated  by  (y,y')  reactions,  yet 
the  greatest  progress  in  the  study  of  this  process  and  the 
underlying  nuclear  structure  has  occurred  in  the  past  de¬ 
cade.  To  a  large  degree  this  renaissance  has  been  driven 
by  research  directed  towards  determining  the  feasibility 
of  a  gamma-ray  laser  whose  pump  scheme  [1,2]  would 
rely  on  the  resonant  deexcitation  of  long-lived  nuclear 
isomers  through  intermediate  states  (IS).  Such  research 
has  required  the  systematic  development  of  a  database 
describing  the  population  and  depopulation  of  metastable 
levels.  This  effort  began  in  earnest  [3-10]  in  1987  and 
has  recently  concentrated  [  1 1  - 16]  on  the  investigation  of 
IS  in  the  largely  unexplored  excitation  energy  range  of 
2-6  MeV. 

Applications  notwithstanding,  the  photoactivation  of 
isomers  is  important  in  its  own  right.  For  example,  the 
integrated  cross  sections  measured  for  IS  can  reveal  in¬ 
teresting  nuclear  structures  [17-19]  since  both  strong 
absorption  transitions  from  the  ground  state  and 
significant  cascades  to  the  isomer  are  needed.  In  this  re¬ 
gard  the  combination  of  photoactivation  and  nuclear  res¬ 
onance  fluorescence  (NRF)  experiments  has  proven  to  be 
a  particularly  useful  tool  and  the  resulting  information 
has  been  successfully  compared  with  microscopic  calcula¬ 
tions  [13,16].  For  deformed  nuclei  a  remarkable 
discovery  [5,11,14]  was  the  appearance  of  low-lying  levels 
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which  efficiently  couple  ground  states  and  isomers 
differing  by  as  many  as  eight  units  in  the  quantum  num¬ 
ber  K,  the  projection  of  the  total  angular  momentum  on 
the  major  nuclear  axis.  The  properties  of  IS  are  impor¬ 
tant  in  nuclear  astrophysics  when  stellar  nucleosynthesis 
of  a  radioisotope  is  expected  to  branch  to  both  ground 
state  and  isomer  which  have  greatly  differing  lifetimes 
[9,11,20].  In  a  practical  field,  nuclides  with  well-known 
IS  can  be  used  to  characterize  intense  bremsstrahlung 
produced  by  linear  accelerators  or  pulsed-power 
machines  [21-23]. 

The  systematic  investigation  of  the  photoexcitation  of 
isomers  was  continued  in  this  work  by  the  study  of  77Se, 
7,Br,  and  ,37Ba  whose  isomers  have  half-lives  on  the  or¬ 
der  of  seconds.  In  addition  to  the  experimental  data,  nu¬ 
clear  structure  calculations  are  presented  for  79Br  and 
137Ba  which  permitted  a  detailed  interpretation  of  the  in¬ 
terplay  between  single-particle  and  collective  properties 
that  provided  the  observed  intermediate  states. 

II.  EXPERIMENTS 
A.  Methods 

The  experiments  were  performed  using  the  injector  of 
the  superconducting  continuous  wave  electron  accelera¬ 
tor  (S-DALINAC)  at  Darmstadt  [24].  Bremsstrahlung 
for  the  experiments  was  produced  by  electrons  which  im¬ 
pinged  on  a  3  mm  thick  tantalum  converter  foil.  The 
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TABLE  I  Summary  of  the  !ilera(ure  values  [25]  for  the  relevant  nuclear  parameters  and  transparencies  .'or  the  escape  of  fluores¬ 
cence  jphotons  from  samples  of  the  nuclides. 


Nuclide 

g  s.  spin 

o 

Isomer  spin 

/, 

Abundance 

Principal 

T  ransparencs 

j ■ 

tkeV) 

J" 

1S> 

<%) 

fluorescence 

l %  i 

IkeVi 

”Se 

i 
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7  • 

T 

17  4 

7.6 
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76 

”'Br 

2 
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9  • 

2 

4  86 

50.69 

207  2 

76 

"7Ba 

i  • 

662 

II 

1 53  12 

1 1  23 

661.66 

85 

electron  energies  were  measured  before  and  after  caoh  ex¬ 
posure  to  an  accuracy  of  50  keV  and  were  varied  from 
E0  =  2  to  7  MeV  in  step  sizes  of  125  and  250  keV.  The 
proper  beam  alignment  was  insured  by  maximizing  the 
dose  registered  in  a  remote  ionization  chamber  which 
viewed  only  the  central  12  mrad  of  the  bremsstrahlung. 
Variations  in  all  beam  parameters  were  recorded  during 
each  exposure,  including  the  total  charge  passed  to  the 
converter. 

Targets  consisted  of  elemental  Se  and  compounds  of 
LiBr  and  BaF,  with  typical  masses  of  3-6  g.  The  materi¬ 
als  were  contained  in  hollow  aluminum  cylinders  with  3.5 
cm  length  and  1 .4  cm  outer  diameter.  The  samples  were 
irradiated  axially  with  their  plane  faces  located  1.3  cm 
from  the  converter  foil.  Nominal  electron  beam  currents 
were  5  fiA  and  the  exposure  periods  were  typically  twice 
the  half-life  of  the  isomer  being  investigated.  At  the  end 
of  each  irradiation  the  samples  were  transported  pneu¬ 
matically  to  a  well-type  Nal(Tl)  detector  for  counting  as 
described  in  more  detail  in  Refs.  [15,16].  To  improve  the 
statistical  accuracy  of  the  data,  measurements  were  re¬ 
peated  for  up  to  sixteen  cycles  for  each  isotope  and  elec¬ 
tron  energy.  A  sample  of  elemental  In  was  also  irradiat¬ 
ed  at  some  energies  for  calibration. 

The  numbers  of  isomers  produced  by  these  exposures 
were  determined  from  the  counting  rates  measured  in  dis¬ 
tinctive  fluorescence  lines.  The  y-ray  signatures  and  oth¬ 
er  relevant  parameters  for  the  isomers  investigated  here 
are  summarized  in  Table  I.  The  raw  numbers  of  counts 
within  peaks  observed  in  pulse-height  spectra  were 
corrected  in  a  standard  fashion  [15]  to  account  for  the 
isomer  half-lives,  the  detection  efficiencies,  and  the  opaci¬ 
ties  of  samples  to  the  escape  of  fluorescence  photons.  A 
typical  example  of  the  data  is  shown  in  Fig.  1(a)  and  was 
obtained  from  a  Se  sample  following  its  irradiation  with 
bremsstrahlung  having  an  end-point  energy  of  £0=4 
MeV.  The  isomeric  decay  signature  was  clearly  identified 
despite  the  relatively  poor  resolution  of  the  Nal  spec¬ 
trometer.  In  addition  to  pulse-height  spectra,  a  mul¬ 
tichannel  analyzer  was  used  to  obtain  decay  curves  like 
that  of  Fig.  1(b).  As  shown,  there  was  excellent  agree¬ 
ment  between  measured  and  adopted  [25]  values  for  the 
isomeric  half-lives. 

B.  Data  analysis 

Recent  investigations  have  confirmed  [11-16,26]  that 
the  population  of  isomers  by  (y,y')  reactions  at  excita¬ 
tion  energies  below  the  photoneutron  threshold  proceeds 
by  the  process  depicted  schematically  in  Fig.  2.  The 


figure  identifies  the  relevant  parameters  ineluding  the 
natural  width  of  the  intermediate  state  T  and  the  branch¬ 
ing  ratios  bu  and  blio  for  decay  from  the  IS  directly  to  the 
ground  state  and  by  unknown  cascade  to  the  metastable 
level,  respectively.  Along  with  the  Breit-Wigner  cross 
section  for  the  absorption  transition,  these  form  the  in¬ 
tegrated  cross  section  for  population  of  the  isomer 
through  the  IS, 


fa  2  Jls  +  I 

■’  1  vT*'t.r . 


where  Jls  and  J0  are  the  angular  momenta  of  the  IS  and 
ground  state,  respectively,  and  E  is  the  energy  of  the  ab¬ 
sorption  transition.  The  integrated  cross  section  is  ex¬ 
pressed  in  usual  units  of  10~29  cm2keV,  equivalent  to 
0.01  eVb.  This  quantity  is  of  primary  importance  for 
photoactivation  studies  since  the  normalized  activation 
Af  giving  the  fractional  yield  of  isomers  produced  by  an 


Channel 


Time  (s) 


FIG.  1.  (a)  Pulse-height  and  (b)  time-decay  spectra  obtained 
from  a  Se  sample  following  its  irradiation  with  bremsstrahlung 
having  an  end  point  of  4  MeV.  The  straight  line  in  (b) 
represents  a  best  fit  (log  scale)  with  Tt/2  =  17.3±0.2  s  which 
agrees  well  with  the  literature  value  [2S]  of  17.45  s. 
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FIG  2  Schematic  representation  of  the  resonant  photoexci- 
laiion  mechanism  by  which  an  isomer  at  an  energy  of  £1V1  is 
populated  through  ihe  hatched  intermediate  slate  t IS>  The  IS 
is  characterized  by  it  excitation  energy  E ,  and  its  natural  width 
T  The  branching  ratios  b0  and  blso  describe  the  decay  of  the  IS 
to  the  ground  state  and  to  the  isomer,  respectively.  The  latter 
quantity  is  the  sum  of  ail  unknown  cascades  through  the  dashed 
les  els  which  lead  to  the  isomer. 


irradiation  per  incident  photon  is 

A>s^a=Vtrr)laoFiE^Eo) '  l2) 

where  ,Vlso  is  the  number  of  isomers  produced,  A'r  is  the 
number  of  target  nuclei,  d>0  is  the  total  photon  flux  for 
the  specified  end-point  energy  E0,  and  F(EJtE0 )  is  the 
distribution  of  energies  within  the  bremsstrahlung  contin¬ 
uum  normalized  to  1.  The  summation  in  Eq.  (2)  allows 
the  participation  of  a  plurality  of  well-spaced  gateways. 

Values  for  <l>0  and  F(Ej,E0)  were  obtained  for  each 
end-point  with  the  well-established  EGS4  electron-photon 
transport  code  [27],  as  described  in  detail  elsewhere 
[11-16].  The  calculated  spectra  were  validated  with  the 
ll5In(y,y')ll5Inm  reaction  which  is  characterized 
sufficiently  to  serve  as  a  calibration  standard  [13,23], 

III.  RESULTS 


End-point  Energy  (MeV) 


FIG.  3.  Normalized  yields  A,  as  functions  of  electron  energy 
for  the  populations  of  (a)  17Se’",  (b)  7‘,Br"’,  and  (c)  IJ,Bam.  The 
asterisks  represent  experimental  results  of  Ref.  [10]  obtained  at 
6  MeV  in  a  completely  different  environment  and  agree  well 
with  the  present  data.  The  solid  lines  correspond  to  excitation 
functions  calculated  using  Tq.  (2)  with  the  intermediate  states 
listed  in  Table  II  and  the  dashed  lines  describe  the  error  bounds 
of  this  approach. 


Figure  3  displays  excitation  functions  of  Af  vs  end 
point  for  the  production  of  77Sem,  79Brm,  and  niBa  m.  For 
comparison,  the  asterisks  represent  the  experimental  re¬ 
sults  obtained  at  E0= 6  MeV  from  Ref.  [10].  The  good 
agreement  observed  between  values  measured  in  com¬ 
pletely  different  experimental  environments  established 
confidence  in  the  procedures.  Intermediate  states  were 
identified  by  sharp  increases  in  the  normalized  activation. 
Integrated  cross  sections  were  then  extracted  by  fitting 
Eq.  (2)  to  the  data  in  such  a  way  as  to  minimize  the  resi¬ 
due, 


RMlE0)=AflE0)-  j  (af) LFIEj,E0) 


Ej-E, 


(3) 


by  successively  removing  contributions  from  the  M 
lowest-lying  IS. 

End-point  energies  below  2  MeV  were  not  accessible  in 
the  present  work  so  previously  measured  experimental 
values  for  IS  at  lesser  energies  were  included  in  the  fitting 
procedure.  In  instances  where  literature  values  were  una¬ 
vailable  or  insufficient  to  reproduce  the  excitation  func¬ 


tions  of  Fig.  3  up  to  the  energy  of  the  first  visible  IS,  a 
single  intermediate  state  was  introduced  at  an  excitation 
energy  below  2  MeV.  Its  position  and  integrated  cross 
section  were  chosen  to  remove  the  baseline  level  of  ac¬ 
tivation  and  did  not  necessarily  correspond  to  that  of  a 
real  level.  Properties  of  IS  which  reproduced  the  excita¬ 
tion  functions  measured  here  are  summarized  in  Table  II. 

The  77Se  excitation  function  at  low  energies  (£0<  2.6 
MeV)  was  partially  reproduced  by  including  the 
30  X  10-29  cm2  keV  IS  at  1.005  MeV  of  Anderson  et  al. 
[4],  but  neglecting  the  small  contributions  of  the  lower- 
lying  levels  listed  there.  The  missing  strength  below  2 
MeV  was  ascribed  to  an  IS  at  1.7  MeV  of  150X  10'29  cm2 
keV  which  served  to  remove  the  remaining  baseline  ac¬ 
tivation.  Boivin  et  al.  [28]  also  examined  the  production 
of  this  isomer  up  to  an  excitation  energy  of  1.7  MeV. 
Their  work  indicated  three  important  IS  at  1.000,  1.190, 
and  1.600  MeV  with  integrated  cross  sections  of 
3.2X10-29,  18X10-29,  and  52X  10-29  cm2  keV,  respec¬ 
tively.  However,  the  yield  calculated  with  these  values 
severely  underestimated  the  data.  It  is  noteworthy  that 
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TABLE  1 1  V allies  of  integrated  cross  sections  i o  I  i;v  and  excitation  energies  ol  ilie  intermediate 
states  most  important  in  the  production  of  these  isomers  by  (y.y'i  reactions  Values  needed  to  fit  the 
data  were  determined  in  this  work  by  minimising  the  residues  of  Eq  (51  V'alues  available  in  the  litera¬ 
ture,  where  available,  were  used  to  reproduce  the  excitation  functions  for  energies  below  that  of  IS 
clearly  identified  in  this  work  Strengths  listed  above  4  MeV  were  introduced  to  account  for  the  excess 
yield  observed  toward  higher  energies  and  did  not  necessarily  correspond  to  single  lev  Is  or  groups  of 
levels  No  error  is  given  lor  the  2  MeV  slate  in  11  Ba  since  this  level  was  simply  used  to  estimate  the 
maximum  strength  corresponding  to  the  baseline  activation  resulting  from  null  measurements  at  ' 

MeV 


Isomer 


E,  (MeV) 


(ef  l;v,  (  1 0  cm2  keV) 


1005 
17*01 
2.6:0  1 
127  0  1 
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0.761 
18+0.  I 
3.2:0  I 
4.410.2 


5.9 
65  +  6 
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7000*  1500 


Refs.  (3,4,29] 


2.0 

3.2+0  I 
4.5±0. 2 


<  15 
220  +  20 
8500+800 


when  their  integrated  cross  section  at  1.005  MeV  was  re¬ 
placed  by  that  of  Ref.  [4],  and  with  their  values  for  IS  at 
1.190  and  1.600  MeV,  good  agreement  was  obtained  be¬ 
tween  the  calculated  and  measured  excitation  function 
below  2.6  MeV. 

In  the  case  of  79Br,  the  well-known  level  [3,4,29]  at 
0.761  MeV  was  included  as  listed  in  Table  II.  The  litera¬ 
ture  does  not  support  the  identification  of  any  other  IS  up 
to  an  energy  of  1.7  MeV.  However,  after  removing  the 
contribution  from  the  above  state  the  shape  of  the  excita¬ 
tion  function  of  Fig.  3(b)  did  not  indicate  a  flat  baseline 
but  instead  followed  the  isochromat  [30]  of  an  IS  posi¬ 
tioned  near  2  MeV.  Therefore  a  level  was  concluded  to 
lie  at  1.8±0. 1  MeV  as  shown  in  the  table. 

The  data  displayed  below  3  MeV  in  the  excitation  func¬ 
tion  for  u7Ba  represent  upper  bounds  on  the  yield  ob¬ 
tained  from  null  measurements  for  fluorescence  from  the 
isomer.  This  small  baseline  was  consistent  with  an  IS  at  2 
MeV  with  less  than  15  X 10~29  cm2  keV  for  its  integrated 
cross  section.  For  comparison  a  level  at  1  MeV  with  less 
than  3  X 10"“ 29  cm2keV  would  have  also  been  allowed. 
The  literature  [31]  describing  excited  states  lying  below 
2.6  MeV  does  not  support  the  identification  of  any  IS  so 
it  was  concluded  that  the  first  intermediate  state  of 
significant  strength  was  located  at  3.2±0. 1  MeV. 

The  results  of  the  application  of  Eqs.  (2)  and  (3)  that  in¬ 
cluded  the  low-lying  levels  discussed  above  are  shown  as 
the  solid  lines  in  Fig.  3.  The  dashed  lines  define  bands 
describing  the  errors  in  the  values  given  in  Table  II  and 
those  quoted  in  the  literature. 

IV.  MODEL  CALCULATIONS 
A.  "Br 

The  photoexcitation  of  79Brm,  an  odd-even  nucleus, 
was  studied  in  the  quasiparticle-phonon  model  (QPM). 


The  method  for  treating  such  nuclei  in  this  theoretical 
framework  is  described  in  detail  in  Ref.  [32]  and  has 
been  successfully  applied  to  the  calculation  of  NRF  and 
isomer  photoexcitation  ]16]  in  89Y.  Calculations  of  the 
properties  of  the  81Br(y,y‘)8lBrm  reaction  are  also  avail¬ 
able  [19]  which  should  show  very  similar  structure  to 
those  for  79Br  since  both  nuclei  have  the  same  g.s.  and 
isomer  spins.  The  two  extra  neutrons  in  81  Br  are  not  ex¬ 
pected  to  significantly  change  the  average  field. 

Calculations  were  performed  with  the  wave  function 

*VUM)  =  CJ  a)M  +  2  Df(Jv)(a]mQl,  )JM  )%  ,  (4) 

Aft  f 

for  states  with  angular  momentum  J  and  projection  M. 
In  Eq.  (4),  a]M  denotes  the  one  quasiparticle  and  Q]^,  the 
phonon  creation  operator  for  angular  momentum  A.,  pro¬ 
jection  fi,  and  random  phase  approximation  root  number 
i.  The  g.s.  wave  function  of  the  even-even  core  is  denot¬ 
ed  by  %  and  v  is  the  number  within  a  sequence  of  states 
of  given  Jv. 

The  effective  QPM  Hamiltonian  induces  the  average 
field,  pairing  interactions,  and  a  residual  interaction  be¬ 
tween  quasiparticles.  The  latter  was  adjusted  to  repro¬ 
duce  the  energy  and  strength  of  the  2*  and  3,+  phonons 
in  the  neighboring  even-even  nuclei,  A  Woods-Saxon  po¬ 
tential  with  the  parametrization  of  Ref.  [33]  was  used  to 
describe  the  average  field.  Phonons  with 
/,r=l±,2±,3±,4±,5-,6+,  including  the  collective  as  well 
as  two-quasiparticle  states,  were  considered  in  the  sum¬ 
mation  of  the  second  term  of  Eq.  (4).  The  complete 
configuration  space  up  to  Ex  =  12  MeV  was  taken  into  ac¬ 
count  to  avoid  truncation  effects. 

Numerical  calculations  were  performed  with  the  code 
PHOQUS  [34].  The  calculations  were  restricted  to  one- 
phonon  coupled  states,  since  the  total  electromagnetic 
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strength  is  known  to  be  unaffected  by  the  inclusion  of 
multiphonon  coupled  configurations.  An  effective  spin  g 
factor  g’.„  =  (0. 8)gfrcr  was  used.  Effective  charges 
ecti  ~  <  A‘  /  4  )e  and  e"g~  —(Z  /  A)e  were  introduced  for 
£  1  transitions  to  separate  the  central  mass  motion 

B.  mBa 

The  single-hole  configuration  of  the  unpaired  neutron 
in  |,7Ba  relative  to  the  A'  =82  shell  closure  suggested  that 
the  hole-core  coupling  model  [35],  called  the  unified  mod¬ 
el  IUMI,  would  be  a  useful  approach.  In  the  comparable 
case  of  "'in  (one  proton  hole  relative  to  Z  =50).  en¬ 
couraging  results  were  obtained  both  for  the  description 
of  NRF  and  for  the  isomer  photoexcitation  [13].  Also  a 
detailed  UM  study  of  l43Nd  (one  particle  relative  to 
N  =82)  has  recently  been  reported  which  successfully  ac¬ 
counted  for  all  major  features  of  the  low-energy  spectrum 
[36]. 

The  configuration  space  consisted  of  lh  (one-hole) 
states  (with  respect  to  the  semimagic  l38Ba  nucleus)  and 
lp-2h  (one-particle-two-hole)  states  (with  respect  to 
l36Ba).  The  g.s.  wave  function  of  l36Ba  relative  to  the 
neutron  closed  shell  at  Ar  =  82  was  calculated  with  the 
pairing  interaction  of  the  Hamiltonian.  A  realistic  ap¬ 
proximation  to  the  ,36Ba  g.s.  wave  function  was  essential, 
because  important  contributions  to  the  electromagnetic 
excitations  would  have  been  missed  using  a  simple  pure 
two-hole  (j,  n  )„♦  configuration  as  suggested  by  the 
independent-particle  model. 

Single-particle  energies  were  taken  from  systematics  of 
the  odd-mass  Ba  nuclei  and  the  single-particle  space  was 
coupled  to  collective  2+  and  3“  phonons  in  the  neighbor¬ 
ing  even  nuclei  [37],  Excitation  energies  and  transition 
strengths  for  those  phonons  were  adopted  from  experi¬ 
mental  data.  The  energy  shift  between  lh  states  and  lp- 
2h  states  across  the  shell  closure  was  estimated  from  the 
neutron  separation  energy  differences  [35]  between  the 
137Ba  and  139Ba  nuclei. 

An  effective  charge  ejff=(1.5)e  and  effective  g  factor 
rff-<0.7lgU  were  used  for  the  calculations  of  elec¬ 
tromagnetic  transitions.  A  collective  limit  gR—Z/  A  was 
assumed  for  the  gyromagnetic  ratio.  As  was  the  case  [13] 
for  llsIn,  strong  transitions  were  hardly  affected  by  using 
other  gR  values,  e.g.,  gR  =0. 

The  results  of  the  UM  calculations  were  validated 
against  the  low-energy  spectrum  [31]  up  to  Ex  —2  MeV. 
The  isomer  energy  was  found  to  be  150  keV  too 

high,  but  this  was  sensitive  to  the  h  n/2  single-particle  en¬ 
ergy  and  the  3f  coupling  strength  in  l3<Ba.  Transitions 
between  low-lying  states  (£2  transitions  only)  were  repro¬ 
duced  within  a  factor  of  2  for  the  scarce  information 
available  in  the  literature. 

V.  DISCUSSION 
A.  "Se 

The  nucleus  77 Se  shows  a  very  complex  structure  at 
low  excitation  energies  [38,39]  and  no  calculations  are 
available  for  electromagnetic  transitions  above  1  MeV. 


The  actuation  function  in  lie  Atai  is  dominated  by  a 
transition  at  2.6  MeV  The  integrated  cross  section  of 
Table  II  may  be  converted  into  a  decay  width  and  the 
corresponding  reduced  transition  probabilities  arc 
fl(£l)>  2.9X10  1  \V  u  .  «  (  Af  1  )  >  0  1 7  W  .u.,  and 

R<E 2)^33. 0  W.u  .  depending  on  the  assumed  nature  of 
the  excitation  step  These  values  are  lower  limits  only, 
since  the  largest  possible  product  of  branching  ratios, 
f>06l5O  =  0.25  was  used.  Typically,  the  branching  ratio 
factor  could  be  expected  to  be  about  2-5  times  smaller 
than  this  value. 

It  is  interesting  to  note  that  yields  obtained  for  Se"1 
are  much  larger  than  those  in  other  isomers  in  the  same 
mass  region  (79Br,|i|Br,1'  Sr, '  A  e  A  large  g.s  deforma¬ 
tion  (£=—0.28)  has  been  derived  [39]  for  Se  from 
Coriolis  coupling  calculations  and  therefore  this  provides 
another  hint  for  the  peculiar  role  of  deformation  in  the 
variations  of  integrated  cross  sections  for  the  population 
of  isomers  in  different  nuclei  [  1 5). 

B.  )9Br 

The  ground  state  and  isomer  have  spin  structures  of 
j-  and  |  +  ,  respectively,  so  the  isomer  can  be  reached  by 
either  an  £2  excitation  to  a  J  ”  =  j  IS  followed  by  an  £  1 
decay,  or  the  reversed  scheme  of  an  £1  excitation  to  a 
7r=2+  IS  with  subsequent  E 2  decay.  Intermediate 
states  found  in  the  QPM  calculations  are  compared  to  the 
experimental  data  in  Fig.  4  and  show  this  behavior  Only 
two  strong  IS  are  observed  and  the  contributions  of  all 
other  levels  are  at  least  two  orders  of  magnitude  smaller. 

A  difficulty  in  the  model  results  was  the  poor  reproduc- 


FIG.  4.  Positions  and  strengths  of  intermediate  states  in  "Br 
measured  experimentally  (upper)  and  calculated  using  the 
quasiparticle-phonon  model  as  described  in  the  text  (lower). 
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non  of  the  isomer  excitation  energy,  1.270  MeV, 
which  strongly  distorted  the  energy-dependent  integrated 
cross  sections.  Therefore,  the  experimental  value  of 
0.127  MeV  was  used  as  the  basis  for  electromagnetic 
transitions.  A  similar  problem  occurred  in  the  MBr  cal¬ 
culations  [19],  There  it  was  shown  that  the  inclusion  of 
two-phonon  coupled  configurations  shifted  the  lowest  -  ' 
State  to  roughly  the  correct  position  without  strongly 
affecting  the  strength  of  the  intermediate  states. 

The  calculated  low-energy  IS  has  ,/r=7~,  in  agree¬ 
ment  with  the  established  0.761  MeV  level  Assuming 
the  experimental  1.8  MeV  state  also  to  be  J  -  l  ,  a  very 
good  correspondence  was  obtained  for  the  total  energy- 
independent  isomer  population  strength  below  Ex—  2 
MeV.  Such  a  fragmentation  would  have  naturally 
emerged  from  the  model  calculations  if  more  complex 
configurations  had  been  considered.  Thus  it  was  conclud¬ 
ed  that  all  IS  below  2  MeV  had  Jr=j~  and  that  the  in¬ 
termediate  strength  had  a  common  excitation-decay 
structure. 

The  microscopic  structure  was  extracted  from  the  ■ 
state  wave  function 

'P,({“)  =  (0.34)2/7/2  +  (0.85)(2/?3/2®2i+)+  ••• 

The  excitation  occurred  through  the  2 9  phonon  com¬ 
ponent  coupled  to  the  main  g.s.  configuration  while  the 
decay  to  the  isomer  took  place  in  a  2/7/2  —  lg,,/2  £1 
single-particle  transition. 

The  strong  experimental  IS  at  3.2  MeV  was  well  repro¬ 
duced  by  the  QPM.  The  calculated  oF  was  a  factor  of  3 
larger,  but  again  inclusion  of  many-phonon  states  would 
have  been  expected  to  bring  the  values  into  closer  agree¬ 
ment.  The  QPM  predicted  Jv~j+  for  the  IS  and  a  re¬ 
versed  excitation-deexcitation  sequence  compared  to  the 
|  -  state.  Analysis  of  the  wave  function 

.  • 

)=(0.38)2dJ/2  -HO. 90 )( lg9/2®2^)  •  •  • 

revealed  that  a  2pin—*2din  £1  transition  was  responsi¬ 
ble  for  the  excitation  and  phonon  coupling  to  the  dom¬ 
inant  isomer  configuration  provided  an  £2  decay. 

The  common  feature  in  the  intermediate  levels  of  the 
QPM  results  is  the  presence  of  single-particle  com¬ 
ponents  in  the  wave  functions  belonging  to  strength  hav¬ 
ing  a  centroid  energy  lying  much  higher.  These 
configurations  from  the  next  major  shell  are  strongly 
damped  by  coupling  to  collective  excitations  [32]  and  this 
pushes  fragments  to  low  energies.  Exactly  these  com¬ 
ponents  were  found  to  be  crucial  for  the  isomer  popula¬ 
tion  in  79Br. 


C.  IJ7Ba 

The  experimental  data  showed  evidence  for  only  a  sin¬ 
gle  IS  up  to  4.3  MeV.  Within  the  experimental  limits,  no 
contribution  was  observed  below  2  MeV  in  agreement 
with  the  literature  [31].  Thus,  photoexcitation  of  137Bam 
has  proven  to  be  particularly  important  for  the  brems- 
strahlung  characterization  technique  of  Ref.  [22]  since 
the  isomer  yield  directly  provides  the  photon  flux  at  a  sin¬ 
gle  energy. 


The  results  of  the  UM  calculations  are  displayed  in  the 
middle  part  of  Fig.  5  and  indicate  a  strong  transition  .it 
2.75  MeV.  Within  the  limitations  of  the  model  u  ap¬ 
peared  justified  to  assign  this  to  the  experimentally  deier- 
mined  level  450  keV  higher  in  energy  at  3.2  MeV  The 
integrated  cross  section  of  the  state  was  ovcrpredicled  In 
a  factor  of  about  2-3,  depending  on  whether  the  experi 
menial  or  model  excitation  energy  was  used.  However, 
similar  to  the  79Br  case  this  was  probably  due  to  the  basis 
truncation  and  larger  configuration  spaces  would  be  ex¬ 
pected  to  bring  the  results  closer  to  experiment. 

An  analysis  of  the  wave  functions  of  the  IS  identified 
2d \/2  — *2/5/2  single-particle  £1  transitions  as  the  com¬ 
mon  excitation  mechanism.  Because  of  the  large  spin 
difference,  AJ  =  4  between  g.s.  and  isomer,  the  decay 
occurs  by  complicated  cascades  which  cannot  be  easily 
summarized. 

A  particular  problem  of  the  UM  calculations  already 
discussed  for  the  M5In(y,y')mIn'”  reaction  was  a  sys¬ 
tematic  deviation  for  the  average  £1  and  £-2  transition 
strengths  with  respect  to  experimental  systematics  in  this 
mass  region  [40],  Good  correspondence  of  the  distribu¬ 
tion  shapes  was  obtained  by  comparing  the  reduced  tran¬ 
sition  probabilities  of  all  possible  transitions  up  to  4  MeV 
in  the  model  with  Ref.  [40].  However,  the  maximum  had 
to  be  normalized  with  factors  of  /(£1)  =  0.2  and 
/  ( £2 )  =  50  while  absolute  A/1  strengths  were  reproduced 
well. 

A  reanalysis  of  the  isomer  population  implementing 
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FIG.  5.  Positions  and  strengths  of  intermediate  states  in 
,,7Ba  measured  experimentally  (upper),  calculated  with  the 
unified  model  (middle)  and  calculated  with  the  unified  model  us¬ 
ing  empirically  introduced  correction  factors  for  the  average 
transition  strengths  (lower).  Those  factors  were  /(£1)=0.02 
and  /(£2)=50. 
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Che  normalization  factors  is  shown  in  the  lower  part  of 
Pig.  5.  Surprisingly  little  change  was  observed  with  the 
main  difference  being  some  additional  strength  at  3-3.5 
MeV.  The  E 2  correction  factor  seemed  unrealistically 
large  considering  how-  well  the  low-energy  transitions 
were  described.  It  is  important  to  note  that  many  of  the 
weaker  E 2  transitions  resulted  from  destructive  interfer¬ 
ence  of  the  collective  and  single-particle  parts  of  the  tran¬ 
sition  operator.  The  renormalization  factor  for  £1  tran¬ 
sitions  seemed  more  realistic  due  to  the  sensitivity  to 
smaller  amplitudes  of  the  wave  functions  which  would 
have  been  smeared  out  in  an  enlarged  configuration 
space.  Accordingly,  the  best  description  might  lie  some¬ 
where  between  the  two  extremes  displayed  in  Fig.  5. 
However,  the  strong  IS  at  2.75  MeV  which  likely  corre¬ 
sponds  to  the  experimental  3.2  MeV  level  was  hardly 
affected  by  these  variations  in  the  UM  calculation. 

VI.  CONCLUSIONS 

Excitation  functions  were  measured  for  the  photoac¬ 
tivation  of  77Sem,  7<)Brm,  and  n7Bam  using  bremsstrahlung 
with  end-point  energies  in  the  range  2-7  MeV.  This  per¬ 
mitted  the  identification  of  intermediate  states  and  the 
extraction  of  their  integrated  cross  sections. 

A  particularly  strong  IS  at  2.6  MeV  was  observed  for 
the  population  of  77Sem.  Its  transition  strength 
significantly  exceeds  10~1 2 3 4 5 6 7 8 9 10  W.u.  for  an  £1  excitation  or 
0.1  W.u.  for  an  Ml  excitation,  both  of  which  are  quite 
large  values.  In  fact  the  magnitude  of  77Sem  yields  com¬ 
pared  to  other  isomers  in  this  mass  region  provide  further 
support  for  a  correlation  between  integrated  cross  sec¬ 
tions  and  the  g.s.  deformation.  Such  a  correlation  was  re¬ 
cently  discovered  in  the  A  —  160-  200  mass  region  [15] 


and  is  important  to  gamma-ray  laser  research  since  many 
candidate  nuclei  are  well  deformed 

The  appearance  of  strong  intermediate  si. lies  icquires  a 
special  combination  of  properties  including  a  large  partial 
g.s.  width  and  significant  mixing  in  the  wave  function 
Thus  the  number  of  important  levels  is  typically  v c i > 
small  at  energies  below  about  5  MeV  as  confirmed  by  the 
model  results.  The  calculations  further  emphasized  t lie- 
importance  of  a  delicate  interplay  of  single-particle  arid 
collective  amplitudes  in  the  IS  wave  functions  In  all  ex¬ 
amples  analyzed  so  far  (Refs.  (13.19)  and  this  work'  the 
decisive  factor  has  been  the  appearance  of  fragments  oi 
outer-shell  single-particle  strength  at  low  energies  b 
would  be  of  interest  to  compare  these  conclusions  with 
findings  in  other  valence-shell  regions  and  in  nuclei  with 
large  g.s.  deformations.  Experimental  and  theoretical 
work  along  these  lines  is  underway  . 
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The  Mossbauer  effect  was  used  to  investigate  films  of  nanophase  diamond  (NPD)  implanted  with 
isolopically  pure  57Fe  at  a  dose  of  5x  101*  atoms/cm2  and  an  energy  of  20  keV,  When  defects  and 
voids  created  by  the  implantation  were  repaired  with  an  overcoating  layer  of  NPD,  the  recoil-free 
fraction  at  room  temperature  for  these  samples  was  found  to  be  /dll= 0.94^:0.06  with  a 
corresponding  Debye  temperature  of  1140  K.  This  relatively  high  value  for  /  makes  NPD  films  a 
promising  host  matrix  for  microgram  quantities  of  Mossbauer  isotopes. 


The  principal  figure  of  merit  for  a  matrix  in  which  Moss¬ 
bauer  nuclei  are  embedded  is  the  fraction  /  of  events  which 
can  occur  without  any  recoil  of  a  nucleus  emitting  or  absorb¬ 
ing  gamma  rays.  The  probability  for  the  recoilless  emission 
or  absorption  of  gamma  rays  is  related  to  the  mean-square 
vibrational  amplitude  of  an  atom  in  the  direction  of  observa¬ 
tion  and  is  given  by  the  relation  f=cxp(-k2(x2))  where 
k=2ir/\.  Since  the  wave  number  k  is  determined  for  a 
given  transition,  an  increase  in  /  can  only  be  accomplished 
by  decreasing  the  mean-square  vibrational  amplitude  of  at¬ 
oms  in  the  host  lattice.  Such  a  change  may  be  accomplished 
by  either  reducing  the  temperature  of  the  lattice,  or  by  in¬ 
creasing  the  interatomic  binding  energy.  In  practice,  the  latter 
can  be  done  by  substituting  the  Mossbauer  nuclei  of  interest 
into  a  rigid  lattice. 

In  the  Debye  theory  of  solids,  rigid  lattices  have  high 
Debye  oscillator  frequencies,  wD ,  and  correspondingly  high 
Debye  temperatures,  0D  related  by  kB0D=hu>D .  The  recoil- 
free  fraction  is  related  to  0D  and  the  temperature,  T  of  the 
sample  by1 
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where  £*  is  the  recoil  energy  for  a  free  nucleus.  Clearly  / 
increases  asymptotically  as  a  function  of  temperature  with 
increasing  0D .  For  pure  samples  of  materials  containing  only 
the  nuclei  of  interest  for  Mossbauer  spectroscopy,  typical 
Debye  temperatures  are  on  the  order  of  a  few  hundred  de¬ 
grees  Kelvin. 

In  contrast,  natural  diamond  has  a  Debye  temperature  of 
00=2230  K,  the  highest  in  nature.2  It  is  therefore  one  of  the 
most  attractive  of  the  host  lattices  with  the  potential  for  sup¬ 
porting  Mossbauer  nuclei.  Because  of  the  convenience  aris¬ 
ing  from  the  use  of  a  single  device  to  accelerate  ions,  sepa¬ 
rate  isotopes,  and  decelerate  the  desired  fraction,  our 
interests  have  been  focused  upon  implantation  as  a  means  of 
introducing  Mossbauer  nuclei  into  hosts.  Numerous  studies 


*>l*reiem  address:  Advanced  Materials,  Vut>  Inc.  IMO  Industries,  Garland, 
IX  75046-9015. 


involving  the  implantation  of  Mossbauer  nuclei  into  natural 
diamond  have  been  reported  in  the  literature.1-7 

Recently,  substantial  progress  has  been  made  in  the  de¬ 
velopment  of  different  types  of  pure  carbon  films  having 
many  of  the  properties  of  diamond.  These  films  appear  to  be 
combinations  of  the  various  polytypes  of  diamond  and  often 
have  mechanical  and  thermal  properties  similar  to  those  of 
natural  diamond. 

Nanophase  diamond  (NPD)  films  are  unique  among  the 
growing  class  of  these  new  diamondlike  materials  and  offer  a 
number  of  important  solid-state  properties.8-12  As  used  here, 
these  NPD  films  were  condensed  from  ions  of  carbon  pro¬ 
duced  by  the  laser  ablation  of  graphite  at  intensities  sufficient 
to  impart  kinetic  energies  in  excess  of  0.5  keV.  They  were 
prepared  in  vacuum  on  substrates  at  room  temperatures  with¬ 
out  hydrogen  or  other  catalysts  and  without  columnar  pat¬ 
terns  of  growth.  Such  layers  of  NPD  nominally  contain 
75%-85%  diamondlike  sp3  bonds  concentrated  into  nodules 
of  20-30  nm  in  diameter  as  determined  by  STM  and  con¬ 
firmed  by  TEM.13  The  bulk  properties  of  these  films  closely 
approach  those  of  natural  diamond  in  the  areas  of  thermal 
management  and  hardness.14  These  films  can  also  be  depos¬ 
ited  on  a  wide  variety  of  different  substrates.  In  this  work, 
we  investigate  the  utility  of  these  films  as  host  matrices  for 
targets  containing  microgram  quantities  of  Mossbauer  iso¬ 
topes  by  performing  comparative  measurements  of  the 
recoil-free  fractions  for  57Fe  implanted  into  a  film  of  NPD 
and  for  iron  foils  enriched  in  57Fe. 

The  57Fe  implanted  sample  was  prepared  by  depositing  a 
1  pm  thick  film  of  nanophase  diamond  on  a  Ti  substrate  with 
a  diameter  of  25  mm  and  a  thickness  of  3  mm.  Titanium  was 
chosen  since  its  bonding  properties  to  NPD  films  have  been 
characterized  and  reported  in  the  literature.13  The  bonding  of 
NPD  to  Ti  is  particularly  strong  because  of  the  formation  of 
an  interfacia]  layer  of  TiC  which  is  20-30  nm  thick.  After 
deposition,  the  film  was  ion  implanted  with  5X1016 
atoms/cm2  of  isotopicaliy  pure  57Fe  at  an  energy  of  20  keV. 
The  iron  atoms  penetrated  to  a  depth  of  approximately  20  nm 
in  the  film  as  determined  by  Rutherford  backscattering  spec¬ 
troscopy  (RBS)  and  represented  a  concentration  of  about 
10%.  With  an  effective  implantation  area  of  3.3  cm2,  a  total 
mass  of  15.5  mg  of  57Fe  was  present  in  the  sample. 


Appl.  Phys.  Lett.  64  (10).  7  March  1994 


0003-6951/94/64(10)11 221/3/96.00 


C 1994  American  institute  of  Physics  1221 


In  order  to  account  for  the  effects  of  57Fe  impurity  at¬ 
oms  in  the  NPD  matrix,  it  is  customary  to  define  an  effective 
Mossbauer  Debye  temperature,  0M  which,  for  a  monatomic, 
regular  lattice  is  given  by7  17 


0U=0, 


m  y 


(2) 


where  8t)  is  the  Debye  temperature  of  the  host  lattice,  and 
y0,  y  are  the  force  constants  for  the  interaction  between  the 
host  atoms  and  impurity  atoms  of  mass  m„  and  m.  respec¬ 
tively.  As  a  first  approximation,  we  assume  that  the  force 
constants  are  nearly  equal.15 16  In  the  case  of  the  implanted 
NPD  sample.  /»„=  12  amu  for  the  C  atoms,  and  m  =57  amu 
for  the  vFe  atoms.  Setting  y0=y,  we  find  that  the  Debye 
temperature  for  the  NPD  lattice  is  a  factor  of  2.18  more  than 
the  average  effective  Debye  temperature  of  the  iron  sites,  as 
would  be  determined  from  the  Mossbauer  effect 

In  addition  to  the  reduced  effective  Debye  temperature 
for  the  57Fe,  the  relatively  low  number  of  atoms  implanted  in 
the  diamond  film,  combined  with  the  presence  of  the  Ti  sub¬ 
strate,  made  the  use  of  standard  spectroscopic  techniques  re¬ 
lying  on  the  transmission  of  gamma  rays  impractical.  To  ob¬ 
serve  the  resonant  effects  produced  by  the  57Fe  nuclei 
required  the  use  of  a  backscattering  technique.  In  this  work, 
both  conversion  electron  Mossbauer  spectroscopy  (CEMS) 
and  photon  scattering  techniques  were  used  to  collect  data 
for  the  /  measurements.  The  CEMS  is  advantageous  because 
of  the  large  internal  conversion  coefficient  for  57Fe.  How¬ 
ever,  it  is  only  sensitive  to  a  small  surface  layer  in  the  sample 
with  a  thickness  on  the  order  of  the  mean  escape  depth  for 
electrons  in  the  material.  For  conversion  electrons  in  metallic 
iron,  this  depth18  is  57  nm  and  is  probably  greater  in  NPD. 
Since  the  thickness  of  the  implanted  layer  is  only  20  nm  the 
mean  escape  depth  is  greater  in  either  case. 

Mossbauer  spectra  were  collected  using  a  He-CH4  gas 
flow  proportional  counter  with  an  approximate  flow  rate  of  5 
ml  per  minute  trf  4%  CH4  in  He.  A  typical  Mossbauer  spec¬ 
trum  of  the  57Fe  in  NPD  sample  is  shown  in  Fig.  1  directly  as 
measured.  The  spectrum  consists  of  two  components;  a 
quadmpole  split  doublet  and  a  singlet.  The  doublet  appears 
to  arise  from  iron  atoms  located  in  low  symmetry  or  intersti¬ 
tial  sites  between  the  sp3  nodules  whereas  the  singlet  corre¬ 
sponds  to  substitutional  or  high-symmetry  sites  in  the  sp3 
hard  nodules  themselves.4,7 

A  relative  measurement  of  the  recoilless  fraction  for  the 
S7Fe  in  NPD  was  made  using  the  photon  scattering  geometry. 
This  arrangement  was  used  because  it  is  sensitive  to  the  en¬ 
tire  volume  of  the  sample,  unlike  CEMS.  The  measurement 
was  performed  by  comparing  the  resonant  scattering  inten¬ 
sity  from  the  implanted  NPD  sample  with  that  from  metallic 
iron  targets  having  a  known  recoil-free  fraction  of  /= 0.69  at 
room  temperature.19  The  same  ^Coffd)  source  was  used  for 
all  of  the  targets20,21  and  made  the  measurements  indepen¬ 
dent  of  its  recoil-free  fraction.  The  recoil-free  fraction  of  the 
NPD  sample  was  determined  by  comparing  the  experimental 
areas  under  the  resonance  peaks  in  the  spectra.  This  compari¬ 
son  is  based  upon  the  assumption  that  such  an  optically  thin 
sample  is  free  from  spurious  additional  absorptions  and  has 
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MG  1.  Conversion  cleclron  Mossbauer  spectrum  of  S  *  Ml’  aloms/em"  of 
"Fc  implanted  in  a  I  /im  thick  him  of  nanophase  diamond  After  implan 
tation  the  sample  was  overcoated  with  47  nm  of  NPD  to  partially  repair 
implantation-induced  damage.  The  spectrum  is  composed  of  three  lines,  a 
symmetric  doublet  with  a  quadrupolc  splitting  of  0  H50  mm/s  and  an  isomer 
shift  of  +  0  062  mm/s  due  to  atoms  in  low  symmetry  interstitial  sites  and  a 
singlet  with  an  isomer  shift  of  -0..726  mm/s  which  can  he  attributed  lo 
nuclei  in  high  symmetry  substitutional  sites  in  the  sp'  hard  nodules  The 
source  was  7Co(Pd)  with  a  nominal  activity  of  2  mCi 


the  additional  advantage  of  canceling  errors  due  to  purely 
instrumental  effects. 

Data  were  collected  for  a  variety  of  calibration  foils  hav¬ 
ing  different  masses  of  57Fe.  The  total  counts  in  the  reso¬ 
nance  lines  were  plotted  as  a  function  of  the  total  mass  of 
57Fe  in  each  sample.  This  calibration  curve  is  shown  in  Fig. 
2.  As  expected,  the  count  rate  increases  linearly  with  the 
amount  of  57Fe  in  each  sample.  The  total  signal  count  rate  for 
the  S7Fe  implanted  in  the  nanophase  diamond  film  is  plotted 
against  the  calibration  curve  and  direct  comparison  indicates 
that  the  fluorescent  signal  from  the  57Fe  in  the  nanophase 
diamond  is  greater  than  would  be  expected  for  a  correspond¬ 
ing  metallic  absorber  having  the  same  number  of  Mossbauer 
nuclei.  From  the  calibration  data  with  /Fe=0.69,  a  value  of 
/dia=0.94±0.06  is  obtained.  This  value  for  the  recoil-free 
fraction  corresponds  to  an  effective  Debye  temperature  of 


He  of  5TFe 


FIG.  2.  Calibration  curve  for  the  resonant  counting  rate  of  the  metallic 
samples  as  a  function  of  the  mass  of  Mossbauer  nuclei  coouined  in  them. 
From  these  data  the  lower  limit  on  the  recoil-free  fraction  of  the  nanophase 
diamond  film  is  /*,“0.94±0.06. 
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0m=523  K  for  (he  '7Fe  atoms  in  NPD.  Using  Eq.  (2),  we 
obtain  a  Debye  temperature  of  1140  K  for  the  NPD  film 
itself.  By  comparison,  the  Debye  temperature  for  metallic  Fe 
is  470  K  - 

Although  Ihe  measured  Debye  temperature  for  the  NPD 
lattice  is  much  greater  than  that  for  typical  metallic  host  ma¬ 
trices,  it  is  still  below  the  Debye  temperature  for  natural 
diamond  A  possible  explanation  for  this  is  that  the  implan¬ 
tation  of  s,Fe  nuclei  into  the  lattice  could  result  in  partial 
graphitization  of  the  v/»'  component  of  the  film.  The  sp' 
graphitic  bond  is  more  thermodynamically  stable  than  the 
sp '  diamondlike  bond  so  that  any  excessive  energy  trans¬ 
ferred  to  the  lattice  from  the  implantation  could  tend  to 
change  the  sp*  bonds  to  sp2  22 

For  natural  diamond,  implantation-induced  damage  "can 
be  partially  repaired  by  heat  treating  the  sample  either  during 
or  after  the  implantation.’’  We  have  found  that  the  deposition 
of  an  additional  thin  layer  of  nanophase  diamond  after  im¬ 
plantation  accomplishes  a  similar  effect,  probably  by  filling 
voids  by  the  condensation  of  additional  material.  In  this  case 
the  ease  with  which  implantation-induced  damage  can  be 
repaired  provides  an  additional  advantage  for  the  use  of  NPD 
films  as  host  matrices  for  small  amounts  of  Mossbauer  nu¬ 
clei. 

The  results  obtained  with  the  ,7Fc  implanted  NPD 
sample  closely  parallel  the  results  obtained  by  Sawicki  et  al.2 
In  that  work,  S7Co  was  implanted  into  natural  diamond  at  a 
dose  of  2x  10n  atoms/cm2  and  an  energy  of  50  keV.  The 
Debye  temperature  for  the  high  symmetry  substitutional  sites 
was  found  to  be  1300  K  while  that  for  the  low  symmetry 
interstitial  sites  was  550  K.  Our  result  of  0D- 1140  K  for  the 
NPD  film  shows  that  its  rigidity  is  comparable  to  that  of 
natural  diamond  in  spite  of  the  orders  of  magnitude  higher 
dose  of  5xl0IA  atoms/cm2  used  in  the  preparation  of  the 
sample. 

Comparing  the  best  of  the  data7  from  the  literature  for 
natural  diamond  with  the  result  for  NPD  shows  the  two  to  be 
reasonably  equivalent  in  terms  of  the  effective  Debye  tem¬ 
peratures  experienced  by  the  implanted  s7Fe  nuclei.  If  the 
NPD  offers  an  advantage,  it  is  to  be  found  in  cost  and  prac¬ 
ticality.  In  the  work  reported  here  a  3.3  cm2  area  of  NPD  was 


doped  to  a  concentration  of  10%  57Fe.  Such  high  concentra¬ 
tions  over  large  areas  of  natural  diamond  would  seem  to  lie 
outside  the  realm  of  practicality.  Moreover,  the  use  of  NPD 
may  support  the  implantation  of  other  Mossbauer  nuclei  such 
as  llvSn  and  l25Te.  If  comparable  successes  were  realized  the 
recoil-free  fractions  would  be  -0  7  at  3(H)  K  in  both  cases 
The  authors  gratefully  acknowledge  the  support  of  this 
work  by  SDiO/TNl  through  NRL  and  preparation  of  the 
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Thu  applicability  ol  the  Monte  t  ,nlo  piogram  GLANT3  for  the  calculation  ol  low  energy,  thick  taiget  htcmss'ir.ihlung  is 
investigated  The  results  ol  (it.ANT.V  version  15.  shoes  very  good  correspondence  to  those  of  the  well  established  f:GS-l  aide  In 
contrast  it  is  shown  that  GLANT3.  version  II.  leads  to  erroneous  angular  distributions,  total  intensities  anil  even  spectial  shapes 
for  angles  far  from  the  critical  angle  The  use  of  the  M'ln(y,-y')u4lnm  reaction  for  an  absolute  calibration  of  bremsstrahlting 
spectra  is  also  demonstrated  This  ntelhod  relies  on  the  reconstruction  of  the  isomer  yield  excitation  I unc lion  liom  a  sei  ol  In 
levels  with  well  characterized  properties  It  can  he  utilized  for  continuous  photon  spectra  up  to  about  3  Me V 


I.  Introduction 

Recently  considerable  interest  has  been  shown  in 
the  field  of  resonant  photoabsorption  experiments  us¬ 
ing  thick  target  brcmsstrahlung  with  energies  of  several 
MeV.  This  has  been  motivated  by  the  discovery  1 1)  of 
low-lying  collective  magnetic  dipole  transitions,  the  so 
called  “scissors  mode’’  [2],  as  well  as  the  observation  of 
enhanced  electric  dipole  transitions  in  heavy  deformed 
nuclei  [31.  For  both  problems  nuclear  resonance  fluo¬ 
rescence  (NRF)  provides  an  ideal  experimental  tool  [4J. 
Another  field  of  experimental  activity  is  the  investiga¬ 
tion  of  resonant  photoexcitation  of  isomers  [5]  which 
has  an  impact  on  fields  as  varied  as  -y-ray  laser  re¬ 
search  [6],  nuclear  astrophysics  [7-10]  and  nuclear 
structure  [11-14]. 

While  the  study  of  ( y,y ')  reactions  has  a  long  his¬ 
tory,  its  usefulness  has  always  been  hampered  by  the 
lack  of  experimental  methods  with  adequate  accuracy 
for  an  absolute  calibration  of  fluxes  and  spectral  distri¬ 
butions  of  intense  photon  continua.  Commonly  these 
continua  are  thick  target  bremsstrahlung  or  Compton 
spectra  of  strong  radioactive  sources.  It  has  therefore 
become  increasingly  popular  to  use  Monte  Carlo  (MC) 
simulations  which  provide  an  appropriate  solution  even 
for  complex  geometries  [15]. 
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In  NRF  experiments  independent  data  of  "refer 
ence”  nuclei  with  well  measured  transition  strengths 
covering  the  relevant  excitation  energy  region  can  be 
used  for  a  sampling  of  the  unknown  photon  spectra 
and  for  a  normalization  of  the  flux.  It  was  found  [16.17] 
that  a  simulation  with  the  MC  program  GEANT3  [  18] 
described  the  shapes  of  the  forward  thick  target 
bremsstrahlung  spectra  rather  well.  However,  investi¬ 
gations  of  the  photoexcitation  of  isomers  usually  in¬ 
volve  more  complex  geometries  which  lead  to  a  sensi¬ 
tivity  not  only  on  the  spectral  shapes,  but  also  on  the 
spectral  distributions.  Additionally,  the  total  flux  must 
be  known  to  calibrate  the  experimental  results  (see  e  g 
refs.  [5,13]).  In  this  work,  the  analysis  was  based  on 
another  electron-photon  shower  MC  program,  EGS4 
f’9].  The  code  is  widely  used  in  medical  physics  and  its 
basic  correctness  has  been  verified  experimentally  [20]. 

A  thorough  comparison  of  the  results  obtained  with 
both  programs  for  identical  geometries  revealed  mas¬ 
sive  differences  which  questioned  the  conclusions 
drawn  in  refs.  [16,17].  Therefore,  one  purpose  of  the 
present  work  is  a  detailed  presentation  of  results  from 
the  calculations  of  thick  target  bremsstrahlung  with  the 
MC  codes  GEANT3  and  EGS4.  In  particular,  results 
from  the  recently  released  version  15  of  GEANT3  are 
compared  to  those  of  the  earlier  version  11,  on  which 
the  calculations  described  in  refs.  [16,17]  were  based. 

In  the  second  part  an  experimental  method  for  the 
absolute  calibration  of  thick  target  bremssstrahlung 
spectra  applicable  up  to  endpoint  energies  E0  »  3  MeV 
is  presented.  It  relies  on  the  measurement  of  the  yield 
excitation  function  of  the  "sIn  isomer  obtained  in 
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hrcmsstruhlung  irrudi.ilions  with  x;u\ing  endpoint  en¬ 
ergy  ]  1 2.2 1 1  The  importance  ol  .leeur.ile  methods  for 
the  absolute  eahhralion  ol  continuous  photon  spectra 
might  he  interred  from  Fable  1  ol  ref  (22)  There,  all 
available  data  on  the  pholoexcilalion  ol  "''In'"  through 
a  single  intermediate  level  at  1078  kcV  arc  summa¬ 
rized  This  represents  the  most  studied  photocxcitation 
reaction  in  the  literature.  Yet.  the  large  scattering  of 
i loss  sections  even  in  this  particularly  simple  ease  is 
obvious 

Furthermore,  the  method  provides  a  valuable  tool 
for  the  characterization  of  nuclear  simulators  where 
the  extremely  high  dose  outputs  in  a  single  cycle  dis¬ 
able  many  conventional  calibration  methods  [23].  A 
similar  application  for  photon  energies  iimiicd  to  £,  == 
1.5  McV  and  using  photocxcitation  of  the  77Sc  and 
7‘‘Br  isomers  is  described  in  refs.  (24.25). 

2.  Monte  Carlo  simulations 

Simulations  of  thick  target  brcmssirahlung  were 
performed  with  the  MC  codes  GEANT3  and  EGS4. 
These  codes  include  all  electron/ positron  and  photon 
interactions  relevant  at  MeV  energies.  The  experimen¬ 
tal  setup  is  represented  by  building  blocks  of  various 
basic  geometries  which  can  consist  of  different  materi¬ 
als.  Because  of  the  significant  differences  between  the 
results  of  EGS4  and  GEANT3,  version  11  (or  briefly 
GEANT3.il),  calculations  with  the  improved  version 
15  are  also  presented.  The  incorporated  geometry  re¬ 
produced  the  bremsstrahlung  facility  (14)  at  the  injec¬ 
tor  of  the  S-DALINAC  electron  accelerator  in  Darm¬ 
stadt.  It  consisted  of  a  100  p.m  Al  vacuum  exit  window, 
an  air  path  of  2  cm  and  a  3  mm  thick  Ta  radiator  which 
is  capabable  of  a  complete  stopping  of  electrons  with 
energies  up  to  9  MeV.  Based  on  experimental  findings 
a  beam  spot  of  2  jnm  diameter  with  a  Gaussian  distri¬ 
bution  was  assumed  for  the  electrons  incident  on  the 
radiator.  Further  details  of  the  experimental  environ¬ 
ment  can  be  found  in  Fig.  1  of  ref.  [14].  The  spectral 
distributions  were  recorded  an  a  circular  plane  at  a 
distance  of  4.77  cm  at  which  a  radius  of  4  cm  corre¬ 
sponded  to  an  opening  angle  of  80°. 

Calculations  were  performed  for  electron  energies 
E0  -  1,  2,  3,  5  and  7  MeV.  In  most  cases  the  spectra 
correspond  to  107  initial  electrons  which  were  sampled 
with  50  kcV  binsize.  Polynomials  up  to  fourth  order 
were  then  fitted  to  the  histogram  spectra  to  average 
out  statistical  fluctuations  and  show  the  main  features 
more  clearly.  For  the  1  MeV  calculations  the  number 
of  triggers  was  increased  to  4  X  107  in  order  to  keep 
statistical  errors  at  a  level  comparable  to  the  higher 
energies.  Below  an  energy  cutoff  Ec  -  600  keV  the 
computations  were  stopped  for  electrons  and  photons. 
This  threshold  was  chosen  since  EGS4  spectra  for  the 


above  energies  were  ahead;,  available  It  run  the  data 
analysis  ol  the  plutloaeln  ation  expeniuents  |N  12  14  ?n| 
where  the  inhumation  al  lowei  eneigies  w,,s 
needed  Foi  a  test  ol  the  piogiams  down  to  the  low, -a 
energies,  a  seperale  inn  was  peiloinied  with  /  In 
keV  for  an  electron  cneigy  of  1  MeV 

Angular  disli ibntions  ol  the  energy  integialed 
brcmstrahlung  flux  arc  displayed  in  Fig.  I  loi  the 
various  electron  energies  In  this  and  the  following 
figures  F.GS4  results  aie  displayed  as  solid  lines. 
GEANT3.il  as  doited  and  GEANT3  15  as  dashed 
lines  Note  that  the  data  for  very  small  angles  (II  2  i 
are  let t  out.  They  show  huge  fluctuations  because  ol 
poor  statistics  combined  with  a  large  areal  weighting 
factor  in  the  conversion  ol  radial  to  angular  distnhu 
lions  The  differences  obtained  with  GEANT3  11  and 
EGS4  arc  dramatic.  The  GF.ANT3  1 1  angular  distnhu 
lion  is  very  narrow  with  a  maximum  roughly  at  tlu 
critical  angle  6F  =  ntcc: / £„  and  a  minimum  under  tl 
in  contradiction  to  experimental  results  for  thick  target 
bremsstrahlung  (15).  On  the  other  hand  the  improved 
GEANT3.15  version  shows  very  good  agreement  with 
EGS4.  Differences  arc  observed  at  small  angles  only 
where  EGS4  intensities  tend  to  be  somewhat  higher, 
but  typically  less  than  10%  except  for  the  £„  =  1  MeV 
results. 

In  Fig.  2  the  total  intensity  radiated  into  the  80° 
opening  angle  is  depicted  as  a  function  of  electron 
energy  £0.  While  at  the  highest  energy  reasonable 
agreement  of  all  three  calculations  is  obtained, 
GEANT3.il  leads  to  significantly  larger  intensities  with 
increasing  deviations  towards  lower  energies.  The 
GEANT3.15  and  EGS4  curves  are  very  close  except  for 
at  £0=  1  MeV  where  the  former  calculation  is  about 
25%  larger. 

Brcmstrahlung  spectra  for  £0  =  3  MeV  correspond¬ 
ing  to  different  angular  regions  are  displayed  in  Fig.  3. 
Again,  excellent  agreement  of  GEANT3.15  and  EGS4 
is  observed,  in  particular  for  smaller  angles  (upper  two 
figures).  At  larger  angles  (lower  two  figures)  EGS4 
predicts  slightly  higher  intensities  near  the  endpoint 
energy.  The  shape  of  spectra  calculated  with 
GEANT3.il  is  much  too  steep  at  larger  angles.  How¬ 
ever,  the  total  spectrum  is  dominated  by  the  most 
forward  angles  and  in  the  region  around  the  critical 
angle  the  spectral  shape  produced  with  GEANT3.il 
coincides  reasonably  well  with  the  other  results.  At 
other  electron  energies  a  behaviour  similar  to  Fig.  3  is 
observed.  Thus,  despite  the  erroneous  description  of 
the  angular  distributions  and  total  fluxes  the  conclu¬ 
sions  of  refs.  [16,17]  are  not  affected,  since  they  de¬ 
pend  on  only  the  angle  integrated  shape. 

The  results  in  Figs.  1-3  demonstrate  a  very  good 
correspondence  of  the  GEANT3.15  and  EGS4  results. 
Because  of  various  inherent  approximations  in  the 
codes  the  comparison  is  most  critical  at  low  electron 
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Fig.  I.  Angular  distributions  relative  to  the  beam  axis  of  the 
energy-integrated  hiemstrahlung  intensity  from  a  3  mm  Ta 
radiator  and  for  electron  energies  £0  =  1,  2,  3,  5  and  7  MeV 
The  EGS4  results  are  displayed  as  solid  lines,  the  GEANT3.1 1 
results  as  dotted  lines  and  the  GEANT3.15  results  as  dashed 
lines. 

energies.  Fig  4  presents  angle  integrated  spectra  for 
E0  =  1  MeV  calculated  down  to  a  cutoff  £c  =  10  keV. 
Because  of  self  absorption  in  the  thick  target  no  pho¬ 
tons  are  registered  below  200  keV.  The  overall  inten- 


Endpoint  Energy  (MeV) 

Fig.  2.  Total  intensity  of  the  bremsstrahlung  emitted  in  a  80° 
forward  opening  angle  as  a  function  of  electron  energy.  The 
EGS4  results  are  displayed  as  solid  line,  the  GEANT3.il 
results  as  dotted  line  and  the  GEANT3.1S  results  as  dashed 
line. 


Fig  3  Bremsirahlung  spectra  lor  different  angular  regions  at 
3  MeV  The  EGS4  results  are  displayed  as  solid  lines, 
the  GEANT3.il  results  as  dolled  lines  and  the  GIANT’  15 
results  as  dashed  lines 


sity  of  the  GEANT3.I5  spectrum  is  about  .3(1'" "<  larger 
as  already  seen  in  Fig.  2  and  in  the  400-8IK)  keV 
region  the  spectrum  decreases  faster  than  the  EGS4 
results.  Close  to  the  endpoint  energy  a  larger  relative 
yield  is  observed  with  EGS4  in  accordance  with  the 
general  trend  in  all  calculations. 

For  comparison,  the  analytical  approximation  of 
Schiff  (27]  for  thin  target  bremsstrahlung  into  forward 
direction  is  shown  as  dashed-dotted  line.  The  thin 
target  assumption  limits  the  comparison  to  energies 
where  self  absorption  is  negligible  (E  s  300  keV).  Simi¬ 
lar  to  the  findings  of  ref.  (16]  at  higher  electron  ener¬ 
gies  the  MC  spectral  shape  at  intermediate  energies  is 
well  reproduced,  but  the  high-energy  part  is  severly 
overestimated. 


3.  Experimental  bremsstrahlung  calibration  with  the 
n5In(7,Y'),l5In"  reaction 

The  photoexcitation  of  isomer  proceeds  via  reso¬ 
nant  photoabsorption  into  higher-lying  intermediate 
states  (IS)  which  subsequently  decay  with  a  certain 
probability  to  the  metastable  level,  either  directly  or  in 
a  cascade.  In  the  case  of  a  continuous  photon  source 
like  bremsstrahlung,  information  about  the  energy  and 
integrated  cross  sections  of  IS  can  be  obtained  from 
the  isomer  yield  excitation  function  (5,8,11-14,26].  Al¬ 
ternatively,  if  the  relevant  properties  of  the  IS  are 
known  from  other  experiments  the  measured  yield  can 
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Fig  4  Angle-integrated  brcmssltahlung  spectra  for  fc„  =  I 
MeV  The  histogram  binsi/e  is  10  keV  The  cutoff  energy  m 
the  calculations  wav  lowered  to  £,  =  10  keV.  The  EGS-l 
results  are  displayed  as  solid  histogram  and  the  GEANT3  I  s 
results  as  dashed  histogram  The  dashed-dotted  line  corre- 
ponds  to  a  calculation  with  the  Schiff  formula  12K)  fot  a  thin 
target 


be  utilized  to  deduce  the  photon  field  0(0,  £R,  £„) 
for  a  given  endpoint  energy  £0.  Here,  0  is  the  angle 
relative  to  the  incident  electron  beam  and  £R  repre¬ 
sents  the  energy  of  the  resonant  IS.  As  described 
below,  the  latter  method  is  applicable  for  the 
115In(y,y')ll5Innn  reaction  at  photon  energies  up  to 
*  3  MeV. 

An  investigation  of  the  ll5In  isomer  photoexcitation 
was  performed  for  electron  energies  £0  =  2-6  MeV  at 
the  injector  of  the  S-DALINAC  electron  accelerator 
[281  in  Darmstadt.  Disks  of  nominally  1  g  of  natural 
indium  with  a  diameter  of  3.8  cm  were  placed  in  the 
beam  direction  at  a  distance  of  6.43  cm  from  the 
converter.  After  an  irradiation  of  typically  4  h  with  an 
average  electron  ‘current  of  20  jaA  the  characteristic 
isomeric  transition  (£  =  336  keV,  Tt/2“  4.486  h)  of 
1,5  In  was  recorded  with  a  Ge(Li)  detector.  These  val¬ 
ues  were  transformed  into  yields  including  the  correc¬ 
tions  described  in  ref.  [8].  A  detailed  description  of  the 
data  analysis  can  be  found  in  ref.  [21]. 

The  isomeric  yield  Nt  can  be  expressed  as 


Alf*AfT/£V(£) 

Ec 


d 4>(£,  £0) 


where  iVx  denotes  the  number  of  target  atoms,  dd>/d£ 
describes  the  spectral  intensity  per  energy  of  the  pho¬ 
ton  field  for  an  endpoint  energy  £0,  and  <r(£)  is  the 
resonant  absorption  cross  section.  Note  that  the  0 
dependence  of  the  spectral  distribution  is  implicitly 
considered  in  the  definition  of  the  target  geometry  in 
the  MC  calculations. 
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with  i , r  /  T  giving  the  integrated  eioss  seelion  ol  the 
; i h  level  We  note  that  noiucson.ml  eioss  seelions 
wlneh  would  mhihil  ihc  use  of  Eg  (2).  are  nol  consul 
eied  (  I. mils  of  importanl  nonresonanl  conliibutions 
|2')]  io  i lie  phoioexeitation  of  isomeis  have  teeeniK 
been  Jispioven  1 30]  and  shown  to  have  resulted  Ironi 
i he  omission  of  intense  contributions  to  the  photon 
fields  in  those  experiments  due  to  environmental 
Compton  scattering 

The  integrated  cross  seelion  of  an  IS  is  lelated  lo  its 
level  properties  by 


2,l|s  are  ihe  spins  of  ihc  g  s  and  ihc  inlermediaie  stale. 

/  describes  the  total  width  and  b„.  h1M>  denote  ihe 
branching  ratios  for  decay  into  the  g.s.  and  isomeric- 
state.  respectively.  If  the  relevant  level  properties  £, 
b, ,.  61V1.  £  and  J  of  alt  intermediate  levels  arc  known. 
Eq  (2)  can  be  inverted  and  information  on  the  spectral 
intensities  al  the  resonance  energies  can  be  extracted 
For  "'In  the  available  data  base  [31]  provides  a  com 
plcte  characterization  of  all  levels  contributing  to  the 
isomer  yield  up  to  £ ,  =  1.65  MeV.  The  adopted  data 
and  derived  integrated  cross  sections  are  summarized 
in  Table  1  A  theoretical  yield  function  can  be  con¬ 
structed  from  the  MC  photon  spectra  and  tes:  -d  against 
the  measured  results. 

This  method  was  used  in  ref.  [12]  to  determine  the 
contribution  of  low-lying  levels  and  to  extract  informa¬ 
tion  on  previously  unknown  IS  at  higher  excitation 
energies.  Here,  we  are  concerned  with  the  quality  of 
reproduction  of  the  experimental  yield  curve  at  low 


Table  1 

Properties  of  ll5In  levels  contributing  to  the  population  of  the 
0.336  MeV  isomer  for  £,s2  MeV  (from  ref.  (311) 


E. 

(MeV) 

J"  h/2 

(ps) 

*0 

(%) 

(%) 

(oT)^ 

(meV  b) 

g.s. 

9/2*  6x  10M  * 

- 

- 

- 

0.336 

1/2’  4.486  b 

- 

- 

- 

0.934 

7/2*  57(7) 

99.3(3) 

0.5(2) 

0.013302) 

0.941 

5/2*  15.1(1.4) 

88.0(6) 

12.0(16) 

0.83(8) 

1.078 

5/2*  0.99(5) 

83.4(6) 

15.5(4) 

118(6) 

1.449 

9/2*  0.35(4) 

84.7(7) 

0.2(1) 

3.3(4) 

1.463 

7/2*  0.063(10) 

93.5(7) 

0.91(8) 

89(14) 

1.486 

9/2*  0.44(6) 

79.6(6) 

0.45(9) 

6.5(9) 

1.497 

7/2  *  0.167(20)* 

78.1(8) 

21.9(8) 

639(80) 

1.608 

7/2*  0.113(18)* 

69.0(7) 

4.8(2) 

160(26) 

*  Partial  g.s.  width  from  ref.  [32),  branching  ratios  from  ref. 
[31). 
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Endpoint  Energy  (MeV) 

Fig  5  Yield  excitation  function  of  the  ll'ln(y.-y')"'ln"‘  rc.u 
lion  The  solid  line  describes  a  calculation  using  the  energies 
and  integrated  cross  sections  of  Table  1  and  GIANT '  I' 
bremsstrahlung  spectra  The  dashed  lines  give  the  relative 
error  of  the  data  summarized  in  Table  1 

energies  using  the  MC  generated  photon  spectra  Fig 
5  displays  the  data  and  theoretical  results  (solid  line) 
using  the  GEANT3.15  spectra  described  in  section  2. 
The  dashed  lines  give  the  relative  error  of  the  calcula¬ 
tion  introduced  by  the  experimental  data  of  Table  I 
Since  the  lowest  technically  feasible  electron  energy 
was  £0  =  2  MeV  and  the  level  information  above  E,  = 
1.65  MeV  is  incomplete,  some  problem  could  arise 
from  the  presence  of  additional,  yet  unknown  IS  with 
1.65  s£,s  2  MeV.  However,  Fig.  5  clearly  shows  that 
the  yield  variation  for  E0  -  2-2.8  MeV  is  well  de¬ 
scribed  by  the  calculations  including  only  the  IS  of 
Table  1.  We  conclude  that  no  additional  IS  with  a 
significant  contribution  to  the  isomer  yield  is  present 
up  to  2.8  MeV  excitation  energy.  The  clear  break  in 
the  excitation  function  at  2.9  MeV  is  attributed  to  a 
new  IS  (see  ref.  [12]). 

The  other  important  criterion  is  the  reproduction 
on  an  absolute  scale.  If  one  performs  a  least-square  fit 
to  the  experimental  data  points  below  3  MeV,  an 
overall  normalization  factor  F  « 1.27  results  for  the 
GEANT3.15  spectra.  The  corresponding  result  with 
EGS4  spectra  is  F  —  1.13.  Due  to  problems  described 
below  a  systematic  error  of  ±30%  must  be  assumed. 
Within  this  error,  calculations  with  EGS4  and 
GEANT3.15  are  able  to  quantitatively  describe  the 
yields  obtained  for  the  U5In(y,y')u5Inm  reaction. 

The  main  contribution  to  the  error  comes  from  the 
efficiency  calibration  of  the  Ge(Li)  detector.  The  1,5  In 
isomeric  decay  was  measured  in  close  proximity  (the 
disks  were  placed  on  the  detector  cap)  which  leads  to  a 
strong  sensitivity  on  details  of  the  detector  geometry. 
The  energy  dependence  of  the  detection  efficiency  was 


nK.ivuicd  with  .i  siM  ol  poult  I  ike  y  v.ilibi.ilion  mhik.- 
in  i  distance  ol  III  cm  fiom  llic  dcu  cioi  i.tp  Sin. , 
most  -.  soul ccs  could  mil  lv  used  in  Ilk  t  lose  c.nimn  . 
without  severe  de.idtinie  cllccis  we  li.nl  to  iclv  on 
M(  simulation  ol  the  absolute  ellicicmv  o<  i lie  '  ■ 
kc\  isomeric  ti.uisiiion  Also  ilicic  w.is  .m  in.icitu an 
knowledge  ol  the  precise  dclectoi  civstal  si/c  and  . >! 
the  distance  lioni  the  crystal  suilace  to  the  Al  window 
and  some  discrepancies  belween  the  nicasmed  ell; 
cieiicies  al  l.uge  and  small  distance  and  the  calculi 
lions  Therefore,  we  estimate  the  accuracy  to  be  m. 
belter  than  30‘T  Unfortunately,  in  the  meantime  the 
detector  has  snllcrcd  from  severe  neutron  damage  and 
no  repetition  ol  the  efficiency  calibration  i'  possible  to 
resolve  this  problem 

4.  Discussions  and  conclusions 

The  present  work  addresses  two  questions  related 
to  the  problem  of  an  absolute  calibration  ol  thick 
target  bremsstrahlung  for  election  energies  ol  seveial 
MeV  The  calculation  of  the  electron-photon  cascade 
with  the  MC'  codes  EGS4  and  GEANT3  are  compared 
While  the  results  of  EGS4  are  reasonably  verified 
against  experimental  results  1 20).  no  such  tests  have 
been  reported  to  dale  for  this  application  of  GEANT3 
The  Figs.  1-4  demonstrate  the  extreme  differences  of 
GEANT3,  version  11.  compared  to  EGS4  in  all  but 
one  aspects:  the  angle  integrated  shape  of  the 
bremsstrahlung  spectra  agrees  reasonably  well  Thus, 
concerns  about  the  use  of  GEANT3.il  spectra  in  the 
calibration  of  NRF  experiments  [16.17]  are  removed. 

On  the  other  hand,  results  obtained  with  GEANT3. 
version  15,  reveal  very  good  agreement  with  EGS4  over 
most  of  the  investigated  energy  range.  At  energies 
E0>  2  MeV  the  remaining  differences  are  small  enough 
to  permit  interchangeable  use  of  the  results  obtained 
with  both  programs  for  most  practical  purposes.  It 
seems  worthwile  to  extend  the  test  of  GEANT3  to 
higher  energies  used  in  medical  applications,  since  a 
recent  comparison  of  the  10-30  MeV  region  reports 
considerable  differences  between  EGS4  and  other  MC 
codes  [33].  The  differences  observed  at  low  energies 
might  be  partly  due  to  the  fact  that  EGS4  neglects  the 
Coulomb  correction  (Elwert  factor  [34]),  which  in¬ 
creases  the  cross  sections  considerably  at  E0  —  1  MeV, 
while  GEANT3.15  is  based  on  the  recent  parametriza- 
tion  of  bremsstrahlung  cross  sections  of  Berger  and 
Seltzer  [35]  using  exact  partial  wave  calculations  [36]  at 
energies  below  2  MeV. 

An  experimental  method  for  absolute  calibration  of 
thick  target  bremsstrahlung  spectral  distributions  with 
endpoint  energies  up  to  about  3  MeV  is  presented.  It 
relies  on  the  measurement  of  the  excitation  function  of 
the  1I5In(y,y')n5Inm  reaction.  Since  the  properties  of 
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all  IS  in  M'ln  with  ,  <  2.N  Mc\  contributing  io  the 
isdiiii  i  popul.ilion  an-  known  tin  viclil  unit  can  he 
calculated  iiMiit  Ml  lesulte  I  lie  ilisti  ibuiion  ol  IS  al 
scvci.il  cxcit.ilion  energies  is  sensitive  to  the  speettal 
slt.ipes  and  the  i.m.itiuii  with  endpoint  energy  tests  the 
eneigv  dependence  ol  the  total  flux  Also  angular 
distributions  can  he  measured  with  the  In  samples 
placed  off  axis 

The  results  described  lieic  indicate  that  the  M( 
programs  tested  above  are  capable  ol  icpiodi icing  the 
experimental  spectral  shapes  and  absolute  intensities 
in  the  McV  tange  Due  to  the  uncertainties  in  the 
efficiency  determination  of  the  tied  a)  detector  used- 
for  detection  ol  the  "'In  isomeric  transition,  a  reta 
lively  large  systematic  error  of  about  fit'/  remains  for 
the  absolute  normalization  In  the  light  of  possible 
applications  like  the  characterization  of  unknown  pho¬ 
ton  fields  of  new  hremsstrahlung  facilities  ot  nuclear 
simulators  and  the  calibration  of  isomer  excitation  ex¬ 
periments  it  would  be  worthwhile  to  repeat  and  extend 
the  above  described  experiments  in  order  to  determine 
an  absolute  calibration  of  the  MC  generated  spectra 
with  improved  accuracy. 

Nevertheless,  these  results  already  constitute  a 
benchmark  for  the  series  of  recent  photoexcitation 
experiments  [5,8,12-14,26],  The  confirmation  of  abso¬ 
lute  values  of  the  integrated  cross  sections  (if  IS  at 
higher  excitation  energies  extracted  from  the  MC  gen¬ 
erated  photon  spectra  can  be  utilized  for  an  extension 
of  the  calibration  method  towards  higher  photon  ener¬ 
gies.  As  a  result  of  the  improved  nuclear  data  this 
method,  called  “X-ray  activation  of  nuclei"  [23-25,37], 
has  now  be  markedly  improved  as  discussed  in  ref.  [38]. 
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Low-energy  conversion  electron  MOssbauer  spectroscopy  using 
a  chevron  microchannel  plate  detector 
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A  detection  system  for  conversion  electron  Mossbauer  spectroscopy  is  described.  A  chevron 
microchannel  plate  assembly  attached  to  a  two-stage  electrostatic  lens  is  used  to  preferentially 
detect  electrons  with  energies  <13  eV.  Mossbauer  spectra  collected  with  these  electrons  can 
provide  information  about  a  variety  of  solid-state  surface  phenomena. 


I.  INTRODUCTION 

Conversion  electron  Mossbauer  spectroscopy1-4 
(CEMS)  is  a  low-noise  spectroscopic  technique  based  on 
the  detection  of  conversion.  Auger,  and  secondary  elec¬ 
trons  emitted  from  a  material  as  a  result  of  the  de¬ 
excitation  of  Mossbauer  nuclei.  It  provides  information 
about  the  environment  in  which  these  sensitive  nuclei  are 
embedded.  However,  the  range  which  can  be  probed  with 
the  method  is  limited  by  the  escape  depth  of  the  electrons 
in  the  material.  For  57Fe  the  mean  electron  range  in  me¬ 
tallic  iron  has  been  reported  to  be  57  and  36  nm  for  the  7.3 
and  5.4  keV  conversion  and  Auger  electrons,  respectively. 1 
It  has  been  shown  that  the  component  with  energies  <15 
eV  comprises  about  50%  of  the  total  signal  and  provides 
information  from  the  top  5  nm  of  the  sample.*"7  Because 
the  range  of  those  electrons  is  so  short,  they  can  provide 
specific  information  about  a  variety  of  solid-state  surface 
phenomena  and  aid  in  the  characterization  of  implanted 
layers.*-10 

In  this  article  we'discuss  the  design  and  construction  of 
a  relatively  simple  spectrometer  that  was  used  to  selec¬ 
tively  and  efficiently  detect  this  low-energy  component.  It 
was  based  upon  a  chevron  microchannel  plate  (MCP)  as¬ 
sembly  with  an  active  diameter  of  25  mm  that  served  as  the 
detector.  A  two-stage  electrostatic  lens  was  used  to  prefer¬ 
entially  focus  and  accelerate  low-energy  electrons  from  the 
target  to  the  detector.  This  acceleration  increased  the  en¬ 
ergy  of  the  electrons  to  a  level  more  compatible  with  the 
sensitivity  of  the  input  microchannel  plate. 

The  detection  of  higher  energy  signal  electrons,  which 
escape  from  and  provide  information  about  deeper  regions 
of  the  target,  could  serve  as  a  source  of  "contamination”  to 
the  signal  in  that  they  carry  little  or  no  information  about 
the  surface  effects  being  investigated.  The  number  of  these 
electrons  detected  was  reduced  to  negligible  levels  by  min¬ 
imizing  the  geometrical  solid  angle  between  the  Mossbauer 
target  and  the  MCP  assembly.  This  was  possible  since  the 
trajectories  of  the  more  energetic  electrons  were  not  signif¬ 
icantly  focused  by  the  relatively  weak  electrostatic  fields  of 
the  tens. 


II.  EXPERIMENTAL  DETAIL 
A.  Apparatus 

Figure  1  shows  the  experimental  arrangement  used  in 
these  measurements.  A  57Co(Pd)  source  with  a  nominal 
activity  of  10  mCi  was  mounted  on  a  Doppler  motor  ex¬ 
ternal  to  the  vacuum  chamber  containing  the  sample  and 
electron  detector  and  provided  the  14.4-keV  Mossbauer 
radiation.  To  optimize  the  passage  of  resonant  y  rays  into 
the  target  chamber  a  thin  entrance  window  was  con¬ 
structed  from  titanium.  The  14.4-keV  y  rays  were  only 
slightly  attenuated  but  the  window  efficiently  filtered  out  a 
significant  number  of  the  lower  energy  x  rays  that  were 
produced  by  the  decay  of  source  nuclei.  Without  filtering, 
these  x  rays  would  have  seriously  degraded  the  perfor¬ 
mance  of  the  detection  system. 

Low-energy  electrons  emitted  from  the  target  were  col¬ 
lected  with  a  two-stage  aperture  lens.  The  MCP,  target, 
electrical  feedthrough,  and  electrostatic  lens  were  mounted 
as  an  integral  unit  on  a  6-in.  vacuum  flange.  This  provided 
for  convenient  mounting  and  removal  of  the  entire  assem¬ 
bly  from  the  vacuum  chamber  for  adjustments  to  the  target 
geometry.  The  angular  position  of  the  Mossbauer  target 
with  respect  to  the  direction  of  the  source  y  rays  was  ad¬ 
justed  by  rotating  it  about  the  vertical  axis.  Data  were 
typically  collected  with  the  target  at  a  45*  angle  with  re¬ 
spect  to  the  detector  axis  and  the  incident  radiation. 

The  chamber  was  initially  evacuated  with  an  Alcatel 
CFF-450  turbomolecular  pump  (TMP)  to  a  pressure  of 
about  3x  10“ 7  Torr.  However,  vibrations  from  the  TMP 
seriously  degraded  the  resonant  absorption  in  the  absorber 
and  could  not  be  used  to  maintain  the  vacuum  necessary  to 
operate  the  chevron.  To  provide  for  vibrationless  pumping, 
a  140  f/s  Varian  diode  Vaclon  pump  was  used  to  maintain 
the  required  pressure  of  2x  10-6  Torr.  It  is  important  to 
note  that  such  pumps  operate  by  ionizing  gas  in  a  magnet¬ 
ically  confined  cold-cathode  discharge.  Thus  proper  pre¬ 
cautions  had  to  be  taken  to  prevent  neutral  and  charged 
particles,  as  well  as  hard  ultraviolet  radiation  from  escap¬ 
ing  from  the  Vaclon  pump  and  degrading  the  performance 
of  the  detector.  To  avoid  this  difficulty  a  line-of-sight  baffle 
was  used  that  consisted  of  six  successive  aperture  plates. 
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FIG.  1.  Schematic  diagram  of  the  experimental  arrangement 

This  was  inserted  between  the  Vaclon  pump  and  target 
chamber  with  the  apertures  staggered  in  a  helical  pattern. 

B.  Lens  design 

The  structure  and  dimensions  of  the  lens  used  in  this 
work  are  shown  in  Fig.  2.  The  steering  of  the  electrons  was 
performed  primarily  by  the  potential  gradients  established 
between  two  annular  rings  mounted  inside  a  wire  cylinder. 
The  rings  had  a  common  outer  diameter  of  12.7  cm  and 
inner  diameters  of  2.5  and  6.4  cm,  respectively.  The  overall 
length  of  the  lens  was  14  cm.  It  not  only  served  to  focus 
and  accelerate  electrons  to  the  detector,  but  also  as  a  Far- 


FIG.  2.  Drawing  of  the  two-element  electron  lens  designed  to  efficiently 
collect  electrons  having  energies  <  IS  eV.  The  steering  of  the  electrons  was 
performed  primarily  by  the  potential  gradients  established  between  A, 
and  X].  With  the  inclusion  of  the  endcaps  on  the  assembly,  the  lens 
system  not  only  served  to  focus  and  accelerate  electrons  to  the  detector 
input  plate,  but  also  as  a  Faraday  cage  in  shielding  the  lens  elements  from 
the  ground  potential  and  noise  electrons  of  the  chamber  walls.  The  opti¬ 
mum  operating  biases  for  the  lens  were  found  to  be  FX|=-28I7  V, 
- 1967  V,  and  —2757  V  when  the  chevron  was  operated  at 

- 1616  V  and  the  MCP  grid  at  - 1870  V. 


Grid  UCP  1  MCP  2  U«Ui  La od« 

I  \  /  / 


FIG.  3.  Circuit  diagram  of  the  supporting  electronics  for  the  chevron 
microchannel  plate  assembly 


aday  shield  isolating  the  lens  elements  from  the  ground 
potential  of  the  chamber  walls  due  to  the  inclusion  of  end- 
caps  on  the  assembly.  This  minimized  the  number  of  noise 
photoelectrons  produced  inside  the  chamber  that  entered 
the  sensitive  region  of  the  detector. 

To  model  the  lens  and  to  determine  the  optimum  po¬ 
tentials  for  the  various  elements  the  program  EGN2C  was 
employed  which  is  traceable  to  SLAC."  This  program  is 
designed  to  calculate  trajectories  of  charged  particles  in 
electrostatic  and  magnetostatic  fields  and  is  capable  of  in¬ 
cluding  effects  due  to  space-charge  and  self-magnetic  fields. 
In  this  work  these  effects  were  not  investigated.  Input  to 
the  program  included  electrode  boundary  conditions  and 
initial  conditions  for  electrons  emitted  from  the  surface  of 
the  target. 

Various  electron  emission  configurations  including 
plane  waves  and  point  sources  were  simulated  for  electrons 
with  energies  ranging  from  15  eV  to  7.3  keV.  Biases  were 
found  that  achieved  an  effective  solid  angle  for  collection  of 
nearly  2ir  steradians  for  electrons  with  energies  below  15 
eV  and  are  indicated  in  Fig.  2.  A  circuit  diagram  of  the 
supporting  electronics  for  the  MCP  assembly  is  shown  in 
Fig.  3.  Simulations  indicated  that  the  lens  had  a  negligible 
effect  on  the  trajectories  of  electrons  with  energies  greater 
than  about  100  eV  for  the  biases  used.  The  detection  effi¬ 
ciency  for  these  electrons  was  then  solely  defined  by  the 
geometrical  solid  angle  between  the  detector  and  the  tar¬ 
get. 

C.  Noise  reduction 

To  obtain  Mossbauer  spectra  with  large  signal-to-noise 
ratios,  proper  collimation  of  the  source  was  critical  to  re¬ 
duce  nonresonant  scattering  inside  the  chamber.  Various 
collimator  designs  were  used  in  these  experiments  but  the 
best  performance  was  realized  with  a  1.25  cm  thick  lead 
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FIG.  4  Conversion  electron  Mossbauer  spectra  for  a  93.55%  enriched 
'7Fe  target:  (a)  Lens  in  a  focusing  or  active  mode  and  collected  for  3  4  h 
with  a  counting  rate  of  2  47  cps.  (b)  With  the  lens  in  a  passive,  field-free 
mode  possessing  no  focusing  or  accelerating  properties  the  signal  count¬ 
ing  rate  dropped  to  0.087  cps.  The  data  were  collect  five  times  longer  than 
in  (a)  The  geometrical  collection  of  the  lens  accounted  for  less  than 
3  5%  of  the  signal  when  the  lens  was  operated  at  optimum  biases. 


plate  that  had  a  circular  aperture  at  its  center.  The  aper¬ 
ture  was  tapered  so  that  its  diameter  increased  in  the  di¬ 
rection  of  propagation  of  the  radiation  to  reduce  scattering 
from  the  collimator  walls. 

Radiation  scattering  inside  the  chamber  was  mini¬ 
mized  by  providing  an  exit  window  diametrically  opposite 
the  entrance  window.  This  window  had  a  diameter  of  10 
cm  and  was  constructed  of  mylar  to  avoid  significant  back- 
scattering  of  photons  and  photoelectrons.  The  beam  size  of 
Y  rays  from  the  37Co  source  was  optimized  by  proper  col- 
limation  to  illuminate  the  entire  target  without  expanding 
to  a  diameter  that  would  impede  its  exit  from  the  output 
window. 

To  further  reduce  scattering  in  the  system  the  lens  was 
constructed  from  a  copper  mesh  which  had  an  open  area  of 
69%  and  a  wire  diameter  of  250  pm.  Perturbations  of  the 
electric  field  resulting  from  the  use  of  a  mesh  instead  of  a 
continuous  sheet  of  metal  have  been  shown  to  be  negligible 
at  distances  greater  than  the  wire  separation.12  The  per¬ 
centage  of  open  area  was  further  increased  by  alternately 
removing  wires  in  the  region  of  the  lens  that  were  in  the 
direct  path  of  the  input  radiation. 

III.  RESULTS 

The  target  foil  used  to  evaluate  the  effectiveness  of  the 
lens  was  a  5.0X5.0X2X  10-4  cm  iron  foil  with  J7Fe  nuclei 
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enriched  to  93.55%.  For  purposes  of  comparison,  data 
were  collected  with  the  electron  lens  operating  in  an  “ac¬ 
tive”  mode  at  the  optimum  potentials  for  low-energy  elec¬ 
trons  and  in  a  passive  mode  where  there  were  no  acceler¬ 
ating  fields.  Data  collected  with  the  lens  in  the  focusing  or 
active  mode  are  shown  in  Fig.  4(a).  The  total  collection 
time  for  the  data  was  3.44  h  and  the  total  signal  counting 
rate  was  2.47  counts  per  second.  The  data  with  the  lens  in 
the  passive,  field-free  mode  are  shown  in  Fig.  4(b).  There 
the  signal  counting  rate  was  0.087  counts  per  second.  The 
data  were  collected  five  times  longer  than  the  correspond¬ 
ing  data  with  the  lens  active  to  obtain  reasonable  statistics. 

Operation  of  the  lens  in  a  passive  mode  helped  deter¬ 
mine  the  geometrical  collection  efficiency  for  the  system.  In 
this  geometry  it  was  concluded  from  the  data  of  Fig.  4  that 
the  lens  increased  the  signal  count  rate  by  a  factor  of  28 
while  the  background  signal  increased  by  only  a  factor  of 
1 1 .  For  the  biases  used  in  the  active  mode  the  trajectories 
of  the  higher  energy  signal  electrons  were  defined  solely  by 
the  geometrical  solid  angle  between  the  target  and  MCP 
assembly  and  were  the  same  for  both  bias  configurations. 
This  contribution  to  the  signal  counting  rate  could  have 
been  reduced  further  by  simply  increasing  the  distance  be¬ 
tween  the  target  and  detector  assembly. 

IV.  DISCUSSION 

Detection  of  electrons  with  energies  <15  eV  can  pro¬ 
vide  detailed  information  about  a  variety  of  surface  phe¬ 
nomena  in  the  top  5  nm  of  a  material.  In  this  article  we 
have  described  the  design  and  construction  of  a  detection 
system  for  conversion  electron  Mossbauer  spectroscopy. 
This  device  preferentially  detects  electrons  with  energies 
<  15  eV  and  has  the  potential  for  providing  a  simple  alter¬ 
native  to  more  complex  and  expensive  electron  detectors 
which  are  currently  used  in  energy-resolved  Mossbauer 
spectroscopy. 

It  was  shown  that  the  detection  efficiency  of  higher 
energy  conversion  and  Auger  electrons  can  be  minimized 
by  choosing  an  appropriate  geometry.  Numerous  computer 
simulations  using  plane-wave  and  point-source  electron 
configurations  showed  that  the  electrostatic  lens  gave  an 
effective  solid  angle  of  almost  2 rr  steradians  for  the  signal 
electrons  having  appropriately  low  energies. 

Although  not  pursued  in  this  work,  it  is  possible  with 
the  current  design  to  collect  virtually  all  the  electrons  emit¬ 
ted  from  the  target.  This  could  be  accomplished  if  the 
potential  difference  between  the  input  of  the  chevron  and 
target  foil  is  made  sufficiently  large.  Furthermore,  if  a  re¬ 
tarding  field  analyzer  is  positioned  between  the  target  and 
MCP  the  lower  energy  electrons  can  be  efficiently  filtered 
out.13  In  this  manner  it  is  possible  to  obtain  energy- 
resolved  Mossbauer  spectra  over  a  wide  range  of  electron 
energies. 
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ABSTRACT 

We  determine  the  intensity  of  the  nuclear  multiphoton 
sidebands  to  the  forbidden  components  of  the  14.4  keV  Mossbauer 
transition  in  57Fe,  generated  by  the  applicatior  f  radio- 
frequency  magnetic  fields  to  ferromagnetic  sample.  in  the 
presence  of  static  magnetic  fields.  Further  we  determine  the 
intensity  of  the  second-order  sidebands  to  the  14.4  keV 
transition  in  the  case  known  as  the  radio-frequency  collapse, 
when  the  entire  hyperfine  magnetic  field  acting  on  the  57Fe 
nuclei  follows  the  oscillations  of  the  applied  radio-frequency 
field.  Then  we  discuss  a  number  of  characteristic  asymmetries 
due  to  quantum  interference,  appearing  in  the  multiphoton  radio¬ 
frequency  sideband  spectra. 
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I .  INTRODUCTION 


We  consider  in  this  work  a  class  of  nuclear  multiphoton 
processes  involving  the  simultaneous  interaction  with  a  nucleus 
of  a  gamma-ray  photon  and  of  radio-frequency  photons  from  an 
applied  oscillatory  magnetic  field,  the  intermediate  states  being 
the  magnetic  sublevels  of  the  nuclear  states.  The  ratio  of  the 
intensity  of  a  two-photon  transition  to  the  intensity  of  the 
single  gamma-ray  transition  between  the  same  nuclear  states  is  of 
the  order  of  the  square  of  the  ratio  of  the  interaction  energy  of 
the  radio-frequency  magnetic  field  with  the  nucleus  and  the 
frequency  of  this  field. 

The  first  experimental  demonstration  of  a  nuclear  transition 
involving  simultaneously  a  gamma-ray  photon  and  a  radio-frequency 
photon  is  due  to  West  and  Matthias  in  1978, 1  for  the  6.2  keV 
state  of  181Ta  under  an  applied  radio-frequency  magnetic  field  of 
23  mT  having  a  frequency  of  0.94  MHz.  A  similar  nuclear  two- 
photon  process  has  been  demonstrated  in  1988  by  Ikonen  et  al.2-4 
for  the  93  keV  state  of  67Zn  under  an  applied  radio-frequency 
magnetic  field  of  13.4  mT  at  the  frequencies  3.1  kHz  and  9.3  kHz. 

When  we  consider  the  possibility  to  observe  such 
simultaneous  gamma  radio-frequency  transitions  in  57Fe,  we  have 
to  take  into  consideration  the  fact  that  the  half-life  of  the 
14.4  keV  transition  in  57Fe  is  98  ns,  compared  to  a  6.8  ns  half- 
life  for  the  6.2  keV  transition  in  181Ta,  and  a  9.2  ns  half-life 
for  the  93  keV  transition  in  67Zn.  Moreover,  the  spins  and  the 
nuclear  magnetic  moments  of  the  ground  and  the  excited  states  are 
respectively  1/2,  0.090/in  and  3/2,  -0.155/<n  for  57Fe,  compared  to 


7/2,  2 . 37/ijj  and  9/2,  5.3 for  181Ta,  and  5/2,  0.87/in  and  1/2, 
0.58/in  for  67Zn,  where  nn  -  5.05  x  10“27  J/T.  Thus,  the  energy 
of  interaction  of  an  applied  radio-frequency  field  with  a  57Fe 
nucleus  is  considerably  smaller  than  the  coupling  of  the  same 
field  to  181Ta  or  67Zn  nuclei.  Therefore,  a  two-photon 
experiment  for  57Fe,  without  any  enhancement  mechanism,  would 
require  applied  magnetic  fields  with  an  amplitude  of  the  order  of 
0.1  T  at  a  frequency  of  a  few  MHz,  which  is  difficult  to  produce. 

However,  if  a  57Fe  nucleus  is  in  a  ferromagnetic 
environment,  the  direction  of  the  hyperfine  magnetic  field  acting 
on  the  nucleus  is  determined  by  the  direction  of  the 
magnetization  in  the  sample,  and  the  orientation  of  the  bulk 
magnetization  can  be  controlled  with  the  aid  of  applied  static 
and  radio-frequency  magnetic  fields.  This  approach  has  been  used 
by  Heiman  et  al . 5-6  who  observed  the  spectrum  of  the  gamma  rays 
scattered  by  finely  powdered  metallic  iron  under  the  action  of  an 
oscillatory  magnetic  field  whose  frequency  of  26  MHz  was  equal  to 
the  hyperfine  splitting  of  the  upper,  14.4  keV  state  of  57Fe. 
The  incident  gamma  rays  were  themselves  resonant  to  one  of  the 
allowed  transitions  between  the  magnetic  nuclear  sublevels,  and 
the  arrangement  has  been  referred  to  by  the  authors  as  an 
excited-state  nuclear  magnetic  resonance  experiment.  As  a  result 
of  the  application  of  the  radio-frequency  magnetic  field  to  the 
iron-powder  scatterer,  lines  appeared  in  the  spectrum  of  the 
scattered  gamma  radiation  at  the  position  of  the  other  known 
transitions  of  57Fe. 


On  the  theoretical  side,  the  effects  of  circularly-polarized 
radio-frequency  magnetic  fields  on  the  Mdssbauer  spectra  of  57 Fe 
have  been  analyzed  by  Gabriel.7 

Initially  Heiman  et  al.  have  tried  to  use  a  metallic  iron 
foil  for  the  scatterer,8  but  they  have  found  that  the  application 
of  the  radio-frequency  magnetic  field  produced  in  this  case 
additional  lines  at  positions  which  did  not  correspond  to  any 
known  transition  in  ^7pe  _  They  have  suggested9  that  the 
observed  sidebands,  occurring  in  a  foil  but  with  negligible 
intensity  in  a  powder,  result  from  the  acoustic  vibrations 
generated  in  the  foil  by  the  applied  radio-frequency  magnetic 
field.  Similar  observations  have  been  reported  by  Asti  et 

al.10,11 

The  fact  that  the  hyperfine  magnetic  field  acting  on  a  57Fe 
nucleus  in  a  ferromagnetic  environment  can  be  made  to  oscillate 
between  opposing  directions  by  the  application  of  radio-frequency 
fields  has  been  demonstrated  by  the  radio-frequency  collapse 
effect  discovered  by  Pfeiffer  in  1971. 12 • 13  He  found  that,  as  a 
result  of  the  application  of  an  oscillatory  magnetic  field  with 
an  amplitude  of  1.5  mT  and  frequencies  of  39  MHz,  61  MHz  and  106 
MHz  to  a  6-pm-thick  permalloy  foil,  the  six-line  hyperfine 
pattern  of  57Fe  disappeared  and  has  been  replaced  by  a  single 
central  line,  accompanied  by  a  number  of  sidebands  displaced  on  a 
frequency  scale  by  multiples  of  the  applied  frequency. 

Pfeiffer  attributed  the  presence  of  the  sidebands  in  his 
collapse  experiments  entirely  to  the  acoustic  vibrations 
generated  in  the  permalloy  foil  by  the  radio-frequency  field. 


However,  Olariu  et  al.14  pointed  out  that  the  oscillation  through 
large  angles  of  the  hyperfine  magnetic  field  acting  on  the  57Fe 
nuclei,  produced  by  the  application  of  the  radio-frequency  field 
to  the  permalloy  sample,  results  necessarily  in  the  generation  of 
nuclear  multiphoton  sidebands  whose  transition  intensities  are 
comparable  to  those  observed  by  Pfeiffer.  The  fact  that  the 
entire  hyperfine  magnetic  field  can  be  made  to  follow  the 
direction  of  the  applied  radio-frequency  field,  and  the  reduced 
amplitude  of  the  mechanical  vibrations  in  the  permalloy  sample, 
due  to  the  small  magnetostrictive  constant  of  permalloy,  show 
that  it  is  possible  to  observe  nuclear  multiphoton  transitions 
for  57Fe  nuclei  in  a  ferromagnetic  environment. 

The  radio-frequency  collapse  effect  has  been  further 
investigated  by  Kopcewicz  et  al.  in  ferromagnetic  amorphous 
metals, 15 • 16  in  FeB03,17  and  in  Fe-Ni  alloys,18  with  the 
conclusion  that  for  sufficiently  intense  radio-frequency  fields, 
of  the  order  of  1.5  mT,  the  hyperfine  magnetic  field  acting  on 
the  57Fe  nuclei  follows  the  direction  of  the  applied  field.  The 
sidebands  observed  in  these  experiment  have  been  in  general 
attributed  to  acoustic  modulation. 

As  a  corollary  of  their  studies  on  the  magnetic  modulation 
of  the  phase  of  67Zn  nuclei,  Ikonen  et  al.2  have  reiterated  the 
previous  prediction14  that  if  the  hyperfine  field  acting  on  57Fe 
nucleus  can  be  switched  between  two  opposite  directions  at  a  rate 
determined  by  the  frequency  of  the  applied  radio-frequency  field, 
this  will  generate  sidebands  without  any  mechanical  motion. 


The  intensities  of  the  nuclear  multiphoton  sidebands  for  the 
six  allowed  transitions  of  57Fe  have  been  calculated  by  Olariu19 
for  57 Fe  nuclei  embedded  in  thin  ferromagnetic  films  under  static 
and  radio-frequency  magnetic  fields.  He  remarked  that  the  use  of 
a  thin  ferromagnetic  film  having  a  thickness  of  a  few  thousand 
angstroms,  instead  of  a  foil,  is  likely  to  minimize  the 
mechanical  vibrations,  and  moreover  would  allow  the  determination 
with  reasonable  accuracy  of  the  time  evolution  of  the 
magnetization  of  the  sample. 

The  role  played  by  the  spin  waves  in  the  generation  of  the 
radio-frequency  sidebands  has  been  emphasized  by  Sinor  et  al.20 
who  showed  that  it  is  possible  to  produce  a  modulation  of  the 
phase  of  57 Fe  nuclei  in  paramagnetic  media  by  transporting  spin 
waves  of  large  amplitude  from  ferromagnetic  sources. 

The  resonance  line  splitting  of  hyperfine  lines  of  57Fe 
nuclei  in  radio-frequency  magnetic  field,  caused  by  Rabi 
oscillations,  predicted  in  a  number  of  papers, 7 • 19 > 21  has  been 
recently  observed  by  Vagizov22  and  independently  by  Tittonen  et 
al.23  Special  precautions  have  been  taken  in  these  works  to 
minimize  the  mechanical  vibrations. 

It  is  apparent  from  this  survey  of  previous  worH,  on  radio¬ 
frequency  sidebands  that,  while  the  ferromagnetic  foil  samples 
can  provide  the  required  large  radio-frequency  fields,  in  the 
form  of  oscillatory  hyperfine  magnetic  fields,  the  observation  of 
the  nuclear  multiphoton  processes  is  often  hindered  by  the  ease 
with  which  mechanical  vibrations  can  be  induced  in  ferromagnetic 
samples  through  magnetostriction,  especially  at  lower 


frequencies .  The  importance  of  nuclear  multiphoton  processes 
becomes  apparent  when  we  refer  to  recent  developments  in  atomic 
multiphoton  processes,  where  problems  like  inversionless 
amplification24'25  or  electromagnetically-induced 
transparency26-28  are  under  active  consideration. 

The  purpose  of  this  work  is  to  brinq  into  evidence  a  number 
of  characteristic  properties  of  nuclear  multiphoton  processes 
involving  gamma-ray  and  radio-frequency  photons.  Thus,  in 
Section  II  we  determine  the  intensity  of  the  nuclear  multiphoton 
sidebands  to  the  forbidden  components  of  the  14.4  keV  transition 
of  57Fe,  generated  by  the  application  of  radio-frequency  magnetic 
fields  to  ferromagnetic  samples,  in  the  presence  of  a  static 
field.  In  Section  III  we  determine  the  intensity  of  the  second- 
order  sidebands  to  the  14.4.  keV  transition  of  57Fe  in  the  case 
of  the  radio-frequency  collapse,  when  the  entire  hyperfine 
magnetic  field  follows  the  oscillations  of  the  applied  radio¬ 
frequency  field.  In  Section  IV  we  discuss  a  number  of 
characteristic  asymmetries  due  to  quantum  interference,  appearing 
in  the  multiphoton  radio-frequency  sideband  spectra. 

II.  Intensity  of  the  first-order  sidebands  to  the  forbidden 

e,  3/2  ->  g,  -1/2  and  e,  -3/2  ->  g,  1/2  transitions  in  57Fe 

The  intensities  of  the  first-order  sidebands  to  the  six 
allowed  transitions  between  the  14.4  keV  excited  state  e  and  the 
ground  state  g  of  57 Fe  nuclei  have  been  determined  in  Ref.  19  as 
a  function  of  the  direction  and  polarization  of  the  gamma  rays 
for  a  linear  and  a  circular  polarization  of  the  applied  radio- 


frequency  magnetic  field.  The  transitions  e,  3/2  g,  -1/2  and 
e,  -3/2  ->  g,  1/2  are  forbidden  as  single  gamma-ray  processes, 
but  these  forbidden  transitions  have  in  general  non-vanishing 
multiphoton  sidebands,  the  intensity  of  which  we  shall  calculate 
in  this  section. 

The  intensity  of  the  sidebands  to  the  forbidden  hyperfine 
transition  of  57Fe  has  been  reported  previously29  in  the 
particular  case  of  a  linearly-polarized  radio-frequency  mangetic 
field  i?*,  and  for  a  direction  of  propagation  of  the  gamma  rays 
which  was  perpendicular  to  the  static  and  radio-frequency 
magnetic  fields.  In  the  previous  paper29  it  has  been  also 
pointed  out  that  sidebands  to  the  forbidden  hyperfine  transition 
of  57Fe  can  be  identified  in  the  M&ssbauer  spectra  of  57 Fe  nuclei 
in  permalloy  recently  reported  by  Tittonen  et  al.23  The  presence 
of  the  sidebands  to  the  forbidden  hyperfine  transitions  of  57Fe 
constitute  the  unequivocal  proof  that  nuclear  multiphoton 
processes  are  inded  taking  place  in  57Fe  nuclei  subject  to  radio¬ 
frequency  magnetic  fields,  as  predicted  repeatedly.2*7/ 34 '19» 39 

We  shall  determine  the  probability  for  the  transition 
e,  3/2  -»  g,  -1/2  in  a  57Fe  nucleus  by  the  emission  of  a  linearly 

polarized  gamma-ray  quantum,  the  orientation  of  the  magnetic 

— ¥  — ► 

gamma-ray  field  Br  and  of  the  wave  vector  kr  being  shown  in 

Fig.  l,  under  the  action  of  a  static  hyperfine  magnetic  field 
Bq  ,  oriented  along  the  z  axis,  and  of  a  radio-frequency  magnetic 
field  having  the  x  component  Br ffX  cos (a^  ft)  and  the  y  component 
Brf,y  cosfcurft  +  arf ) .  The  positive  direction  of  the  z  axis  is 


chosen  to  coincide  with  the  direction  of  the  static  hyperfine 
magnetic  field  effectively  acting  on  the  57Fe  nuclei,  so  that  the 
magnetic  sublevels,  in  order  of  increasing  energy,  are 
g,  1/2;  g,  -1/2;  and  e,  -3/2;  e,  -1/2;  e,  1/2;  e,  3/2. 

It  can  be  shown  that  the  transition  amplitudes  b %  are 
solutions  of  the  Hamiltonian  equations 

=  (1) 
n 

The  non- vanishing  matrix  elements  are  given  by 
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where  A  is  a  reduced  matrix  element ,  and 

Z  =  sin  <p  sin  L>-cos^cos#coso+/(cos<psin  u+sin  pcosflcosu)- 


(5) 


In  Eqs.  (2)  and  (4),  o)e,  o)„  >  0  represent  the  separations,  on  a 
frequency  scale,  of  the  magnetic  sublevels  of  the  excited  state  e 
and  respectively  of  the  ground  state  g  due  to  the  static 
hyperfine  field  5*,  and  T  =  1/r,  where  the  mean  life  is  r  =  1.4  x 
10”7  s  for  the  14.4  keV  level  of  57Fe.  Ee  and  Eg  are  the 
energies  of  the  excited  and  ground  states  in  the  absence  of  any 
applied  magnetic  field. 

The  Eq.  (1)  can  be  solved  by  a  perturbation  series,  and  the 
second-order  contribution,  which  gives  the  amplitudes  of  the 
first-order  sidebands  to  the  gamma-ray  transitions,  is 
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Substituting  in  Eg.  (6)  the  expression  of  the  appropriate  matrix 
elements  from  Eqs.  (2)  -  (5) ,  and  neglecting  non-resonant  terms, 
we  obtain  an  expression  for  the  intensity  of  the  first-order 
sidebands  for  a  linear  polarization  of  the  radio-frequency  field, 
when  BT£'y  -  0,  that  is  given  by 


|(g =  3^2p2 (sm2  u+cosJ  0c°f  u)^r) 
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where  the  alternate  signs  should  be  read  along  the  upper  line  or 
the  lower  line. 

If  the  radio-frequency  field  has  a  right  circular 
polarization,  so  that  Brf,x  ~  Brt ,y  ~  Brf  and  arf  ~  ~nf2> 
intensity  of  the  sidebands,  obtained  from  Eq.  (6),  is 
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If  the  radio-frequency  field  has  a  left  circular 
polarization,  so  the  Br f/X  «  Brf,y  *  Brf  an<*  arf  “  w/2,  the 
intensity  of  the  sidebands,  obtained  from  Eq.  (6),  is 


\(g-i\e  =M^-(siii2  v+coi  tfcos2  )  , 

\\S,  *  2  f  +  |  %h2T2  K  0)a-ajJ 

Et-E  3a)t  (og 

for  <or= — ^-+®rf  » 


Mm>!“T-o  * 

(9b) 

^  £.-£i,3  <Q.  <Og 

for  • 

It  can  be  shown  that  the  intensities  of  the  first-order 
sidebands  to  the  forbidden  transition  e,  -3/2  ->  gr,  1/2  are 
connected  to  the  previous  intensities  by  the  relations 


lu.ik.-cf = . 


(10a) 


(10b) 


for  all  polarizations  of  the  applied  radio-frequency  field,  where 
the  superscript  (2)  refers  to  the  second-order  contribution  in 
the  perturbation  series. 


As  was  the  case  with  the  sidebands  of  the  allowed 
transitions, 19  the  intensity  of  the  sidebands  towards  higher 


frequencies  of  a  parent  line  is  in  general  different  from  the 
intensity  of  the  sideband  towards  lower  frequencies  of  the  same 
parent  line.  This  asymmetry  becomes  particularly  striking  in  the 
case  of  a  circular  polarization  of  the  radio-frequency  field, 
when  not  only  the  parent  forbidden  line  has  a  zero  intensity,  but 
the  intensity  of  one  of  the  sidebands  is  also  equal  to  zero. 

By  combining  the  results  of  Ref.  19  and  the  present  work,  we 
obtain  the  intensities  of  the  sidebands  towards  higher 
frequencies,  labelled  by  the  subscript  +,  which  are  shown  in  Fig 
2  for  linear,  right  circular  and  left  circular  polarizations  of 
the  radio-frequency  magnetic  field.  The  intensity  of  the 
unperturbed  transition  e,  1/2  ->  g,  If  2,  averaged  over  the 
polarization  of  the  gamma  rays,  and  observed  at  the  angle  0  = 
7i/2,  was  chosen  as  unity.  The  remarkable  aspect  of  the  pattern 
in  Fig.  2  is  the  appearance  of  the  supplementary  lines  which  are 
the  sidebands  of  the  forbidden  parent  transitions.  The 
appearance  of  these  lines  in  a  spectrum  is  characteristic  of  the 
multiphoton  generation  of  the  radio-frequency  sidebands. 

If  the  57Fe  nucleus  is  under  the  action  of  a  circularly 
polarized  magnetic  field  but  the  static  field  Bj  is  equal  to 
zero,  the  intensity  of  the  first-order  sidebands,  at 
o)r  =  (£.-£,)/#  ±  <urf,  can  be  shown  to  be 
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where 


The  ratio  of  the  intensity  of  the  first-order  sideband  to  the 
parent  line,  at  the  angle  of  observation  0-0,  is  then,  for  a 
circular  polarization  of  the  radio-frequency  field. 


+  3/+/!). 


(12) 


The  results  in  Eqs.  (11)  and  (12)  differ  from  the  corresponding 
expressions,  Eqs.  (32)  and  (33)  of  Ref.  19,  by  the  contribution 
of  the  sidebands  to  the  forbidden  lines,  evaluated  in  this 
section. 


III.  Intensity  of  the  second-order  sidebands  in  the  case  of 
radio- frequency  collapse  for  57Fe 

While  the  second-order  perturbation  calculations  are  in 
general  sufficient  to  describe  the  generation  of  radio-frequency 
sidebands  at  frequencies  of  tiie  order  of  100  MHz  or  higher,  the 
third-order  contribution  to  the  perturbation  series  is  needed  to 
describe  the  second-order  sidebands  which  become  observable  at 
lower  frequencies.  In  this  section  we  shall  determine  the 
intensity  of  the  second-order  sidebands  for  the  case  of  the 
radio-frequency  collapse,  when  the  hyper fine  magnetic  field 
acting  on  the  57Fe  nuclei  follows  the  oscillations  of  the  applied 
radio- frequency  field,  so  that  the  static  hyperfine  field  is 


The  third-order  contribution  to  the  transition  amplitude  is 


given  by 

•  (13> 

"  *  0  0  0 

We  shall  assume,  as  in  the  previous  section,  that  the  applied 
radio-frequency  field  is  in  the  x,y  plane,  and  shall  evaluate  the 
right-hand  side  of  Eq.  (13)  for  gamma  rays  emitted 
perpendicularly  to  the  x,y  plane,  so  that  we  shall  consider  that 
9-0,  <p  =  o. 

The  non-vanishing  matrix  elements  for  the  present  problem 

are 


Vtw=-^(o, 
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where  Y (t)  is  the  function  defined  in  Eq.  (3),  and 
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where  <o0  =(e, -Eg)/h  .  Since  according  to  Eqs.  (14)  and  (15)  the 

matrix  elements  are  non-vanishing  only  for  Am  -  ±  1,  there  is  a 
non-vanishing  third-order  contribution  to  the  amplitude  b(3) gja;em 
for  the  14.4  keV  transition  only  when  Im-mgl  *  1,  so  that  the 
transitions  to  be  considered  are  e,  3/2  ->  g,  1/2;  e,  1/2  ->  q, 
-1/2;  e,  -1/2  — >  q,  1/2  and  e,  -3/2  ->  g,  -1/2.  From  the 
symmetry  of  the  matrix  elements  in  Egs.  (14)  and  (15)  it  results 
that  it  is  sufficient  to  evaluate  b(3)^  i/2;e,3/2  and 
Jb  ( 3 )  gf  -\/2;e,  1/2  at  a>y=a)(yt2<D rf.  There  are  four  channels 

contributing  to  g#  1/2  ;e,  3/2»  and  there  are  six  channels 

contributing  to  *^3^g,-i/2;e, 1/2*  Thus,  we  obtain  for  the 
intensity  of  the  second-order  sidebands,  at  &^-<2>0±2a>rf , 
averaged  over  the  polarization  u  of  the  gamma  rays, 

h  =  [fi^5*‘(19+30/+20/I+6/s+/<)+2[(^51-)J+(^;)I](3+3/+/1)]>  (16) 


where 


B,  =  BAjc  +iBttyeia« 


B-i  -  Brf  x  +  iBH  ye  ** 
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For  a  linear  polarization  of  the  radio-frequency  field,  when 
Brf,y  ~  °>  the  ratio  of  the  second-order  sideband,  averaged  over 


gamma-ray  polarization,  to  the  unperturbed  partent  intensity  is, 
for  <Oy  *  <dq  ±  2<urf , 

^>=t^(tS:),(31+4J/+24/'+6/,+/<)>  (l8> 

where  f  =  -HgJe/neJg.  This  expression  is  equivalent  to  the 
result  reported  previously  in  Eg.  (17)  of  Ref.  14. 

If  the  radio-frequency  field  acting  on  the  nucleus  can  be 
described  as  a  superposition  of  two  independent  orthogonal  modes, 
so  that  Br fiX  »  Brf,y  ~  Brf  and  the  relative  phase  ar f  is 
random,  then  the  ratio  of  the  intensity  of  the  second-order 
sideband,  averaged  over  gamma-ray  polarization,  to  the  parent 
intensity,  is 


=  V31+42/+24/2+6/j+/4)  .  (19) 

256  \J§n 

If  the  radio-frequency  field  is  circularly  polarized,  so 
that  BrffX  =  Brffy  =  Br f  and  l«rfl  *  n/2,  the  ratio  of  the 
intensity  of  the  second-order  sideband,  averaged  over  the 
polarization  of  the  gamma  rays,  to  the  unperturbed  parent 
intensity  is 


j*c)  _  J_  M.Bh 
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If,  as  a  result  of  the  application  of  an  oscillatory 
magnetic  field  of  frequency  a^f,  the  hyperfine  field  assumes 


alternatively  two  opposite  values  ±  Bs,  the  Fourier  component  of 
this  hyperfiene  field  at  the  frequency  will  have  the 
enhanced  amplitude  (4 /n)Bs.  The  amplitude  of  the  component  of 
the  afore-mentioned  hyper  fine  field  at  the  frequency  2a)rf  is 
equal  to  zero,  but  the  amplitude  of  the  component  at  the 
frequency  3a>rf  will  be  4Ss/3ic,  and  can  produce,  at  low 
frequencies,  a  sizable  sideband.  These  considerations  lose  their 
applicability  if  the  passages  between  the  opposing  values  of  the 
hyper fine  field  are  rounded  off. 

Combining  the  results  of  Ref.  14  and  the  results  of  the 
present  work,  we  have  represented  in  Fig.  3  the  emission  gamma- 
ray  spectrum  with  sidebands  in  the  case  of  radio-frequency 
collapse  for  57 Fe  nuclei  in  a  permalloy  sample,  for  a  linear 
polarization  of  the  applied  oscillatory  magnetic  field  at  several 
frequencies,  assuming  that  the  amplitude  of  the  oscillating 
hyperfine  field  is  equal  to  its  static  value.  Moreover,  in  Fig. 
4  we  have  represented  the  emission  gamma -ray  spectrum  with 
sidebands  at  a  fixed  frequency  for  several  configurations  of  the 
applied  radio-frequency  field.  The  intensities  in  parts  (a)  and 
(b)  of  Fig.  4  are  calculated  according  to  the  Bessel-f unction 
formula  of  Ref.  14,  and  the  intensities  in  part  (c)  and  (d)  are 
calculated  according  to  the  perturbation  expressions  of  Ref.  14 
and  the  present  work.  It  is  apparent  from  Figs.  3  and  4  that  if 
the  hyperfine  field  follows  the  applied  radio-frequency,  as  it  is 
supposed  to  be  the  case  for  the  radio-frequency  collapse  spectra, 
the  generation  of  nuclear  multiphoton  sidebands  can  be  very 
intense . 


IV.  Quantum  interference  effects  in  the  radio-frequency 
multiphoton  spectra  of  57Fe 

While  the  intensity  of  the  acoustic  sidebands  of  a  given 
order  is  a  decreasing  function  of  frequency,  the  intensity  of  the 
multphoton  sidebands  of  a  given  order,  in  the  presence  of  static 
magnetic  fields,  can  increase  in  certain  frequency  ranges.  The 
increase  with  frequency  of  the  intensity  of  the  multiphoton 
transition  can  arise  from  the  fact  that  some  of  the  transition 
amplitudes  have  resonances  when  the  applied  frequency  becomes 
equal  to  the  separation  in  frequency  units  between  the  magnetic 
sublevels  of  the  ground  or  excited  states,  i.e.  when  a  real 
intermediate  state  in  the  nucleus  coincides  with  the  position  of 
the  virtual  multiphoton  level. 

This  pattern  is  encountered  for  example  to  the  first-order 
sideband  towards  lower  frequencies  of  the  transition  e, 
1/2  -»  g,  1/2,  in  the  presence  of  a  static  magnetic  field  5*  and 
of  a  linearly  polarized  radio-frequency  magnetic  field.  If  the 
emitted  gamma  rays  are  observed  along  a  direction  perpendicular 
to  the  plane  of  the  applied  fields,  0  **  n/2,  <p  *  n/2,  then  the 
ratio  of  the  first-order  sideband  intensity  /ff  to  the  intensity 

of  the  unperturbed  parent  transition  Ib,  both  averaged  over  the 
gamma-ray  polarization,  becomes19 


iff 

fM2 

3  (0.  ,  1  <0.  1  > 

h 

1*  J 

8ffl.  +  (yrf  4a>,-<yrf  8 
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(21) 


The  resonances  of  the  relative  sideband  intensity  /lb , 
occurring  at  0% f  -  <oe  and  a^f  *  a>g,  are  shown  in  Fig.  5  for  57Fe 
nuclei  in  iron. 

Another  circumstance  which  can  lead  to  an  increase  with 
frequency  of  the  intensity  of  the  sidebands  is  the  cancellation 
for  certain  frequencies  of  the  partial  amplitudes  contributing  to 
the  total  transition  amplitude.  This  pattern  is  encountered  for 
example  to  the  first-order  sidebands  towards  lower  frequencies  of 
the  transition  e,  1/2  ->  g,  1/2  in  the  presence  of  a  static 
magnetic  field  and  of  a  left  circular-polarized  radio¬ 
frequency  field  of  rotating  amplitude  £rf.  The  ratio  of  the 
first-order  sideband  intensity  lf&'}  to  the  intensity  Ia  of  the 

unperturbed  transition  e,  3/2  ->  g,  1/2  is19 

/g1  !|%W  3*>. 

1. 

for  all  gamma-ray  polarizations  and  directions  of  observation. 
The  zero  in  the  relative  sideband  intensity  /la ,  occurring  at 

(Qj- f  =  2<oeOg/(3©e  +  (Og) ,  can  be  seen  in  Fig.  6  for  57Fe  nuclei  in 
an  iron  environment,  when  the  zero  occurs  at  o> rf  =  19.1  MHz, 
assuming  that  £rf  is  small  compared  to  B* . 

An  interesting  example  of  cancellation  of  amplitudes  is 
encountered  for  the  first-order  sidebands  of  the  forbidden  gamma- 
ray  transition  e,  3/2  ->  g,  -*  1/2.  In  this  case  the 
cancellation  of  amplitudes  occurs  at  a^f  -  0,  so  that  the 
intensity  of  the  first-order  sidebands  at  o^f  *  0  is,  according 


to  Eqs.  (7)  -  (9),  equal  to  zero  for  all  polarizations  of  the 
applied  radio-frequency  field  while,  on  the  other  hand,  the 
generation  of  acoustic  sidebands  is  supposed  to  increase  in 
intensity  at  lower  frequencies.  For  example,  the  ratio  of  the 

first-order  sideband  intensity  |(g,-X|e.K)^)|  to  the  intensity  Ia 

of  the  unperturbed  transition  e,  3/2  ->  g,  1/2,  in  the  presence 
of  a  static  magnetic  field  B *  and  for  a  linear  polarization  of 
the  ratio-frequency  field,  is,  according  to  Eq.  (7) , 


for  all  gamma-ray  polarization  and  directions  of  observation. 
The  intensity  of  these  sidebands  are  represented  in  Fig.  7  for 
57Fe  nuclei  in  iron. 

V.  CONCLUSIONS 

In  this  paper  we  have  demonstrated  a  number  of  asymmetries 
characteristic  of  the  nuclear  multiphoton  generation  of  radio¬ 
frequency  sidebands  to  gamma-ray  transitions  in  57Fe  nuclei.  We 
have  pointed  out  that  these  are  sidebands  of  observable  intensity 
to  the  forbidden  14.4  keV  gamma-ray  transitions  e,  3/2  — >  g,  -1/2 
and  e,  -3/2  ->  g,  1/2,  and  have  calculated  the  intensities  of  the 
first-order  sidebands  for  a  linear  and  a  circular  polarization  of 
the  applied  radio-frequency  field.  Then  we  have  determined  the 
intensity  of  the  second-order  sidebands  in  the  case  of  the  radio¬ 
frequency  collapse,  and  have  evaluated  the  sideband  intensities 
for  a  number  of  configurations  of  the  applied  magnetic  field. 


Finally,  we  have  remarked  that  unlike  the  decreasing  patterns  of 
the  intensity  of  the  acoustic  sidebands,  the  intensity  of  the 


nuclear  multiphoton  sidebands  can  increase  with  frequency 
certain  frequency  ranges,  this  being  a  consequence  of 
resonances  extant  in  the  multiphoton  transition  amplitudes, 
also  of  the  cancellation  for  certain  frequencies  of 
interfering  amplitudes. 
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Figure  Captions 


Fig.  1.  Orientation  of  the  magnetic  field  Br  and  of  the  wave 
vector  kr  of  a  quantum  of  gamma  radiation  emitted  by  a  57 Fe 
nucleus. 

Fig.  2  Intensities  of  the  sidebands  towards  higher  frequencies 
for  the  gamma-ray  emission  spectrum  of  57Fe  nuclei  under  a  static 
magnetic  field  B *  and  of  an  oscillatory  magnetic  field  of 
frequency  Jcoe  =  4,  averaged  over  the  polarization  of  the 

gamma  rays  for  different  configurations  of  the  radio-  frequency 
fields.  (a)  The  radio-frequency  field  is  linearly  polarized, 
Brf.r/B*  =  1/10 ,  and  the  direction  of  observation  is  perpendicular  to 

the  plane  of  the  applied  field,  6  »  ic/2,  <p  -  n/2.  (b)  The 
radio-frequency  field  has  a  right  circular  polarization,  ar f 
=  n/2,  and  B^/B^- 1.  The  spectrum  is  observed  along  the 
direction  of  the  static  field,  0  «  0.  (c)  The  radio-frequency 
field  has  a  left  circular  polarization,  ar f  *  k/2,  and  B^/B^  =  \. 
The  spectrum  is  observed  along  the  direction  of  the  static  field, 
6  =  0. 

Fig.  3.  Gamma-ray  emission  spectrum  with  sidebands  in  the  case 
of  radio-frequency  collapse  for  57Fe  nuclei  in  a  permalloy  sample 
for  a  linear  polarization  of  the  applied  magnetic  field  for 
several  frequencies  a^f  of  the  oscillatory  field.  It  is  assumed 
that  the  amplitude  of  the  oscillating  hyperfine  field  is  equal  to 
the  static  hyperfine  field.  The  direction  of  observation  is 


perpendicular  to  a  plane  containing  the  oscillatory  field.  These 
intensities  are  given  by 

where  J#  is  the  Bessel  function  of  order  N. 

Fig.  4.  Gamma-ray  emission  spectrum  with  sidebands  in  the  case 
of  radio-fequency  collapse  for  57Fe  nuclei  in  a  permalloy  sample 
for  different  configurations  of  the  applied  radio- frequency  field 
at  the  frequency  a^f  -  61  MHz.  (a)  The  hyperfine  oscillatory 
field  is  polarized  linearly  and  its  amplitude  is  equal  to  the 
static  hyperfine  field  Bs.  (b)  The  hyperfine  field  is  polarized 
linearly  and  its  amplitude  is  4Bs/ic.  (c)  The  hyperfine  field  is 
a  superposition  of  two  orthogonal  components  of  amplitude  BSl 
uncorrelated  in  phase.  (d)  The  hyperfine  field  is  polarized 
circularly,  and  the  rotating  amplitude  is  equal  to  Bs.  In  all 
cases,  the  direction  of  observation  is  perpendicular  to  a  plane 
containing  the  oscillatory  field. 

Fig.  5.  Relative  first-order  sideband  intensity  lf£} /lb  for  a 

static  magnetic  field  and  a  linearly  polarized  magnetic  field 

of  amplitude  Br ftX,  for  57Fe  nuclei  in  iron.  The  relative 
intensity  has  resonances  at  a^f  ■  ©e,  a^f  *  ojg.  The 

direction  of  observation  is  perpendicular  to  the  plane  of  the 
applied  field.  The  gamma-ray  intensities  are  averaged  over 
polarization. 


a 


Fig.  6.  Relative  first-order  sideband  intensity  /*^')//- ,  for 
static  magnetic  field  B*  and  a  left  circular  polarized  radio¬ 
frequency  field  Br f,  for  57Fe  nuclei  in  iron.  The  relative 
intensity  /fc(CL)//a  has  a  zero  when  -  2a>e<0g/  (2<oe  +  <Og) ,  or 

=  19.1  MHz.  Br f  is  assumed  to  be  small  compared  to  B* . 

Fig.  7.  Relative  first-order  sideband  intensities 

|(g,-Kk,K)(+CL)|  /la  (full  line)  and  |<^,-X|e,^)!a')|2//.  (dotted  line),  for 

a  static  magnetic  field  B 0*  and  a  linearly  polarized  radio¬ 
frequency  field  Br fiX  for  57Fe  nuclei  in  iron.  The  intensity  of 
these  sidebands  are  equal  to  zero  for  *  0. 
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ABSTRACT 

The  analysis  of  overall  self-consistent  scenarios  for  the  process  of 
stimulated  gamma-emission  is  now  required  because  of  the  nuclear  pumping 
experiments  demonstrated  at  The  University  of  Texas  at  Dallas.  Two  separate 
schemes  will  be  discussed  here:  (1)  The  stimulation  of  radiative  transitions  in 
beams  of  free  isomeric  nuclei  without  overall  population  inversion.  The  central 
concept  is  the  "cooling"  down  of  particles  in  the  beams  to  the  submicroKelvin 
range.  (2)  The  stimulated  annihilation  of  electrons  and  positrons  in  colliding 
particle  and  photon  beams  by  external  low-frequency  radiation.  Numerical 
estimates  for  both  schemes  are  presented. 


1.  INTRODUCTION 

The  study  of  the  problems  inherent  in  achieving  induced  emission  of 
gamma-rays  began  more  than  three  decades  ago  and  interest  has  frequently 
peaked  in  the  intervening  years.  The  very  first  maximum  of  effort  arose  in  the 
early  sixties  when  preliminary  treatments  were  made  independently  in  the  USA 
(1-2)  and  in  Russia  >3-51. 


The  second  peak  arose  in  1974.  In  Russia  this  peak  was  initiated  by  the 
efforts  of  a  brilliant  scientist,  Professor  Rem  Khokhlov,  the  President  of  the 
Moscow  State  University.  Unfortunately,  this  promising  period  ended  in  1977, 
when  Khokhlov,  an  enthusiastic  mountain-climber,  perished  in  the  Pamirs 
range. 


Today  we  can  record  a  new  peak  of  interest,  perhaps  the  most 
significant  in  terms  of  experimental  investigations.  This  period  has  been 
marked  by  the  world's  first  experimental  approach  to  one  of  the  most  important 
part  of  the  problem,  the  pumping  of  populations  of  nuclear  isomers.  This  has 
been  undertaken  by  the  group  of  Professor  C.  B.  Collins  from  the  Center  for 
Quantum  Electronics  of  The  University  of  Texas  at  Dallas,  and  in  collaboration 
with  some  other  Institutes. 

Now,  additional  emphasis  is  needed  upon  other  parts  of  the  problem  of 
induced  gamma-ray  emission  that  are  worthy  of  examination.  Of  particular 
importance  at  this  time  are  those  which  are  aimed  at  - 

*  extending  fundamental  understanding  of  the  stimulated  boson 
emission,  successfully  applied  in  the  optical  laser  physics,  to  a 
new  class  of  quantum  oscillators  -  nuclei  and  antiparticles; 
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*  discovering  ways  of  utilizing  in  modern  science  a  new  energy 
range  of  coherent  photons  -  keV  and  perhaps  even  MeV; 

*  introducing  into  theoretical  and  experimental  practice  a  new  type 
of  nuclear  reactions  -  stimulated  radiative  transitions.  In  particular 
this  would  include  stimulated  radiative  chain  reactions,  or  in  other 
words,  gamma-ray  lasing. 

The  realization  of  these  ideas  may  also  lead  to  important  progress  in 
theoretical  work  that  could  launch  a  new  branch  of  science  and  technology, 
QUANTUM  NUCLEONICS,  as  an  analog  and  continuation  of  the  development  of 
Quantum  Electronics. 

The  problems  in  achieving  this  challenging  goal  are  very  diverse  and 
complex.  The  first  successful  experimental  results  mentioned  above  have 
motivated  the  analysis  of  different  self-consistent  scenarios  for  the 
pumping/stimulated  emission  process  as  a  whole.  In  this  context  the 
discussion  presented  here  will  often  take  more  the  character  of  a  task 
statement  giving  a  direction  for  further  investigations  than  that  of  a  final 
solution  for  the  whole  problem.  However,  it  seems  that  these  efforts  may  be 
truly  well-timed.  Here  two  separate  approaches  to  the  problem  will  be 
discussed:  stimulated  radiative  transitions  in  populations  of  isomeric  nuclei 
without  overall  inversion,  and  stimulated  annihilation  of  electrons  and  positrons 
by  external  low-frequency  radiation. 


2.  CONCEPTS 

2.1  Stimulated  Gamma  Transitions 

Most  proposals  for  a  feasible  system  in  which  to  observe  stimulated 
gamma-ray  emissions  have  as  an  essential  condition  the  necessity  for 
Mossbauer  (recoilless)  transitions  for  nuclei  situated  in  a  solid  matrix.  At  this 
point  it  is  best  to  review  some  well-known  facts  about  such  transitions.  When 
the  temperature  of  the  crystal  containing  the  active  nuclei  is  high  compared  to 
the  Debye  temperature,  the  gamma-transition  line  is  spectrally  split  into  two 
lines.  These  are  the  emission  line  for  excited  nuclei  and  the  absorption  line  for 
unexcited  ones.  Both  lines  are  inhomogeneously  broadened  by  thermal  motion 
and  other  factors.  If  the  line  splitting,  equal  to  twice  the  nuclear  recoil  energy, 
exceeds  the  inhomogeneous  broadening  a  spectrally  local  nuclear  population 
inversion  may  occur,  even  without  the  presence  of  an  overall  inversion  of 
excited  nuclei.  Therefore,  gamma-ray  gain  is  positive  but  small  because  the 
ratio  of  the  homogeneous  (natural)  radiative  line  width  to  the  inhomogeneous 
width  is  very  small.  In  the  low  temperature  case  the  line  splitting  does  not 
exist  and  recoilless  Mossbauer  emission  and  absorption  can  take  place.  Then 
the  line  width  ratio  is  equal  to  unity  but  an  overall  inversion  of  excited  nuclei 
must  be  present  to  achieve  gain.  Thus  the  gain  is  again  small  when  the 
inversion  is  not  sufficiently  large  or  even  negative  when  an  inversion  is  absent. 

Unfortunately,  the  fault  of  the  solid  state  approach  is  the  formation  of 
some  kind  of  vicious  circle.  In  order  to  achieve  an  inversion  the  lifetime  of  the 
excited  nuclei  must  be  long  to  reduce  the  requirements  for  intense  pumping. 
However,  an  unbroadened  Mossbauer  line  of  natural  width  can  not  be  realized 
for  isomers  with  too  long  lifetimes  (i.  e.,  for  a  lifetime  exceeding  one 


microsecond)  because  of  crystal  structure  inhomogeneity.  Because  of 
additional  damage  to  the  lattice  this  problem  grows  catastrophically  under  the 
influence  of  intense  pumping  as  well  as  in  other  circumstances. 

The  result  has  been  that  all  known  attempts  to  use  the  solid  state 
approach  have  so  far  proven  unsuccessful  At  this  time  there  are  wide 
disagreements  about  the  further  benefits  of  this  approach.  The  question 
remains  whether  it  is  possible  make  use  of  the  advantages  of  the  physical 
picture  mentioned  above  while  avoiding  its  negative  features? 

A  very  simple  alternative  approach  involves  a  departure  from  the 
necessity  for  a  solid  matrix.  What  will  happen  when  an  ensemble  of  free 
excited  nuclei  is  cooled?  The  recoil  effects  together  with  the  line  splitting 
remain  and  so  it  is  again  possible  to  obtain  a  local  inversion  over  parts  of  the 
spectral  line  without  an  overall  inversion.  Simultaneously  the  line  width 
decreases  drawing  nearer  to  the  natural  homogeneous  value  and  the  gain 
grows. 

Thus  one  can  imagine  a  possible  scenario  based  on  the  use  of  deeply 
cooled  atomic  or  ion  beams*7'81.  The  realization  of  this  scenario  will  require  the 
following  sequence  of  operations: 

1)  pumping  of  nuclei  (i.e.,  by  (y,y‘)  reactions  or  by  radiative  neutron 
capture)  and  preparing  the  preliminary  cooled  ion  or  atom 
ensemble; 

2)  forming  a  nuclear  beam  and  rapid  supplying  them  from  a  "hot" 
pumping  zone; 

3)  deep  "cooling,"  that  is  to  say,  longitudinal  monokinetization  of  the 
nuclear  beam; 

4)  beam  compression  (if  needed); 

5)  stimulated  gamma-ray  emission  in  the  prepared  "cold  and  dense" 
nuclear  beam. 

In  addition  all  the  steps  must  be  achieved  consistently  with  each  other. 

On  the  final  step  the  threshold  condition  for  amplification  for  gamma 
transitions  is  reached  when  the  rate  of  stimulated  emission  exceeds  that  for 
losses  due  to  the  photoeffect  and  Compton  scattering.  The  temperature 
dependence  of  this  threshold  value  shows  that  deep  "cooling"  or 
monokinetization  down  to  the  submicroKelvin  range  is  needed.  In  addition  the 
maximum  value  of  the  total  gain  over  the  optimum  length  of  a  beam  of 
spontaneously  relaxing  nuclei  is  inversely  proportional  to  the  Doppler  width. 
For  this  reason  deep  "cooling"  is  also  needed  to  increase  the  total  gain. 
Therefore  the  beam  "cooling"  is  the  central  point  of  this  whole  scenario. 

To  achieve  this  may  mean  the  use  of  known  and  promising  methods  of 
laser  light  pressure  that  can  provide  an  effective  temperature  in  the  nanoKelvin 
range,  but  these  need  to  be  improved  in  order  to  accommodate  an  enormous 
rise  in  the  particle  concentration.  In  further  studies  new  methods  to  increase 


the  density  of  deeply  cooled  beams  should  be  proposed  and  the  limiting  factors 
should  be  analyzed. 

It  may  be  useful  to  also  consider  some  combined  methods  of  particle 
monokinetization  involving  an  electrical  acceleration  of  ions  accompanied  by 
optical  laser  velocity  selection.  Optical  pressure  is  not  used  here,  but  the  laser 
light  plays  the  role  of  a  Maxwell's  demon  by  accurately  determining  the 
properties  of  particles,  whereas  the  electric  field  executes  the  task  set  by  the 
demon.  The  idea  is  as  follows*91. 

Negative  ions  are  accelerated  up  to  some  threshold  velocity  which 
considering  the  Doppler  effect  is  sufficient  for  photodetachment  of  electron  by 
a  collinear  optical  laser  beam.  The  resultant  monokinetized  neutral  atoms 
continue  to  move  without  acceleration  at  the  same  constant  velocities  equal  to 
a  threshold  value. 

A  reversed  variant  of  this  method*101  involves  the  compensation  of  the 
initial  random  spread  of  the  particle  velocities  by  controlling  the  duration  of  the 
positive  ion  acceleration  in  a  static  electric  field.  Doppler  velocity 
discrimination  is  performed  by  frequency  scanning  of  a  laser  ionizing  the  parent 
atom  beam. 

Another  group  of  very  simple  methods*7-81  is  derived  from  the  wide 
experience  available  for  radiofrequency  amplifiers  with  low-noise  electron 
beams.  When  ions  are  simultaneously  accelerated  by  an  electrical  field,  their 
longitudinal  random  velocity  spread  drops  drastically  due  to  the  quadratic 
dependence  of  the  ion  kinetic  energy  on  the  velocity.  For  instance,  a  voltage 
of  100  kV  reduces  the  effective  longitudinal  temperature  from  80  Kelvin  down 
to  1  microKelvin.  Consequently  the  Doppler  line  width  also  decreases  as 
considered  from  the  longitudinal  direction.  Obviously  the  degree  of  beam 
"cooling"  as  a  result  of  acceleration  is  not  unlimited.  The  limitations  should  be 
studied  further  in  the  light  of  known  theories  of  plasma  instabilities  and  noise  in 
electron  beams. 

The  gain  in  induced  gamma  photons  is  proportional  to  the  concentration 
of  excited  nuclei.  So  the  problem  of  the  density  enhancement  arises  in  parallel 
to  the  "cooling"  problem.  Leaving  aside  the  known  and  useful  methods  of 
electron  optics  and  laser  assisted  density  enhancement  let  us  discuss  another 
method,  the  ballistic  focusing  of  neutral  particles.  (8>  To  some  extent  this 
combines  both  the  other  approaches.  A  technique  of  electron  optics  is  used 
for  setting  the  initial  conditions  for  the  parent  negative  ions.  Then  laser 
photodetachment  forms  the  ensemble  of  neutral  atoms  with  almost  the  same 
unperturbed  initial  conditions  and  leads  to  further  focusing  unhindered  by 
Coulomb  repulsion. 

A  simple  example,  but  by  no  means  the  best,  of  a  way  of  implementing 
this  method  may  involve  focusing  a  fragment  of  a  negative  ion  beam  having  an 
annular  cross  section  within  coaxial  electrodes  of  circular  cylinder  geometry *8). 
A  pulsed  radial  electric  field  imparts  to  all  the  ions  an  axial  velocity  component. 
A  laser  beam  then  performs  photodetachment  and  motion  of  the  neutral  atoms 
toward  the  axis  takes  place  without  any  influence  from  Coulomb  repulsion. 
This  ballistic  motion  is  accompanied  by  an  increase  in  the  beam  density.  The 
highest  degree  of  density  enhancement  is  achieved  in  the  region  of  the  beam 


waist  and  depends  on  a  small  transverse  random  velocity  components. 
Estimates  show  an  expected  density  enhancement  of  1000. 

To  sum  up  this  analysis  let  us  consider  a  numerical  example  for  a 
hypothetical  isomer  with  a  nuclear  charge  of  20,  a  transition  energy  of  2.5 
keV,  and  a  lifetime  of  1  microsecond.  Longitudinal  monokinetization  by 
electrical  acceleration  up  to  the  ion  energy  of  100  keV  reduces  the  effective 
temperature  from  4  Kelvin  down  to  2.5  nanoKelvin  and  the  Doppler  velocities 
down  to  1  mm/s.  Positive  gain  over  an  optimum  beam  length  of  180  cm  arises 
when  the  parent  ion  beam  current  density  exceeds  1  A/cm2  (taking  into 
account  a  density  enhancement  factor  of  1500). 

This  result  is  promising,  but  the  concepts  need  further  development 
along  lines  including: 

*  searches  for  possible  isomers  with  suitable  nuclear  and  atomic 
properties; 

*  development  of  efficient  pumping  schemes  and  efficient  sources  of 
metastable  nuclei; 

*  improvement  of  methods  of  particle  monokinetization  and  density 
enhancement  relevant  to  the  present  problem  and  critical  analysis 
of  physical  and  technological  limitations. 


✓  2.2  Stimulated  Annihilation  Radiation 

Antimatter  is  the  perfect  source  of  states  with  negative  temperature*11*. 
Nevertheless  there  are  grounds  for  skepticism  as  to  the  feasibility  of  a  simple, 
self-maintaining  photon  chain  reaction,  in  other  words  direct  annihilation  lasing. 
It  is  more  realistic  to  consider  a  method  for  initiating  electron-positron 
annihilation  in  which  the  emission  (or  absorption)  of  one  of  the  photons  in  a 
multiquantum  event  was  produced  from  an  external  electromagnetic  field,  in 
particular  an  intense  optical  laser  field*12'13*.  This  concept  in  the  three-photon 
case  is  based  on  a  continuous  redistribution  of  the  energy  of  the  antiparticles 
between  the  emitted  quanta,  so  that  there  is  a  finite  probability  for  the  emission 
of  photons  of  any  specific  low  frequency*121.  However,  soft-photon  stimulation 
of  two-photon  annihilation  is  possible  only  when  the  corresponding  motions  of 
particles  and  stimulating  photons  are  such  that  the  energy  of  Doppler  shifted 
photons  becomes  equal  to  the  particle  rest  mass  energy'*3*.  In  the  case  of 
optical  stimulation  this  requires  an  inadmissibly  high  particle  energy 
(approximately  equal  to  50  GeV). 

The  new  scenario  for  inducing  the  annihilation  of  parapositronium  in 
colliding  photon  and  particle  beams  (namely,  two  positronium  beams  and  two 
infrared  photon  beams  moving  in  opposite  directions)  consists  of  three  steps: 

1)  exciting  1s-2p  transitions  of  moving  positronium  atoms  by  an 
infrared  beam  propagating  in  the  opposite  direction; 

2)  lasing  in  the  beam  of  these  excited  atoms  (2p-1s  transition),  the 
inversion  being  achieved  as  a  result  of  a  reduction  of  lower-state 
population  due  to  spontaneous  annihilation;*14*  and 


3)  stimulating  the  annihilation  in  the  opposite  positronium  beam  by 
Doppler  shifted  photons  from  the  moving  atomic  laser  realized  in 
the  previous  step. 

The  Lorentz  transformations  on  each  step  leads  to  the  resonant 
condition,  determining  the  particle  energy  of  80  MeV  and  the  wavelength  of 
the  external  radiation  source  of  80  micrometers.  The  resulting  stimulated 
annihilation  radiation  produces  characteristic  photon  energies  of  160  MeV  in 
the  forward  direction  and  1.6  keV  in  the  backward  direction  in  the  laboratory 
frame.  It  should  be  noted  that  in  this  process  the  population  inversion  of  the 
antiparticles  is  not  necessary  and  the  competitive  processes  of  pair  creation  are 
absent,  but  to  achieve  a  good  experimental  contrast  of  stimulated  radiation 
against  a  spontaneous  background  a  very  high  density  of  the  positronium 
beams  is  needed  (an  equivalent  current  density  of  more  than  10  kA/cm2). 

This  work  was  supported  in  part  by  the  Russian  Foundation  of 
Fundamental  Research  under  grant  #  93-02-3768. 
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Abstract 

A  gamma-ray  laser  would  provide  stimulated  emission  at  the  shortest  possible  wavelengths  (<  1  A)  from  transitions  between  excited 
nuclear  states.  The  unique  advantages  provided  by  those  transitions  promise  the  integration  of  pump  energy  over  nanosecond  to  microsecond 
times  and  thus  to  larger  values  than  are  possible  with  x-ray  lasers  which  must  be  pumped  on  the  order  of  picoseconds  1  lowcver,  at  the  end  of 
the  first  cycle  of  study  (1963  -  1981)  it  was  concluded  that  a  gamma-ray  laser  was  impossible  using  the  traditional  approach  to  pumping  nuclei 
with  particle  fluxes.  Since  then  research  has  taken  an  entirely  new  character  with  the  introduction  of  concepts  for  "optical"  pumping  of  a 
gamma-ray  laser,  in  its  most  straightforward  form  as  a  nuclear  analog  to  the  ruby  laser.  Over  the  past  seven  years  the  United  Stales’  main-line 
project  at  the  Center  for  Quantum  Electronics  has  experimentally  demonstrated  the  components  of  those  concepts  Computations  based  on  the 
new  data  have  shown  that  a  gamma-ray  laser  is  feasible  if  some  isotope  has  its  properties  sufficiently  close  to  the  ideal  Now  breakthroughs 
realized  in  the  study  of  the  systemalics  for  the  optical  pumping  of  nuclei  have  permitted  the  identification  of  the  first-ranked  candidate  nuclide 
for  a  gamma-ray  laser.  Problems  currently  being  addressed  arc  the  acquisition  of  macroscopic  samples  for  testing  and  the  demonstration  of 
appropriate  instrumentation. 


Introduction 

Studies  of  the  feasibility  of  a  gamma-ray  laser  have  been  motivated  by  the  promise  of  achieving  stimulated  emission  of  photons  with 
wavelengths  below  1  A  using  the  unique  advantages  presented  by  nuclear  systems.  Among  these  is  the  appearance  in  many  nuclei  of 
extremely  long-lived  excited  states  called  isomers  with  energies  above  the  ground  state  of  keV  to  MeV.  Populations  of  these  metastable 
levels  can  store  teraJoules  per  liter  at  solid  densities  for  lifetimes  up  to  thousands  of  years.  This  makes  upconversion  possible  since  it  would 
not  be  necessary  to  provide  all  the  energy  for  a  gamma-ray  laser  in  situ.  Even  more  important,  electromagnetic  transitions  in  nuclei  provide  a 
variety  of  transition  moments,  the  most  common  being  the  magnetic  dipole,  in  contrast  to  electric  dipole  transitions  which  dominate  our  usual 
experience  with  lasing  transitions.  Thus  the  matrix  element  linking  spontaneous  lifetime  with  the  usual  v3  factor  is  orders -of-magnitude  more 
favorable  for  output  transitions  of  a  gamma-ray  laser  than  might  otherwise  be  expected  at  such  short  wavelengths.  Populations  can  be 
transferred  from  isomeric  storage  states  to  lasing  levels  having  lifetimes  ranging  from  nanoseconds  to  microseconds,  integrating  the  energy  in 
a  pump  pulse  for  much  longer  times  than  for  x-ray  lasers.  Also,  gamma-ray  transitions  routinely  have  natural  line  widths,  as  in  the  MOssbauer 
effect,  insuring  large  cross  sections  for  stimulated  emission.  Without  broadening,  values  for  1  A  nuclear  transitions  typically  exceed  the  cross 
section  for  the  stimulation  of  Nd  in  YAG. 

Despite  this  impetus,  there  was  little  real  progress  towards  a  gamma-ray  laser  prior  to  1982  as  a  result  of  a  concentration  of  research  on  the 
traditional  approach  to  pumping  nuclei  with  intense  particle  fluxes.  As  discussed  in  the  encyclopedic  review  of  Baldwin,  et  al.,1  the  general 
impossibility  of  those  schemes  was  concluded  since  particle  fluxes  would  deposit  unmanageable  amounts  of  heat  in  the  host  matrix, 
destroying  the  MOssbauer  effect  That  review  effectively  signaled  the  end  to  the  traditional  approach  to  a  gamma-ray  laser  and  emphasized 
that  the  critical  test  for  any  new  approach  must  be  the  efficacy  of  pumping  to  release  the  energy  stored  in  isomers  while  maintaining  a  natural 
line  width  for  the  lasing  transition. 

A  new  interdisciplinary  approach2  to  a  gamma-ray  laser  emerged  in  1982  following  precursive  work2'9  that  appeared  near  the  end  of  the 
first  cycle  of  research  but  which  were  largely  ignored.  The  basic  concepts  called  for  the  "optical*  pumping  of  nuclei  by  the  absorption  of 
photons  from  intense  coherent  or  incoherent  fields.  Coherent  optical  pumping  would  use  multiphoton  processes  to  release  the  energy  stored  in 
nuclear  isomers  into  freely-radiating  levels  located  nearby  in  energy.  Incoherent  pumping  would  employ  intense  pulses  of  x  rays  to  pump 
nuclei  as  in  a  direct  analog  to  the  ruby  laser.  In  either  case,  the  use  of  photons  as  the  pump  entity  would  allow  the  preservation  of  the 
MOssbauer  effect  for  nuclei  within  a  thin  film  host  of  low-Z  material  like  carbon  or  beryllium. 

The  blueprint  Cor  a  gamma-ray  laser  now  appears  in  the  textbooks,10  and  research  following  this  approach  at  the  Center  for  Quantum 
Electronics  has  been  rich  with  achievement.  The  fundamental  concepts  of  coherent  pumping  have  been  demonstrated11*13  but  not  surprisingly 
the  greatest  rate  of  progress  has  been  realized  in  the  direction  of  incoherent  pumping.14-25  Major  breakthroughs  have  been  experimentally 
demonstrated  in  the  pumping  and  dumping  of  populations  of  nuclear  isomers  with  x  rays. 

The  greatest  impediment  to  a  gamma-ray  laser  has  remained  the  lack  of  materials  Of  the  1886  distinguishable26*27  isotopes,  29  had  been 
identified  as  first-class  candidates  but  only  (he  two  poorest,  lt0Ta  and  ^Te  were  available  in  sufficient  amounts  for  testing.  Still  the 
discovery  of  resonances  with  giant  cross  sections  for  optically  pumping  those  nuclides1*-20  strengthened  the  feasibility  of  a  gamma-ray  laser  by 
ordera-of-magnitude.  Without  further  materials  for  study,  research  has  concentrated  on  the  systematica  for  pumping  simulation  isomers  from 
their  ground  states31*35  to  aid  in  ranking  the  candidates  and  to  motivate  production  of  the  best  nuclides.  The  findings  have  now  made  it 
possible  to  identify  the  31-year  isomer  of  the  nuclide  17*Hf  as  the  first-ranked  candidate  for  a  gamma-ray  laser. 
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Background 


Pump  Model 

lhc  ciUicj!  key  to  llic  development  of  the  ruby  laser  was  the  identification  and  exploitation  of  a  bandwidth  funnel  in  which  a  broad 
absorption  band  was  linked  lltiough  efficient  cascading  to  a  narrow  laser  level  Likewise  this  is  an  essential  element  m  the  nuclear  analog  to 
the  ruby  laser  depicted  in  Fig  1  Hie  principal  figure  of  merit  for  populating  an  upper-laser  (or  fluorescence)  level  is  the  integrated  cto . . 
section,  (of)  =  7tbtboaor/2,  whether  the  uulial  stale  is  the  ground  stale  or  an  isomer.  The  figure  shows  only  one  pump  band,  sometimes 
called  a  gateway  or  intermediate  state,  but  there  could  be  more  at  higher  excitation  energies.  The  probabilities  that  the  gateway  wall  decay  in 
one  step  back  to  the  initial  level  and  in  any  cascade  to  the  laser  level  arc  bt  and  b0,  respectively  These  quantities  together  with  the  natural 
width  of  the  absorption  band,  T  give  the  useful,  partial  width  btboT  for  the  transfer  of  population  from  initial  to  fluorescence  stales  The  cross 
section  for  the  absorption  step,  aQ  is  given  by  the  Breit-Wigner  formula 
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where  X  is  the  wavelength  of  the  pump  photon,  Ij  and  Ig  are  the  angular  momenta  of  the  intermediate  and  ground  states,  respectively,  and  o(  is 
the  internal  conversion  coefficient  corresponding  to  the  pump  transition.  The  cross  section  for  useful  absorption  of  pump  energy  by  nuclei  is 
much  larger  for  photons  than  for  particles  of  the  same  energy. 

In  the  context  of  nuclear  physics,  optical  pumping  processes  like  that  of  Fig.  1  have  been  known  for  more  than  SO  years  28  They  arc 
commonly  designated  as  (y,/)  reactions  with  y  and  y*  referring  to  the  incident  and  scattered  photons,  respectively.  Typical  values  of  integrated 
cross  section  for  these  reactions  have  been  known  to  range14-17-25  from  I's  to  10's  in  the  "usual  units"  adopted  here,  uu.  =  10-29  cm2  keV 
However,  until  the  past  decade  studies  of  these  processes  never  reached  the  level  of  precision  seen  in  research  using  particle  reactions  due  to 
the  difficulty  in  accurately  calibrating  the  output  spectra  of  the  intense  photon  sources  required.  In  particular  the  most  difficult  to  characterize 
sources  of  x  rays  proved  to  be  radionuclides  for  which  all  intensity  away  from  the  spiectral  lines  resulted  from  environmentally  sensitive 
Compton  scattering.  Little  quantitative  agreement  could  be  found  in  the  literature  between  integrated  cross  sections  obtained  in  different 
experiments  and  as  late  as  1987  there  were  serious  contentions30  over  the  mechanism  by  which  these  reactions  occurred.  However,  in  recent 
years  extensive  experiments  based  on  new  technology  and  computer  codes  for  the  simulation  of  x-ray  spectra  have  resolved  this  controversy 
and  the  basic  model  for  the  optical  pumping  reaction  shown  in  Fig.  1  is  now  on  a  firm  quantitative  basis.  This  was  one  of  the  earliest 
challenges  to  the  gamma-ray  laser  project. 

When  a  population  of  nuclei,  N;  in  some  initial  state  is  irradiated  with  a  continuum  of  photons  extending  up  to  a  maximum  "endpoint" 
energy,  E„,  the  lime-integrated  yield  of  final-state  nuclei,  Nf  is  given  by 

Nf  =  //,0>0  Jo°  <*(E)  F(E,E0)dE  .  (2) 

where  G>0  is  the  total  incident  fluence  and  F(E^,j)  is  the  distribution  of  energies  within  the  continuum  for  the  given  endpoint.  The  initial-to- 
final  state  transfer  is  represented  by  the  energy-dependent  cross  section,  o(E).  It  is  well-known  that  (y,/)  reactions  for  incident  photon 
energies  below  the  threshold  for  photoneutron  production  resonantly  excite  discrete  states  like  that  shown  in  Fig.  1.  Each  resonance  is 
described  by  a  Lorentzian  line  shape  with  a  width  which  may  be  large  an  the  nuclear  scale,  but  is  narrow  compared  to  the  structure  of  most 
incident  photon  continua.  Furthermore,  if  gateways  are  well-spaced  compared  to  their  widths,  Eq.  2  becomes 
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in  which  the  amount  of  final-state  population  has  been  expressed  as  a  fractional  yield,  or  activation,  per  unit  photon  fluence,  A^fEg).  The 
(aT)j  is  the  cross  section  for  population  of  the  final  state  through  the  j*  gateway,  integrated  over  the  Lorentzian  line  shape,  and  the  summation 
includes  all  pump  bands  with  excitation  energies  Ej  <  Eq. 

Experimental  Confirmation 

The  model  of  Eq.  3  was  validated  by  experiments  using  bremsstxahlung  from  six  different  accelerators.  These  devices  included  two 
nuclear  simulator!,  DNA/PITHON  and  DNA/ Aurora,  two  medical  linacs  and  two  research  linacs,  the  injector  to  the  superconducting 
Darmstadt  linear  accelerator  (S-DALINAC)  and  the  Texas  X  at  the  Center  for  Quantum  Electronics.  Oram-sized  samples  containing  the  few 
calibration  isomers  whose  gateways  were  well-kncrwn  from  the  literature  were  exposed  for  times  ranging  from  seconds  to  hours  for  the  cw 
machines  and  to  single  flashes  from  the  pulsed-power  devices.  The  yields  of  isomers  produced  by  each  irradiation  were  obtained  from  the 
fluorescence  detected  after  transferring  the  samples  to  a  quieter  environment  Standard  corrections  were  mmt*  for  the  isotopic  abundance, 
isomer  lifetime,  detection  and  emission  efficiencies,  and  transparency  of  the  samples  to  the  fluorescence.  Photon  spectra  were  calculated  with 


the  EGS4  clcctron/pholon  Udiis|K>rt  code31  which  was  adapted  for  each  experiment  to  give  both  <X>0  and  F(E,E0).  hi  some  cases  <X>0  was 
verified  by  in-line  dosimetry 


Figure  1:  Schematic  representation  of  the  bandwidth  tunneling  of  Figure  2;  Fractional  activation.  A,  for  the  reaction  »7Sr(y.y),7Sr“  as 

population  to  a  laser  level  through  a  gateway  state  with  natural  ,  fUnction  0f  brexnsstrahlung  endpoint,  E0  obtained  with  5  different 

width,  r.  The  incident  photon  y  is  absorbed  with  a  cross  section  o0  accelerators.  The  solid  curve  plots  values  computed  from  Eq.  3 

from  the  initial  level  which  may  be  either  the  ground  state  or  an  ^  gateway  p^ct^  found  m  the  literature  and  calculated 

isomer.  The  cascade  from  gateway  to  laser  level  is  shown  by  the  ph(Jton  ^  agreement  0f  these  values  with  the 

dotted  line  with  promptly  emitted  photons  being  the  y.  Only  one  was  excellent,  validating  the  procedures, 

gateway  is  shown  but  there  could  be  more. 

The  results  of  those  experiments14"17-19  were  in  excellent  agreement  with  the  predictions  of  Eq.  3  using  literature  values76-77-29  for  gateway 
energies  and  integrated  cross  sections.  A  typical  example  is  given  in  Fig.  2  which  shows  a  composite  excitation  function  obtained  for  the 
calibration  isomer  ,7Sr™  using  five  of  the  accelerators.19-21  It  is  important  to  ndte  that  at  an  endpoint  of  4  MeV,  the  yield  of  the  ,7Srn>  isomer 
was  enhanced  16  orders-of-magnitude  by  the  process  shown  in  Fig.  1  over  that  which  could  have  been  produced  by  direct  excitation  of  the 
metastable  level  from  the  ground  state.  This  is  bandwidth  tunneling  at  the  nuclear  level. 

Those  experiments  established  that  there  could  no  longer  be  any  reasonable  doubt  of  the  model  or  the  procedures  for  quantitatively 
measuring  nuclear  fluorescence  efficiencies  using  bremsstrahlung.  In  addition  they  led  to  the  development  of  a  new  method37  for  providing 
absolute  characterizations  of  intense  x-ray  fields  using  isomeric  photoexcitation  reactions.  The  model  was  then  used  to  examine  candidate 
nuclei  for  a  gamma-ray  laser. 


Many  of  the  candidate  isotopes  for  a  gamma-ray  laser  are  deformed  nuclei  whose  energy  levels  are  built  upon  rotational  bands  like  those  of 
diatomic  molecules.  For  a  spheroidal  nucleus  an  additional  quantum,  K  must  be  conserved  which  is  the  projection  of  the  total  angular 
momentum,  J  on  the  major  axis  of  nuclear  symmetry.  Within  each  band  the  resulting  system  of  energy  levels  follows 

E„ ( K.J )  =  Em(K)  +  BmJ(J+ 1)  ,  (4) 

where  J>K>0,and  J  =  K,K  +  l.K+2, . . .  forK*0andJ  =  K,K  +  2,K  +  4,  . . .  forK  =  0.  The  B„  is  a  constant  and  En(K)  is  the  energy 
of  the  lowest-lying  level  in  the  band,  the  *bandhead.*  Other  quantum  numbers  are  identified  by  n.  The  selection  rules  for  electromagnetic 
transitions  require  |  AJ  |  <Mand  |  AK  |  <M  where  Mis  the  multipolarity  of  the  radiation. 

Deformed  nuclei  possess  the  longest-lived  isomeric  levels  because  those  states  lie  in  bands  whose  K  differs  greatly  from  that  of  freely- 
radiating  states  to  which  they  could  energetically  decay.  Such  isomers  would  be  ideal  as  storage  levels  in  an  upconversion  scheme  since  most 
of  the  energy  is  already  in  the  nucleus,  reducing  trigger  requirements.  Of  the  1886  distinguishable  nuclear  isotopes,  29  first-class  candidate 
nuclides  were  identified  by  a  search  of  the  nuclear  data  base  in  1986,  but  an  exact  ranking  could  not  be  determined  due  to  the  paucity  of 
nuclear  fluorescence  data. 


Gianl  Puinnine  Resonances 

Although  attractive  from  the  perspective  of  shelf  life  and  availability,  die  nuclide  1!<>Ta  was  nevertheless  ranked  among  die  poorest 
candidates  since  its  isomer  stores  only  75  keV  and  requires  AK  =  8  for  transitions  to  die  ground-state  band.26-”  A  bandwidth  funnel  would 
also  have  to  span  a  similar  AK  end  theoretical  estimates  of  the  integrated  cross  section  for  die  pumping  were  c\]Kcted  to  be  small,  jvih.-.ps  lu 
15  uu  Still,  this  isotope  exists  naturally  in  100  %  inversion  since  its  isomer  has  a  lifetime  in  excess  of  1012  veins  while  its  puue.J  state  has  a 
halflifc  against  transmutation  of  cnly  8  1  hours  A  macroscopic  sample  was  readily  available  so  18t>]'am  was  die  fust  candidate  tested 

An  experiment  was  conducted18  in  1987  in  which  1.2  mg  of  ,80Tam  was  irradiated  with  brcmsstrahlung  from  a  6  MeV  medical  linac  and  a 
large  fluorescence  yield  was  observed  following  the  exposure  That  was  the  first  experimental  demonstration  that  the  cncrgv  stoicd  in  a 
nuclear  isomer  could  be  dumped  to  freely-radiating  states  as  needed  for  a  gamma-ray  laser.  Simply  the  observation  of  any  fluorescence  from  a 
milligram-sized  target  when  grams  of  material  were  needed  for  the  calibration  isotopes  showed  an  unexpected  efficiency  loi  die  process 
Analysis  of  the  data  indicated  that  the  partial  width  for  the  dumping  of  the  l8t>Tam  population  was  about  0  5  cV  with  an  integrated  cross 
section  on  the  order  of  42,000  uu. 

The  energies  of  the  pump  levels  were  determined  from  the  excitation  function  of  Fig  3  measured  by  a  series  of  irradiations11,  using  die 
vanablc-cnergy  S-DALINAC.  As  was  the  case  for  the  calibration  isomer  ,7Srm,  sharp  increases  in  the  normalized  yield,  called  activation 
edges,  indicated  those  energies.  Computed  photon  spectra  were  used  with  Eq.  3  to  iteratively  fit  the  measurements  by  using  dial  values  for 
the  integrated  cross  sections.  The  lowest-lying  pump  band  was  located  at  2.8  MeV  and  with  an  integrated  cross  section  of  12,000  uu  it  was 
classified  as  a  "giant  pumping  resonance"  compared  to  usual  (y,Y)  reactions.  The  demonstration  that  populations  of  nuclear  isomers  can  be 
dumped  by  x  rays  with  much  larger  integrated  cross  sections  than  could  have  been  expected  stood  as  the  first  major  breakduough  towards  a 
gamma-ray  laser,  enhancing  its  feasibility  by  at  least  four  orders-of-magmtude. 

Svstematics 

Giant  pumping  resonances  were  also  observed20  for  the  second  poorest  candidate,  123Te.  However,  better  candidates  were  not  available  m 
sufficient  amounts  for  testing  so  research  focused  on  the  pumping  of  simulation  isomers  to  study  the  systematics  of  their  gateways  A 
preliminary  survey21  of  19  isomeric  nuclei  was  conducted  over  a  coarse  mesh  of  bremsstrahlung  endpoints  using  the  nuclear  simulators  and 
the  medical  linacs.  As  shown  in  Fig.  4,  giant  pumping  resonances  were  indicated  for  nuclei  in  a  region  of  masses  near  A  =  180,  where  the 
most  deformed  candidates  lie.  A  refined  survey22  performed  with  the  S-DALINAC  found  the  gateway  locations  and  integrated  cross  sections 
for  four  isotopes  neighboring  1S0Ta  that  had  various  degrees  of  ground-state  deformation.  The  results  are  also  given  in  Fig.  4  and  indicate  that 
both  the  excitation  energies  and  integrated  cross  sections  for  the  pump  levels  vary  slowly  amongst  neighboring  nuclides  within  each  mass 
island,  despite  dissimilar  single-particle  structures.  The  strengths  of  those  giant  pumping  resonances  were  correlated22  with  the  ground-state 
deformation,  a  nuclear  core  property. 
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Figure  3:  Fractional  activation,  Aj-  for  the  reaction 
l*°Tam(Y,y),,0Ta  as  a  function  of  bremsstrahlung  endpoint,  E„ 
obtained  with  the  S-DALINAC.  Gateway  parameters  were 
determined  by  fitting  the  measured  excitation  function  using 
Eq.  3  with  trial  values  for  the  (oT)j.  The  integrated  cross 
sections  were  found  to  be  104  burger  than  those  usually 
measured  for  (y,f)  reactions.  These  data  demonstrate  the 
dumping  of  an  isomeric  population  with  x  rays. 


Figure  4:  Systematics  for  excitation  energies  and  integrated  cross 
sections  of  gateway  states  providing  optical  pumping  of  nuclear 
isomers.  The  groupings  correspond  to  mass  islands  between  magic 
numbers  for  neutrons  and  protons.  The  best  candidates  for  a 
gamma-ray  laser  he  in  the  mass-180  island  where  the  largest  values 
of  (ol^  are  found.  Within  each  island  the  integrated  cross  sections 
and  excitation  energies,  plotted  by  the  right-hand  ordinate  and  the  X 
symbols,  vary  slowly  with  changing  mass  number,  A.  This  is  a  sign 
of  the  core  nature  of  the  giant  pumping  resonances. 


The  systcinatics  for  optically  pumping  nuclear  slates  were  further  investigated  for  the  excitation  of  laser-like  levels  with  lifetimes  on  the 
order  of  tens  of  ps.  To  accomplish  this  an  experimental  technique  was  developed53  which  employed  a  gated  CsF  scintillator  to  detect 
fluorescence  from  samples  irradiated  with  tlic  pulsed  Texas  X  Irnac  Tlic  PK1T  was  gated  off  during  the  x-ray  pulses  to  prevent  the  system 
from  being  blinded  by  real-time  Compton  scattering  from  tlic  sample  and  the  environment  A  total  of  lour  additional  nuclei  were  studied 
using  that  technique,  8lBr,  l76llf,  ,8lTa  and  20lllg  Figure  S  shows  a  typical  decay  spectrum  for  l8l'la  obtained  with  a  time-to  amplitude 
converter  (TAC)  and  was  collected  over  7  x  105  cycles  (linac  pulses)  Fluoicsccncc  was  detected  at  energies  of  133,  345,  476  and  4S2  keV 
from  tlic  decay  of  a  level  at  61 5  keV  and  the  agreement  between  the  measured  halfiifc  and  the  literature  value26  of  18  ps  was  excellent  The 
data  analysis33  for  all  four  nuclides  indicated  that  the  integrated  cross  sections  for  optically  pumping  laser-like  levels  were  comparable  to 
those  previously  found  for  the  population  of  long-lived  isomers,  as  shown  in  Fig  6. 
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Figure  5:  Fluorescence  detected  from  the  18  ps,  615  keV  level  in 
mTa  following  its  excitation  by  bremsstrahlung  from  the  Texas  X. 
These  data  were  collected  for  nearly  106  cycles  (linac  pulses)  using  a 
recently  developed  gated  CsF  detector  system.  The  PMT  was  gated 
ofT  during  each  x-ray  pulse  and  then  returned  rapidly  to  a  stable  gain 
to  observe  fluorescence  from  the  sample. 


in 


v*  0 


0 


Pumping 
laser -like 
levels 
IfUS  hfeliru 


50  100  150  200  250 

MASS  NUMBER.  A 


Figure  6:  Integrated  cross  sections  for  pumping  laser-like  levels 
along  with  those  of  Fig.  4  for  the  pumping  of  long-lived  isomers 
(shown  by  the  dotted  lines).  Again,  the  groupings  correspond  to 
mass  islands  between  magic  numbers  and  the  best  candidates  arc 
identified  as  neighbors  of  the  nuclide  U0Ta. 


NRF  Measurements 

The  discovery  of  giant  optical  pumping  resonances  significantly  unproved  the  feasibility  of  a  gamma-ray  laser  by  lowering  pump 
requirements.  However,  the  resolution  of  those  experiments  only  permitted  gateways  to  be  located  within  ±0.1  MeV.  This  left  open  the 
possibility  that  the  reaction  strength  might  not  be  concentrated  in  a  few  transitions,  but  might  be  spread  over  a  great  number  of  closely  spaced 
levels.  If  this  were  the  case  the  cross  section  for  each  participating  state  could  be  small,  seriously  degrading  the  ratio  of  useful,  resonant 
absorption  to  that  absorbed  nonresonantly  and  contributing  to  heating  of  the  sample.  This  was  of  particular  concern  for  the  nuclei  favored  in 
Figs.  4  and  6  which  have  large  excited-state  densities,  and  even  more  so  for  which  is  a  rare  odd-odd  nucleus.  A  series  of  experiments 
was  performed  to  explore  this  possibility. 

Isomeric  excitation  functions  like  that  shown  in  Figs.  2  and  3  reflect  only  one  specific  exit  channel  for  the  pumping  process  depicted  in 
Fig.  1.  That  channel  consists  of  a  decay  cascade  leading  to  the  metastable  level  following  an  absorption  event  in  a  transition  from  the  ground 
state  to  a  gateway.  However,  a  complimentary  exit  channel  corresponds  to  direct  decay  of  the  gateway  back  to  the  ground  state,  with  the 
prompt  emission  of ‘elastic*  photons  of  the  same  energy  as  those  absorbed.  It  can  be  seen  from  Fig.  1  that  in  the  event  of  a  large  (oT^  for  the 
population  of  an  isomer,  there  must  also  be  a  large  integrated  cross  section  for  the  elastic  channel,  (aT)eI  =  xb AvJTfl.  The  square  of  the 
branching  ratio  from  gateway  to  ground  state  occurs  because  the  entrance  and  exit  channels  are  the  same.  This  integrated  cross  section  can  be 
determined  by  measuring  the  energies  and  numbers  of  photons  scattered  in  real  time  during  an  irradiation.  The  observation  of  those  photons 
in  the  presence  of  Compton  background  forms  the  basis  of  nuclear  resonance  fluorescence  (NRF)  experiments. 

An  experimental  arrangement  designed  for  NRF  was  used  to  investigate23’24  the  pumping  of  two  isomers,  *9Ym  and  ,15Inm  which  have 
large  natural  isotopic  abundances.  Gram-sized  samples  were  sandwiched  between  aluminum  and  boron  calibration  targets  and  were 
irradiated  with  well-collimated  bremsstrahlung  from  the  S-DALINAC.  Photons  elastically  scattered  by  strong  transitions  were  detected  with 
carefully  shielded  germanium  detectors  placed  off-axis  and  appeared  in  pulse-height  spectra  as  peaks  superimposed  on  a  large  Compton 
continuum.  Isomeric  excitation  functions  were  also  measured  for  these  nuclides.  Detailed  spectra  are  available  in  the  literature23*34  but  the 
relevant  point  is  that  only  a  few  strong  elastic  transitions  were  detected  in  the  NRF  measurements  over  the  energy  range  from  1.5  to  4.5  MeV 
and  these  were  only  found  within  the  experimental  resolution  of  their  pumping  resonances.  In  the  case  of  n5In,  these  were  located  at  2.8  ± 
0.1  MeV  and  3.3  +  0.1  MeV,  as  shown  in  Fig.  7.  Thus  a  subset  of  those  few  transitions  corresponded  to  the  gateways  responsible  for  the 


optical  pumping  of  the  isomer.  That  conclusion  was  confirmed  qualitatively  by  umfied-inode!  calculations  which  showed  that  the  absorption 
strength  was  concentrated  into  a  small  number  of  El  or  E2  transitions  with  excitation  energies  near  those  found  experimentally 

K  Mixing  -  A  Speculation 

The  surprising  magnitudes  of  integrated  cross  sections  for  getaways  in  the  mass- 180  region  indicate  an  ease  in  pumping  and  dumping 
isomers  through  large  AK  that  is  difficult  to  interpret  in  terms  of  a  single-particle  model  Levels  capable  of  stiong  "K  breaking"  or  "K 
mixing"  have  also  been  recently  observed  for  other  nuclei54  in  which  electric  dipole  transitions  to  levels  near  2  5  MeV  were  found  to  be 
enhanced  despite  substantial  change  in  K.  Most  striking  was  the  measurement55  of  a  level  at  2  65  MeV  which  provides  AK  =  14  in  the 
spontaneous  decay  of  the  3.312  MeV,  3.7  ps  isomer  of  l74Hf  That  K-mixing  slate  is  remarkably  close  to  the  gateway  at  2  8  ±  0  1  MeV 
discovered  for  lire  dumping  of  180Tam  A  full  understanding  of  this  pervasive  behavior  is  not  yet  available,  but  there  is  a  suggestive 
speculation 

Calculations56  have  confirmed  that  some  nuclides  develop  multiple  minima  in  nuclear  potential  as  a  function  of  core  deformation.  Second 
minima  can  contain  fission  isomers  for  higher  masses,  but  for  lighter  nuclei  the  Coulomb  barrier  prevents  fission  Shape  isomcis  result  when 
states  bound  within  a  second  well  cannot  easily  decay  to  freely-radiating  slates  of  significantly  difTcrcnt  nuclear  slia|>c  corresponding  to  the 
primary  well.  If  instead  of  a  second  minimum  a  plateau  occurs  in  the  potential  energy,  as  depicted  in  Fig  8.  shape  isomers  could  not  occur 
However,  within  the  narrow  band  of  energies  shown  excited  states  would  not  correspond  to  a  fixed  core  deformation  Thus  K  would  not  be  a 
good  quantum  number  for  such  levels  and  unhindered,  low-multipolarity  transitions  could  connect  them  to  states  with  widely  differing  K 
values.  This  would  provide  the  large  degree  of  K-mixing  observed  experimentally.  It  is  interesting  to  note  that  shape  isomers,  and  therefore 
second  minima,  arc  predicted  to  lie  at  excitation  energies  of  2  to  3  MeV. 
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Figure  7:  Schematic  representing  integrated  cross  sections  for 
a)  population  of  the  isomer  nsInm  obtained  from  a  measured 
excitation  function  and  for  b)  elastic  scattering  by  absorption 
transitions  obtained  from  NRF  measurements.  Gateway  locations  of 
a)  were  determined  within  ±0.1  MeV  and  strong  elastic  scattering 
levels  were  only  observed  within  those  bands.  Thus  a  subset  of  the 
transitions  of  b)  are  responsible  for  the  population  of  the  isomer  as 
confirmed  by  unified  model  calculations. 


DEFORMA  TION 

Figure  8:  Schematic  representation  of  a  dependence  of  potential 
energy  upon  nuclear  deformation  that  would  tend  to  produce  a 
second  local  minimum  Shape  isomers  have  been  predicted  to  lie  in 
auxiliary  minima  and  for  large  enough  mass  numbers  spontaneous 
fission  can  occur  when  populations  tunnel  out  of  a  secondary  well. 
Here  a  plateau  rather  than  a  well-developed  local  minimum  is 
shown. 


The  First-Ranked  Candidate 

The  studies  of  the  systems  tics  for  optically  pumping  isomers  conducted  over  the  past  six  years  have  yielded  the  following  conclusions: 


1)  Nuclear  populations  trapped  in  long-lived  storage  isomers  can  be  optically  pumped  to  freely-radiating  states  with  x  rays, 

2)  Optical  pumping  reactions  for  the  deexdtation  or  excitation  of  nuclear  isomers  (and  laser-like  levels)  occur  resonantly  by 
bandwidth  fismeling  through  a  small  number  of  gateway  states, 

3)  Gateways  lying  in  the  mass-180  island  are  giant  pumping  resonances  and  can  break  K.  with  integrated  cross  sections  on 
the  order  of  10,000  uu.  compared  with  1  to  10  uu.  as  is  usually  found  for  (y,Y)  reactions,  and 

4)  Within  the  mass-180  island,  the  magnitudes  of  the  integrated  cross  sections  vary  slowly  among  neighboring  nuclides;  the 
excitation  energy  for  the  first  giant  pumping  resenanx  is  between  2.6  and  2.8  MeV. 

On  this  basis  the  31 -year  isomer  '^Hf®2  has  been  identified  as  the  first-ranked  candidate  for  a  gamma-ray  laser  from  the  perspective  of 
triggering.  As  can  be  seen  in  the  energy-level  diagram  of  Fig.  9,  the  metastable  level  stores  2.4S  MeV,  one  of  the  highest-lying  isomers 
known26  that  does  not  exceed  the  energy  of  a  K-breaking  gateway.  The  31-year  half  life  of  the  isomer  results  from  its  being  the  head  of  a  band 
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with  K  -  16  ll  must  first  decay  by  a  severely  hindered  E3  transition  to  the  J  =  1 3*  member  of  the  K  =  8  band  Ihis  AK  -  8  transition  is 
followed  by  a  cascade  that  ends  on  the  J  =  8  isomer  (l78lllTnl)  which  decays  with  a  halflifc  of4  0  s  to  the  ground-stale  Riii.l 
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Figure  9:  Energy  level  diagram  of  the  nuclide  178Hf.  The  31 -year 
isomer  stores  2.45  MeV  and  forms  the  head  of  a  band  with  K  =  16. 
The  natural  decay  occurs  through  the  K  =  8  and  K  =  0  bands  and 
cannot  excite  either  the  higher-lying  K  =  14  bandhead  with  a 
lifetime  of  68  ps  or  members  of  the  K  =  6  band  on  the  left-hand 
side  whose  bandhead  has  a  halflife  of  78  ns.  A  K-mixing  gateway 
level  inferred  from  systematics  is  shown  near  2.8  MeV  and  could 
populate  those  laser-like  states. 


The  hatched  level  in  the  figure  indicates  a  giant  ;niu,|  mg  re  .  -nance 
whose  position  is  infeiTcd  from  systematics  Tins  assignment  is  al 
supported  by  recent  Coulomb  scattering  mea  uieineiit.'  winch 
identified  a  level  at  an  excitation  energy  of  353  I.eV  above  the  isomer 
and  which  could  not  exclude  other  levels  lying  even  closer  The 
gateway  could  be  reached  by  the  absorption  of  a  pump  photon  of  less 
than  “iAn  kcV,  shown  by  the  large  upward  anew  in  the  figure  Due  to 
the  >.2  dependence  for  the  pump  step  shown  in  I:q  I,  this  reduced 
photon  energy  may  enhance  the  integrated  cross  section  by  a  factor  of 
nearly  SO  over  that  observed  in  18l>ra  The  presence  of  two  possible 
output  channels,  one  through  the  K  =  14  band  and  the  oilier  through  the 
K  =  6  band,  greatly  increases  the  probability  of  success  Pragma’icalK, 
the  lifetimes  of  both  laser-like  levels  arc  well-suited  to  the  tspical 
output  durations  of  nanoseconds  to  microseconds  of  pulscd-powcr 
devices  likely  for  use  m  optically  pumping  a  gamma-ray  laser 

Decay  from  the  gateway  to  the  K  =  14  band  would  require  only  AK 
=  2  and  would  populate  the  bandhead  having  a  halflife  of  68  ps 
Observation  of  this  event  followong  the  irradiation  of  a  sample  would 
require  the  detection  of  140  keV  photons,  as  shown  in  the  figure.  This 
energy,  as  well  as  the  68  ps  halflife,  is  already  well  within  the  present 
capabilities  of  the  CsF  scintillator  system 

In  the  event  of  a  stronger  degree  of  K  mixing,  the  K  =  6  bandhead 
with  a  halflife  of  78  ns  could  be  reached.  Requiring  AK  =  8,  this  is  no 
larger  than  that  transferred  through  the  gateway  in  the  dumping  of 
I80ratn  and  iess  than  that  observed  in  the  decay  of  ,74Hfm.  It  is 
expected  that  the  next  generation  gated  detection  system  will  allow 
measurements  on  a  time  scale  approaching  the  78  ns  halflife.  In 
addition,  the  fluorescence  signatures  are  at  the  relatively  high  energies 
of  922  keV  and  1247  keV.  Those  are  particularly  convenient  energies 
for  detection  if  bremsstrahlung  with  a  lower  endpoint  is  used  to  pump 
the  sample. 


Unfortunately,  from  the  viewpoint  of  availability  is  far  inferior  to  18t7Tam.  While  180Tam  exists  at  0.012  %  of  all  tantalum  and  can 

be  prepared  simply  by  mass  separation,  the  31 -year  isomer  of  178Hf  must  be  prepared  by  nuclear  reactions  such  as  spallation  cr  a 
bombardment  At  present  the  entire  world  inventory  ofnHfm2  appears  to  be  about  1015  nuclei.  Procurement  of  a  sample  sufficient  for  testing 
remains  a  problem  being  addressed. 


Conclusions 

The  discovery  of  giant  K-breaking  resonances  fot  pumping  and  dumping  nuclear  populations  and  the  development  of  new  technology  have 
enhanced  the  feasibility  of  a  gamma-ray  laser  by  orders-of-magnitude  over  that  projected  in  the  blueprint  of  1982.  All  of  the  component 
concepts  have  now  been  demonstrated  experimentally,  and  a  recent  model31’39  supported  by  the  new  data  indicated  that  there  is  a  significant 
region  in  parameter  space  for  which  a  nuclear  system  can  be  optically  pumped  to  threshold  without  the  need  to  melt  a  low-Z  host  lattice.  In 
addition  simulation  hosts  have  been  created  by  implanting  isotopically  pure  nuclei  of  37Fe  into  a  thin  diamond  film  produced  using  well- 
established  technology.40  At  this  time  here  are  no  a  priori  obstacles  to  the  realisation  of  a  gamma-ray  laser.  However,  exact  pump 
requirements  depend  sensitively  on  the  structure  of  the  actual  isotope  and  the  feasibility  of  a  gamma-ray  laser  cannot  yet  be  determined. 
Research  by  the  United  States'  main-line  project  has  now  identified  the  3 1-year  isomer  of  17*Hf  as  the  first-ranked  candidate  and  procurement 
is  currently  being  pursued.  Testing  is  planned  for  the  coming  year. 
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